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HAUPTAUFSÄTZE 


Die Drehströmung über festem Grunde. 
Von U. T. Bödewadt in Göttingen 

$ 1. Die Aufgabe. Von den Bewegungen einer sich drehenden Flüssigkeit über einer 
ebenen Platte war bisher der Fall noch nicht mathematisch bearbeitet worden, in dem es sich 
um den stationären Endzustand einer im Unendlichen mit gleichmäßiger Winkelgeschwindigkeit 
kreisenden Flüssigkeit über rullender Grundebene handelte. Den stationären Endzustand einer 
im Unendlichen ruhenden Flüssigkeit über drehender Ebene behandelte W. G. Cochran'); 
die Anfänge der Bewegung bei plötzlichem Einsetzen der Plattendrehung relativ zur Flüssigkeit 
hat in beiden Fällen K. H. Thiriot untersucht, der auch alle vier Fälle nachgemessen hat’?). 
Hier soll der zuerst erwähnte Fall (angenähert im Kerngebiet eines Zyklons verwirklicht) 
unter Anlehnung an die Methode und die Bezeichnungen bei Cochran behandelt werden. 

Die nicht zusammendrückbare Flüssigkeit habe die Dichte o und die Zähigkeits- 
konstante »; sie stehe unendlich hoch über der (x,9)-Ebene z2=0. Bei Einführung von 
Zylinderkoordinaten r, g,2: 2=rcosp, y=rsinp, 2=z Seien u, v, w die radiale, tangen- 
tiale und vertikale Komponente der Strömungsgeschwindigkeit d; die Navier-Stokesschen 
Gleichungen 


Bgrad)p=— grad p+r: Av; Ge. ul a 
gehen dann über in 


wur + Zug D) 5 ev” +nwmu,=— Ip r+ »lurr+ n + upp+ U) 0 
v, + I o+luntmw = — „pet vlorr in ır + vop-+ 0.) u A, 
+ ngpo +nw,= — SP: +rlwrrt n “+ Im wog +.) or 
tu +04 EN ET T 


Die Strömung soll nun drehsymmetrisch um die Achse r=0 sein, d. h. die Ableitungen 
nach x verschwinden. Setzt man 


u=rfe, verg@, w=hl) ...:......(13), 


1) The flow due to a rotating disc. Cambridge Proc. Bd. 30 (1934), S. 365 bis 275. 
2) K.-H. Thiriot: Über die laminare Anlaufströmung einer Flüssigkeit über einem rotierenden Boden bei 
plötzlicher Änderung des Drehungszustandes. Z. angew. Math. Mech. Bd. 20 (1940), S. 1 bis 13. 
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so vereinfachen die Gl. (1.2) sich zu 


Pt ht len nn. . (14.1), 
2fg+hg: = ee. 
+an=—Ip, a 
2f+h, = 0 (1.4.4). 


Im Grenzfall verschwindender Reibung (»=0) soll eine reine Drehung mit fester 
Winkelgeschwindigkeit ® um die 2-Achse (r=0) vorliegen: f=0, g=w, h=0; dann folgt 


oa’ A P,. Für sehr kleine Zähigkeiten » nimmt man nach der Prandtlschen Grenz- 


schichttheorie diesen Wert der horizontalen Ableitung des Druckes als überall gültig an und 
erhält so unter Fortlassen der Gl. (1.4.3) das System 


PIE EFF : 8 nn. 2 rn si 
2fg+hg:= EEE > ° 
a ae 

Um dimensionsfreie Größen zu erhalten, setzt man 
2—= Vz: f@)=v»-F(Ü), ge)=w-G(d), h(2)=yYwvH(£) . ..(1.6). 

So entsteht schließlich 

N „en 
2FG+H@G= a 
2F+H = 0 Be a 


Randbedingungen sind: Haften der Flüssigkeit an der ruhenden Grundfläche (©=0) und 
gleichmäßige Drehung in großem Abstande von der Grundebene ( =»): 


FO)=0, 60)=0, HO=0 ........ (81, 
BER: We ee ee 


Cochran hat das bei ihm auftretende, dem vorliegenden ’°) (1.7), (1.8) ähnliche System 
dadurch gelöst, daß er einerseits Potenzreihen, andererseits auch eine asymptotische Ent- 
wicklung benutzte und beide bei einem endlichen Wert von £ aufeinander abglich. Bei der 
Übertragung dieses Verfahrens bestand die Hauptschwierigkeit darin, die asymptotische Ent- 
wicklung aufzufinden. Diese zeigte sich dann aber für das ganze Gebiet OS{<x so stark 
konvergent, daß für eine grobe Lösung die Potenzreihen noch nicht benötigt wurden. Nur 
weil die Konvergenz der asymptotischen Reihen nicht ausreichte, um mit den Gliedern bis 
zur 6. Ordnung einschließlich (mehr zu nehmen, verbot sich aus Zeitgründen infolge ihres 
stark zusammengesetzten Aufbaues) den Fehler unter 0,03 herabzudrücken, wurden zuletzt doch 
noch Potenzreihen verwendet und bei {=1 mit den asymptotischen Reihen abgeglichen. 

Es sei bemerkt, daß diese Methode (einschließlich der zur Auffindung der asymptotischen 
Reihen) sich ohne weiteres auf alle Gleichungen vom Typus 


y'"=F(y,y) (F ein Ausdruck 2. Grades mit konstanten Koeffizienten) 


übertragen läßt, ferner auf entsprechende Gleichungen mit Ausdrücken höheren Grades und 
auf Systeme solcher Differentialgleichungen, wenn die Randbedingungen von ähnlicher Art 
wie (1.8) sind — vorausgesetzt, daß die gefundenen Reihen genügend stark konvergieren. 


$ 2. Der Lösungsansatz. Fülrt man in die Gl. (1.7) für die Funktionen F, G, H die 
Potenzreihen 


& [” 
FO=a +5 + t+mnt 

= “ 
BOCH + he 
H()= rt 








3) Vgl. die G]. (16) bis (18) bei Thiriot, a. a. 0.8.6. 
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ein, so ergeben sich Rekursionsformeln, aus denen für die ersten Koeffizienten die in Zahlen- 
tafel 1 zusammengestellten Werte folgen. Darin bleiben die beiden Konstanten A, B noch 


Zahlentafel 1. Die Koeffizienten der Potenzreihen. 


a=—A db, =+B eu =d 
a,=+1 b,=® .=+24 
4 =0) b, =D g=-2 
a =—2B? b, =—-24B 4 =—0 
Pe y b, =+4B os —=+4B: 
.=—2 u =® GG =—44 
a, =—4AB? b, = — WU 4?B — 16 B® co —=+4 
ds = + 40 4? — 4 B? bs =+112 AB Ce =+84AB: 
a = — 144 b, == — 136 B Cs = —804?+8B? 
Aa =-+ 144 — 128 A? B? +88 B* bu = — 80 4? B— 484 AB? Co=-+ 288 4 
bu=+ 8272 4?B-+ 1472 B® 
bi = — 27520 AB 


unbestimmt, durch deren passende Wahl die Grenzbedingungen (1.82) zu erfüllen wären. 
Doch werden die Reihen (2.1) für große Werte von £ unbrauchbar. Es muß daher eine bei 
=» konvergente Entwicklung gefunden werden. 

Aus (1.7.3) und (1.8.2) folgt H’(x)=0. Für H(£) ist daher ebenfalls ein fester Grenz- 
wert zu erwarten, den es noch zu bestimmen gilt: 


A A a 
Damit wird (1.7) bei &=» befriedigt. : 
Für F, @, H werden jetzt Reihen angesetzt aus noch unbestimmten Funktionen, die 
bei ©&=» verschwinden müssen; p ist eine frei wählbare Konstante: 


FO=-0+r,. OHM O+--, 
e=l+rv, tw t+PWwÜ)+t:-, RE 
H)=h+p. +. tr O+t 
Diese Reihen sollen für die folgenden Umbildungen soweit wie notwendig als absolut und 
gleichmäßig konvergent vorausgesetzt werden; ob diese Voraussetzung zutrifft, wird sich nach- 
träglich zeigen '). 
Setzt man die Reihen (2.3) in die Gl. (1.7) ein und faßt gleiche Potenzen von p zu- 
sammen, so wird das erhaltene Gleichungssystem sicher dann befriedigt, wenn die Koeffizienten 


jeder Potenz von p einzeln verschwinden. So erhält man zur Bestimmung der unbekannten 
Funktionen die Gleichungen 


BREI EBEN ati er 

yv—h pr -29,= TE ee a ee ee A ni (2.4.2); 
k—1 

PR" how +2 yn= 2 (WnPR-a— YaYyR-at Pa Zum en ARE, 
k—1 

’ ’ < WU ’ - 

vr —hyr -2pr—= I @pıyr-a +vin-) : -- : @88; 
4=1 

Be a 2 re sr Zn or 


Der Wert von p spielt dabei keine Rolle, sondern erst bei der Erfüllung der Anfangsbedin- 
gungen (1.8.1). Dann kann man aber den Faktor p* auch in x, Yx, X» aufnehmen oder, 
was auf dasselbe hinauskommt, jetzt 


ER ey N 
setzen. Dadurch gehen die Reihen (2.3) über in 
FO=-0+P,I+R.I+MI+:-, | 
Eed=el+v, ++ O+-, ee 
H9)=h+,.0+. ++: 


4) Da alle diese Funktionen denselben linearen Differentialgleichungen mit konstanten Koeffizienten (2.5), (2.6), 
nur mit verschiedenen rechten Seiten gehorchen, so folgt ihre Ab- bzw. Zunahme mit der Ordnung k einem exponentiellen 
Gesetz, das noch von Schwingungen überlagert sein kann. Die numerisch augenscheinliche Konvergenz der ersten 
Glieder (Zahlentafel 2) genügt daher, um die Konvergenz im ganzen zu sichern. 
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Um sofort mit diesem Ansatz zu beginnen, hätte man bei der Produktbildung die Cauchy- 
sche Reihenmultiplikation anwenden müssen, doch würden die Gl. (2.4) bis (2.6) dabei sehr 
viel willkürlicher erscheinen. 


$ 3. Die erste Lösungsstufe: Die Funktionen p,(d), y,(d), x.(£). Setzt man für das 
System (2.4) eine Grundlösung in der Gestalt 


Yı ()=e":, Y; (= (er: 


an, so wird der Exponent « komplex: u=—a-+iPß; in reeller Schreibweise hat man 
Be I ar ER, 
RN hd + BREMER SEERGEENEEEEEEE EZ SE - = 2 
Damit wird (2.4) zu 
F"—(2a+h)F’+ala+h)F,+26,=0. ..... . (821), 
G"—Za+h)@, +ala+hh@a,—2F,=0. ..... . (82.2), 
und die Funktionen F,, @, haben die Form 
P Alab Lu 9 9 7.9 VW 7; © | 7 
O0=4,08BT+°8,sinßl. . » » 2: 2 2.20. (832). 


Zur Bestimmung der Konstanten FC,, FS,, *C,, 8, erhält man durch Einführen von (3.3) in 
(3.2) vier homogene Gleichungen. Ihre Determinante muß verschwinden; sie hat den Wert 


(aa +) — BP +4—(2a+n? "+4 [a(a+n) — B?|@a+n) pP. 


Da die beiden Summanden Quadrate reeller Zahlen sind, müssen sie einzeln verschwinden. 
Der letzte Faktor des zweiten Summanden kann nicht Null sein, da dann der erste Summand 
positiv wäre; also muß der vorhergehende Faktor im zweiten Summanden verschwinden: 


a ae ae nie 
Damit der erste Summand ebenfalls zu Null wird, muß außerdem 
EEE San en en ER 


gelten, wobei das Vorzeichen der linken Seite noch unbekannt ist. 
Aus (3.4) folgt 
h as h\? 
Verl) 
Nun sollten die Funktionen , (2), y, (£) für {>x verschwinden. Darum muß «a positiv sein: 


a=- + #+(b) Be u Se 


Nach (3.4) und (3,5) sind a und a-+h die Wurzeln einer der beiden quadratischen 


Gleichungen 
2’ + 7 2+P=0, 


deren Lösungen sich als 


zi:V! » 7 
+71) Wr Gehirn nie A 


angeben lassen. Da die beiden Lösungen a und a-+h reell sind, so muß der Radikand 
in (3.7) positiv sein, also #°?<1. Ferner kann  0.B.d, A. positiv genommen werden, dem- 


nach ist 
0<PB<1. 


Der halbe Unterschiedsbetrag der Wurzeln « und «+h kann auch aus (3.7) entnommen 


werden: 
'h 1 
VRR 


Quadriert man diese Gleichung, so läßt sich daraus £ als Funktion von %h finden: 


Ba er > h 
P=+V—ga+yi+Q mit ı=8 a Er EEE . (3.8). 








XUM 
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Aus (3.6) und (3.8) folgt die Formel 


1 h 
a Bk ı Br De SUN EEE 
man kann wegen (3.8) « auch nur durch oder nur h darstellen. — Aus (3.9) ist das in (3.5) 
noch unbestimmt gebliebene Vorzeichen zu erkennen: 
BE le a dr 


Bei den weiteren Rechnungen konnte manche Vereinfachung durch die aus (3.4) und (3.10) 
fließende Beziehung 


6) 
m te Ar Bu Sahara RE SE 7 © 
erhalten werden. 


Die Erfüllung des homogenen Gleichungssystems für die in (3.3) auftretenden Kon- 
stanten erfordert nur noch die zwei Gleichungen 


Fr =—#S8 78,=0C,. 
Darum sollen hierfür zwei neue Bezeichnungen f, g eingeführt werden: 
4, S=+g, 
| ei ir Tr ee 
c=+9, Ss=—f 
Aus diesen drei durch (2.2), (3.12) erklärten Konstanten f, g, h lassen sich alle Konstanten 
der weiteren Lösungsstufen berechnen. 
Hinsichtlich der Funktion z, (£) folgt aus (2.6) für k=1, daß sie die Form 
DE a ern e 2 
DEE BB. eier a 
haben muß; die Bestimmung der Konstanten soll unten im Zusammenhange mit denen der 
höheren Stufen erfolgen. 
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$ 4. Die weiteren Lösungsstufen. Für k—=2 lauten die Gl. (2.4): 
% hoyy+2y,=P7?—- vw th: -» :::.. 2. AL), 
vw hy, 20, =23o, m tYV/ı - :» : : 2.0.0.0. 12). 
Die rechten Seiten sind Produkte aus dem Faktor e-?«° und trigonometrischen Polynomen 
der Frequenzen 28 und 0. Es gibt also ein partikuläres Integral von (4.1) der Gestalt 
O=eEF(O0; wO=e-?.:0,(), 
mit 
F,(J=F(C} +FC}cos2ßTC+FS}sin2 PL, 
G,()=0} +05 c0os2 BC +°S3sin2 PL. 
In dem allgemeinen Integral ist außerdem ein Vielfaches der Lösung des zugehörigen homo- 
genen Systems enthalten. Dieser Anteil werde mit der ersten Stufe vereinigt, hier also (und 
ebenso für k=3,4,5,...) null gesetzt. 
Die Fortführung dieser Überlegungen ergibt, daß allgemein angesetzt werden kann 


pr(d)=ert*iFr(d), 
ET an. > 7 FUHREN — © 
1(d)=e-ke!H,(d) | 


wobei für k=0 3 
F,(O)=0, Gm, H,(Ö)=h 


ist und für k=1,2,.. 
Fr (O= I ("C),cosaßc+ "S}sinAßL), 


OL Ceosaßc+"ssinß), 0... 0.8); 


4. 





Hr, (= 2 ("C),cosA BL +" 87, sinA ß?}) 


4 


der Summenzeiger A hat darin die Werte k, k—2, k—4,...>0 zu durchlaufen. Für A=0 
ist stets 
ee a ae ee : ° 
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Die Gl. (2.5), (2.6) werden durch (4.2) zu 
F," —(2ka+h)Fy +kalka+h)F, +26; 


Sry ne 5 je ee 
= {F, Fe_n— On Gk-nt+(Fn —_na F,) Hx-n} | 


n=1 


Gr — (2ka+h) Gy +kalka+h) Gk —2F; 


B55 Fa i ; a ne a: 0er SE 
u ” (2 F„ @k-n + (an — na G,) Hr_n} 

n—=1 
Hy kaH,+2F,=0 LE ER BT re re Fe (4.6). 


Aus (4.6) entsteht durch (4.3) ein Gtichungeau, das ne. und AS: liefert : 


Ad Betr 5 lkac, +Ap”S)). 
; 9 . (4.7); 
re 
für A=0 kommt einfach 
ua. a Er Ta nee 3 A 


Die Koeffizienten von H,(£) in (3.14) ergeben sich aus (4.7) für k=1, A=1 wegen (3.4), 
(3.11), (3.12) zu 
AC,=ßPf+(a+hlg, | 


& . (4.9). 
"S,—-ßg-(atnf | es 

Im Fall A=k>1 verschwinden übrigens alle diese Koeffizienten: 

Fok—FSk_6 ck _6gk Fk _FSk_ 0 für k>2, 
wie unten noch gezeigt wird. 
Die rechten Seiten der Gl. (4.5) nehmen bei Einsetzen von (4.3) die Formen 

= I YCgeosaß? TEmHB.....:... daR 
= IV Op eosAßt+"S,sinAßl} Be 2 


an; A durchläuft wieder die Werte (k), k—2, k—4,...0. Die Berechnung der Koeffizienten 
Y GG: (SL, n, %S7 geschieht im folgenden Paragraphen. Links treten dieselben Frequenzen 
auf; die Gl. (4.5) spalten sich dadurch in mehrere Systeme von je vier linearen Gleichungen 
folgender Art: 





+24’7G-2B’S, +2 °C le, 
+2B,"C, +24," 8, +2 ='S}, 
u 7 +2 4°, —2B'S=°ch, arg 
2 "+23 +2 45-95} 
Hierin sind Ah, B, PR für Ausdrücke in a, ß: 
i 2 d=kalkath)— ?R—=k(k—1)a +(k— 42) ß®, 
a an .. 
Die Auflösung von (4.11) geschieht nach dem Schema 
FC -I4+ M'/CH-N!S—P’C—Q°S), 
FS=—(—-N/C+M!S+Q°’C—P°S), 
(4.13), 


G öde 


-(+PfC+OIS+M’C+N*S), 


og 


(- @/C+P/!S— N’C+M’S) 
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in dem M, N, P, @, 4 die Bedeutung 


M=A(#?+B-+1), N=B(4?’+B’—1), 
P=4-—B’+1, Q=2AB, ge 
=2(P’+ Q°) 
haben; 84 ist die Determinante des Systems (4.11). Allen Größen in (4.13), (4,14) sind je- 


weils dieselben zwei Zeiger k, A zu geben. Für A=0 verschwinden die S-Koeffizienten, 
es bleiben 


FW 1 fr gpd 
aus Fa nl 
1 (4.15) 
imo gr od 
"= p,En .+'C,) 
mit 
k 
BE 2 SE 2 
Er = (a? (k +): | RE" ER TE 


D;=2(1+ Ey?) 


$ 5.- Die rechten Seiten von (4.5). Um die Koeffizienten /O, (S,’C,’S der rechten Seiten 
von (4.5) zu finden, muß man zunächst Produkte der Form F, @G;_„ entwickeln. Aus (4.3) 
gewinnt man dafür die Formel 
u 


a) - 1 
F, G-n= B les» BL: 5 Ye 


% 


> Zee 
run 
vi (€ N 2% Gc)\ 


+sinaßL- 7 | 


u 


er. a E Fo: 
+ j \eos2 BE: „ rin +"s°8) 


n u) 
4b. v 
u rr u | Bm Ian 
+sini ßL- > Zrc 0 —- iR c)\, 

in der x die Werte der Reihe k, k —2, k —4, >0 durchläuft («=0 nur bei geradem n), 
a die zwischen n — 2x und x liegenden Werte der - Reihe n, n — 2, n—4,.. 20; A nimmt 
die zwischen k— n und —n befindlichen Zahlen der Reihe k, k — 2, k—4,... an, und für v 
kommen die zwischen k— n—4 und —/4 liegenden der Zahlen n, n— 2, n—4, ...20 in 


Betracht. (Soweit in der Summe dieselbe Zeigeranordnung sich wiederholte, ist sie nur beim 
ersten Gliede hingeschrieben.) 
k=]1 
Man bilde jetzt den Ausdruck — F,„@xk_-n, vertausche die Summationsreihenfolge und 


a1 
wende das Ergebnis auf (4.5) an. Nach Änderung einiger Bezeichnungen erhält man die Formel 


fer er 2 ©2 (’ ° BB —"0-°c+°8 ( 


o—( 


+lo B"S— u aFe)-"C+loßC+u a "ss 





k—4 
+3 > a Aa ER A Pi 
—— 
2 I v 
PS RR 58 Fa Fr 
= (7 BFs va C) C+ C-(rß S—vra‘C) 


- (B’C+ra Fs) . "S-"S-(rß’C+va”S)! 
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und drei ähnliche für /S;, ’C;, ”Sx. In diesen Formeln (k=Z22, 0O</sk—2) durchläuft « 
die Werte aus der Reihe 0, o+2,0+4, ..., o+k-—4, für de 1<usk=1 gilt, und v 
die Werte 7, +2, r+4,...., für welche die Ungleichung 1<£r<sk — A — r zutrifft. 


pri . .. .. F ge 
Für A=0 entstehen entsprechende, doch kürzere Formeln für ’C, und °C;. 


Bei der Auswertung ist auf die Gl. (3.12), (4.4), (4.9) und (4.10) zu achten. Die letzte 
ist so zu beweisen, daß man zunächst die den obigen entsprechenden allgemeinen Formeln für 


'%&,... aufstellt und dann durch Nachrechnen bestätigt, daß Ü=...='S,=0 ist wegen 
(4.9). Durch (4.13), (4.7) folgt daraus n.. ==. Die allgemeinen Formeln für 
'C, ... enthalten für k=3 in jedem Gliede mindestens einen Faktor dieser Reihe . al 
so daß auch die ’C,, ... verschwinden, usf. 


Um die Funktionen F(d), G@(&), H(£) nach (2.7), (4.2), (4.3) bis k<S6 aus f,g, h zu 
berechnen, braucht man zunächst folgende Hilfsgrößen: 


S=f’+g; te ae 5 

v—= 5 (2 3 h ß):; w=— n,(«+2m; u=wa. 2 Me 
55 1 

r=zlf-Pgp; d=z(ag+Pf) (5.3), 

P=af+Pßg; a A SER 

,=—ßPB; u, = — BY f 

=, —(k- Bf; K=t4,—(k—3)g ER vu. 


Aus diesen lassen sich schnell 


H S; H t ” 
e=—- Ss=-- (5.6) 
a ı a 
und 
Fr" be» Gp® H,;"! 
C‚=uS; 0,‚=vS; C‚=wS (5.7) 


finden, mit denen weiter folgende Hilfsgrößen zu bilden sind: 
„= — 2°C); 9=(hB—1)*c}; ir 2 ee | 


Für k=3 ist die Anwendung des Schemas (4.13) nicht mehr zu umgehen. Man erhält nach 


möglichster Vereinfachung für die /C, /S, ’?C, ?S bis k=6 die nachstehenden Formeln: 


E=+99n+®9; = —-,n+ Pc; | in 
. 2.J3). 


ER BRNO a war Lg 


(Hier wie im folgenden sind zur besseren Übersicht die ebenso wie in der ersten Zeile wieder- 
kehrenden Faktoren in der zweiten und den weiteren Zeilen fortgelassen.) 


2 5} 


< Y < "a! rn FC y o- in 
« ge + 28, ’C, +2 1,"8, I C+f S, Y °C # "Ss, ö °C Pr c | (ö 10) 
o). 
=+g ft Be ee + 
drin 
Szı =+k + 8% +f 9 — N —_y ; 
°C71 =+g +f + 5x — ba ö —y s . &.11). 
S=-f +9 3 m +r ii 
(gültig für k=3,4,5,...) 
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= +3, +3,79 + Sr" —-5"S); 








S=—3t, +, —f g Pe 
IG=+g —f +3s, +3t, ö +y 
S,=+f +9 De 
(= —29 CA+2Y-ı19”’C; 
IS,=+2f +2(— 1%) 
(5.12). 
=+(4s,+2f)) +2(+ 9 
Seth +2d) —-20+ 9 ; 
= + lc+lc)-2a+c”Cc,+ n+fSl.e; 
S=lir+rls (rt) + 5 -FC, 
eirHil- U Head) +°8 
1S,= IS/= IS} — es, +(@s)-a-”’8) —°C, 
Gert N EHE HE 
Par ee A 2. zu; Haie... Bat, An A 
fe? _ fe? Fe F7® Gp® P? (ra? | Fat? 6.18); 
G= /d-3’c”ch Fr RER EM 
= 10, +°C 2°C, -BI-"C,°8, + 8,0) 
. . Hy! Hs 9a? = p? 
y bedeutet hier den bei der Berechnung von "(,, , schon aufgetretenen Ausdruck . 
Cr, rag Cr‘ „a-s"s,; 
Beh, + — g 2 au 
IG=+9 -f +4, + —5 +7 
St 0 +0 0 4, +, —y RE 
d=+2,’7%,-2,"5—9 °c —f Ss, Pc +6#® 
IS= +21, en +f g ö ER 
Id=+g 57 +28, DE ö y 
Sul EEE. d BEREE, 3 U SEBBE 2 7 SEE > Zen Zur 
(5.14). 
Get °G -n+ S-28; 
Selten) + 5 -'G 
(= 0, .g® c} N Bi +%E 
g =—2°5, ® "Ss, Ss?) °C ; 
1 2 DE, De Be 
lH rer +8‘); 


'S za Is, ” is, + Is, —"C, ’S, Sr “, °S;; 
=: +30, +10; 
SE + is, 
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$ 6. Lösung der Anfangsbedingungen durch die asymptotischen Reihen. Aus (2.7), (4.2) 
und (4.3) ergeben sich die Anfangswerte bei {=0 als Summen der C-Koeffizienten: 


FO)=0+ f +’ÜU+ "Hl" +H)+---, | 
+ rar +. (6.1) 
HW=r+®c +++ +Fc)+--- | 

Für die Anfangswerte der Ableitungen von F(£) und @(£) bekommt man 

+a@’d+s8’ ++ "Ü)+--) 

pe "5, +2"5, + +38) +); | 
B=+6'(0)= - Bi +3 ++ +: ) | 6.3) 
+ Bl, +) 


Die Abschnitte der Reihen (6.1), (6.2) bis zu den Gliedern der k-ten Stufe einschließlich sollen 
durch Fir), @k), ---, Bim bezeichnet werden. — Auch für den Quotienten A/B wird eine 
Bezeichnung gebraucht: 


A=—-F(O= 2y 


4=m; entsprechend: Be m) IR en ne a 

Die Werte (6.1) sind nach $ 5 Funktionen der drei Parameter f, g,h. Diese sind nach 
(1.8.1) so zu wählen, daß die drei Funktionen verschwinden, d. h. das Gleichungssystem (1.8.1) 
ist nach f, 9, k aufzulösen, was hier durch fortschreitende Näherungen geschehen soll. 

Zu diesem Zwecke wird zunächst das Gleichungssystem F„=@,=H„=0 ge- 
löst. Ausgangswerte dafür erhält man aus dem System F„=f=0, G»=1+9=0, 
H„=hll+M+26=0. Aus den beiden ersten folgt f=0, 9g= —1; das dritte fordert dann 
a=0. Damit « Null oder doch möglichst klein wird, muß nach (3.8), (3.9) h>0 sein; zu- 
nächst genügt h—=2. Mit diesen Ausgangswerten wird das System F,,=:-':=0 nach dem 
Newtonschen Verfahren gelöst; für die dabei benötigten Ableitungen der drei Funktionen 
erhält man brauchbare analytische Ausdrücke. Nach einigen Näherungsschritten erhält man mit 


f=-031; 9=—0%8; h=+121 
die Werte 
F»=+003; G„=-—001239; Hu. = + 0,0008. 

Die daraus folgenden Verbesserungen von f, g,h liegen in der Größenordnung 0,01, so daß 
nun die höheren Ordnungen berücksichtigt werden müssen. 

Mit den zuletzt verwendeten Werten von f,g,h kommt z. Be H „= — 0,0815, 
H „= — 0,0%0. Die hinzugekommenen Glieder sind also keinesfalls zu vernachlässigen. 
Formeln für die Ableitungen von F,,,... aufzustellen, lohnt nicht bei der zu erwartenden 
Unlandlichkeit. Versucht man ersatzweise die formelmäßig bekannten Ableitungen von 
Fa, ... zu verwenden, so zeigen sich nur langsame Fortschritte. Es empfiehlt sich daher 
-— und zwar sogleich für F\„,,-..—, Durchschnittswerte der Ableitungen numerisch zu be- 
stimmen, indem man f, 9, h nen um passende kleine Beträge ändert und die dadurch er- 
zeugten Änderungen der Funktionswerte feststellt. Die so einmal ermittelten Werte für die 
Ableitungen werden dann bei den weiteren Näherungen unverändert weiter verwandt. 

Bestimmt man so einmal zu den Werten 


f= - 0316; g= 10106; h=-+ 1,3252 


die Funktionswerte F,,, ..-, M,.„, (auf der vollautomatischen Mercedes-Euklid-Rechenmaschine 
waren dazu etwa 10 Stunden erforderlich), so zeigt eine Betrachtung der Weıtfolgen 
Fa, Fin, ---» Fin; --- Mo, daß die Glieder der sechsten Ordnung unentbehrlich sind, wenn 
der Fehler der Funktionswerte nur die Größenordnung 0,01 haben darf, daß aber anderer- 
seits eine Weiterführung der Rechnung bis zur 7. Stufe, wodurch sich die Rechenarbeit etwa 
verdoppeln würde, keinen wesentlichen Gewinn an Genauigkeit verspricht. 

Anhaltswerte für die zweckmäßige Größe der Variationen von f, 9, h ergaben sich aus 


F 
den Ableitungen von F\,,.... Aus den so erhaltenen Durchschnittswerten für > f 135 zn 


wurden durch Umkehrung der Matrix die Werte für ni =fp, ...., hy. Diese liefern zu 
(6) 











XUM 
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den Werten der Funktionen F,,, ... die Verbesserungen von f,g,h. Eine Wiederholung 
des Rechnungsganges führte zu den Werten 
f=-02475; g=-10569; h=+ 1,3548, 
deren Durchreehnung die Ergebnisse 
Fo =+0,0003; Ga = + 0,0000; Ha = + 0,0002; | 
Au=+0%3; Bu=+07668; mu=+1359 | 


(6.5) 


brachte. Aufzeichnung der Reihen F%), ..., mix) gestattete eine, wenn auch rohe Schätzung 
der Restglieder: ihre Größenordnung liegt für F(0), @(0), H(0) bei etwa — 0,01, für m bei 
etwa — 0,025. Der Fehler der aus den Größen (6.4) berechneten Funktionswerte F(£), @(£), H(d) 
kann demnach bis zu etwa + 0,03 ansteigen. 


$ 7. Genauere Lösung durch Mitverwendung der Potenzreihen. Eine größere Genauigkeit 
läßt sich erreichen, indem man für O<{<1 die Potenzreihen (2.1) verwendet, für I<{<x 
die asymptotischen Reihen; denn wegen der Faktoren e-*«> konvergieren diese bei {=1 
schon wesentlich besser als bei {=0. Als Ausgangswerte wurden nicht die zuletzt erhaltenen 
Werte (6.4), (6.5) genommen, sondern benachbarte Werte auf Grund einer Schätzung des Ein- 
flusses der vernachlässigten Restglieder. Für die Exponentialreihen wurden gewählt: 


f= 0,3; g=—10; h=+ 1,34; 
damit wurden an der Stelle ©=1 die Werte 
F(l), (1), H(l), F’(1), @' (1) 
berechnet. Für die Potenzreihen wurde ausgegangen von 
A=+0,%; B=+0,77. 
Damit wurden für ©=1 dieselben Größen berechnet; die erhaltenen Werte wurden zum Unter- 
schied mit einem Stern versehen: 
F (1), @(1)*, H(1)*, F’(1)*, @’(l)*. 

Die Differenzen zwischen beiden Wertegruppen 

Zum: ..:  KBeaiu- u... ... 2... ih 
beliefen sich bei den oben für f, g, h, A, B angenommenen Werten auf: 

L,=+ 0.021852; L, = + 0,001417; L, = — 0,031872; 

L,= + 0,013329; L,= + 0,004559. 

Sobald diese Differenzen verschwinden, müssen die Potenzreihen wie die Exponential- 
reihen dieselbe, allen Randbedingungen (1.8) genügende Lösung der Differentialgleichungen (1.7) 
darstellen. Es war also das Gleichungssystem ,=1,=..=1L,=0 nach f, 9, h, A, B auf- 
zulösen. Das geschah wiederum nach dem Newtonschen Verfahren; Durchschnittswerte für 
die Verbesserungskoeffizienten wurden gewonnen, indem die oben gewählten Werte der 5 Para- 
meter einzeln um — 0,01; — 0,01; + 0,01; — 0,01; + 0,01 abgeändert wurden. Dreimalige An- 
wendung des Newtonschen Verfahrens lieferte die Zahlen 

A=+0,941%6;  B= + 0,773242; | 
f=-0,38891; g=-—1,048332; h=+1,348831 | 
die zu den Differenzen 
L, = — 0,000000 ; L, = + 0,000349; L, = — 0,000009 ; 
L, = — 0,000016;; L, = — 0,000022 
führten. Dabei haben die letzten berücksichtigten Reihenglieder (von den Potenzreihen wurden 


die in der Tafel 1 mitgeteilten Koeffizienten verwendet) zu diesen Differenzen noch Beiträge 
von etwa der Größe 


L, : + 0,00017; L,: + 0,00030 ; L,: + 0,00016; 
L,: + 0,00008;  L,:+0,00066 
gegeben. Die Restglieder der Reihen werden höchstens dieselbe Größenordnung haben. Die 
nächsten Verbesserungen wären für 
A : — 0,000069 ; B : + 0,000297 ; 
f: + 0,000129; 9: — 0,000500 ; h:: + 0,000363 , 


(7.2), 
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so daß die Werte (7.2) schon bis auf etwa eine Einheit der dritten Stelle richtig sind. Da- 
mit haben die aus ihnen berechneten Reihenkoeffizienten, die in Zahlentafel 2 zusammen- 
gestellt sind, die entsprechende Genauigkeit, und die daraus berechneten Funktionswerte (Zahlen- 
tafel 3, Bild 1) haben ebenfalls höchstens einen Fehler der Größe 0,002. Zu größeren Be- 
trägen können sich die Ungenauigkeiten der Anfangswerte nicht summieren, wei) alle sonstigen 
Abrundungsfehler vermieden worden sind: die Sinus- und Kosinuswerte sind einer sechs- 
stelligen Tafel entnommen, und im übrigen wurde die ganze Rechnung achtstellig durchgeführt. 








Zahlentafel 2. Zahlentafel 3. 
Die Koeffizienten der Exponentialreihen. Die Funktionen F(£), G(2,H (2). 
ce. 1 Er G(£) H (2) 
i= — 0,2389 g9= - 1.0483 h=-+ 1,3488 

0.0 0.000 0.000 | 0.000 
a= + 0,4449 8= + 0,8934 0.5 | — 0.349 | + 0.383 | + 0.194 
: 1.0 | — 0.484 | +0.736 | + 0.624 
HC, = — 2,0939 HS, = — 0,5080 15 — 0.450 | + 1.014 | + 1.099 
2.0 | — 0.329 | 41.198 | + 1.498 
FC = +0,4057 60 = —0,558 #03 = 0,9118 2.5 0.176 | +1.272 | + 1.746 
3.0 | —0.086 | 41.272 | + 1.849 
FG = — 0,1622 ©C5=+ 0,1024 #05 = — 0,1553 35 | +0.067 | +1.218 | + 1.830 
Ft —=+0082 63 —=+03556 H5=-+ 0,1311 4.0 | +0.123 | +1.141 | + 1.732 
45 | +0.137 | + 1.064 | + 1.599 
FC = — 0,0366 60} = — 0,0316  HCI= — 0,0412 50 | +o.121 | +1.008 | + 1.468 

FG —+0,0135  6Ci=+0,0242 FC} — + 0,0106 jap u ru ei: 

i n En er % i 5.5 | + 0.088 | + 0.961 | + 1.362 
"5 =+0004 6% = 0,075 HS = — 0,0046 6.0 | +0.050 | 40.943 | 41.294 
FC —=+0,0162 «CI = 400111  #C1 = -+.0,0189 .- ı+ 2. + nn y Zu 
rs, —+0,0203 61 —=--0,058 HS} = + 0,0107 a Tess | H res 
FC; = +0,0013 603 =+-0,013 HC} = — 0,0003 8.0 | — 0.027 | +0.986 | + 1.300 
FSs — — 0,0017 €S3—= — 0,0019 45} = — 0,0012 85 | -—- 0.024 | +0.999 | + 1.326 

x 90 | —0.018 | +1.008 | + 1.347 
FC = + 0,0065  6C5=+ 0,051 HC; = + 0,0049 95 | — 0.010 | +1.012 | + 1.361 
FC; = +0,0002 603 = — 0,0059  H#HC}= — 0,0013 10.0 | — 0.003 | + 1.012 | + 1.368 
FS3 — — 0,000 653 —=— 0,007 #453 —= — 0,021 105 | +0.002 | + 1.010 | + 1.368 
FC = — 0,001  C5=+ 0,0000 BC; = — 0,0001 11.0 | +0.005 | + 1.007 | + 1.365 
FSt — — 0,0008 #54 —= — 0,0002  HSt = — 0,0001 115 | +0006 | + 1.008 | + 1.359 


12.0 | -+0.005 | + 1.000 | + 1.354 
12.5 | +0.004 | +0.998 | + 1.349 


© | 0.000 \ + 1.000 | + 1.349 


Zu den Werten der Zahlentafel 3 ist noch zu bemerken: Die ersten sechs Reihenglieder 

geben nicht genau den Sollwert 0; ihre Summen sind vielmehr 
Fo=+0.057;  Go=+003; Ha =+ 0.0030. 

Daher wurde für das Restglied die negative 6. Teilsumme als Wert angenommen, so daß sich 
durch seine Hinzufügung bei {=0 die Funktionsweıte 0 ergaben. Nun enthalten die Sum- 
manden, aus denen die Restglieder sich zusammensetzen, alle den Exponentialfaktor e-«° in 
7. oder höherer Potenz; daher wurde bei den folgenden Argumentwerten als Ersatzwert für 
das Restglied sein mit e-?@° vervielfachter Wert bei {=0 hinzugefügt, wodurch aber für 
=2 die vierte Stelle nicht mehr beeinflußt worden ist. Alle Werte der Zahlentafel 3 wurden 
dann auf drei Stellen abgerundet. Die drei Funktionen entfernen sich für £> 12.5 nur noch 
weniger als um 0.005 von ihren Endwerten 0, 1, h. 


$ 8. Erörterung der Ergebnisse. Die horizontale Geschwindigkeitskomponente ist nach 
(1.3), (1.6) durch u=eroF(}), v=rwG(£) gegeben, der Winkel zwischen Strömungsrichtung 
und Radius ist also nur von der Höhe 2=yr/»[ abhängig; die Projektionen der Strömung 
sind daher in jeder Höhe logarithmische Spiralen vom Steigungsmaß “-76 gegen die 
Tangente. Die mit der Höhe sich ändernden Strömungsrichtungen sind in Bild 2 dargestellt. 
Am Boden ist die Abweichung von der Tangentenrichtung einwärts am größten: die Steigung 

0 

beträgt dort An = 
weichung auswärts findet sich bei {=4,63; ihr Steigungsmaß ist + 0,129, was einem Winkel 
von 7,4° entspricht. Demnach sind die Strömungsrichtungen in höheren Schichten bis zu 58° 
gegen die Strömungsrichtung am Boden nach auswärts gedreht. — In noch größerer Höhe 


—m=— 1,218, entsprechend einem Winkel von 50,6° Die größte Ab- 
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Bild 2 (rechts). 

Die Horizontalkomponenten der 
Geschwindigkeit über einem 
Punkte P der Grundebene als 
Funktion der Höhe S. 
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dreht die Strömung wieder einwärts und erreicht bei etwa © —=8,06 noch einmal einen Extrem- 
wert mit einer Abweichung von 1,6° einwärts. Fortan wechseln mit weiterem Anfsteigen um je 


A2=Yrlo A= Yrlo 73,60 Yrlo REG EEE A PRFZTF Re 


größte Abweichungen auswärts mit solchen einwärts ab, doch nimmt ihr Betrag jedesmal mit 
a 


dem Faktor e ?— 0,21 ab, so daß alle diese weiteren Abweichungen zu vernachlässigen 
sind. In diesen großen Höhen hat die Strömung praktisch schon überall die durch (1.8.2) 
geforderte Richtung der Tangente an die Kreise um die Achse r—=0. 

Auch die senkrechte Geschwindigkeitskomponente ist vom Ort auf der Grundfläche nach 
(1.3), (1.6) nicht abhängig, sondern nur von der Höhe. Die zu einer Zeit f, in einer waage- 
rechten Ebene in der Höhe 2, =yr/»Z, befindlichen Flüssigkeitsteilchen erreichen daher alle 
zu gleicher Zeit t, die Ebene in der Höhe z, yv/o» £,; die darüber verstrichene Zeit ist durch 
1 | dt 


o) 
> 
gegeben. Die stärkste Aufwärtsbewegung findet in der Höhe Z=3,31 statt: dort ist 
w=1853y vo. 
Während des Aufsteigens von der Höhe z, zur Höhe 2, bewirkt die Strömungskompo- 
nente » eine Drehung um die Achse von der Größe 
da 
LTR 
a, = \ id! ee sa na ee a © 
A 
Hier wie auch in (8.2) kann für £, nicht der Wert 0 eingesetzt werden, weil am Boden nach 
(1.7.3), (1.8.1) H’(0)—=0 ist, die Integrale also divergieren würden. — Setzt man in (8.3) für 
G(£) und H(£) ihre Grenzwerte ein, so erhält man den in großen Höhen bei einem vollen 
Umlauf um die Achse zu überwindenden Höhenunterschied zu 


dz2=yrla Al=yrlo2ah=867 Yılm . .» 2» 222.2. 084; 

die währenddessen vergangene Zeit ist At=2xr]/w. Diese während eines Umlaufs um die Achse 
> 7 

erstiegene Höhe ist zu unterscheiden von der oben angegebenen Höhe 7 Volo = 7,% Vrio, 


um die sich ein Teilchen gehoben hat zwischen zwei gleichsinnigen stärksten Abweichungen 
seiner Bewegungsrichtung von der genauen Tangentenrichtung; unterdessen vergehen also 


0,83 = rund : Umläufe. 185 








252 Bödewadt, Die Drehströmung über festem Grunde z.un ee 





so daß die Werte (7.2) schon bis auf etwa eine Einheit der dritten Stelle richtig sind. Da- 
mit haben die aus ihnen berechneten Reihenkoeffizienten, die in Zahlentafel 2 zusammen- 
gestellt sind, die entsprechende Genauigkeit, und die daraus berechneten Funktionswerte (Zahlen- 
tafel 3. Bild 1) haben ebenfalls höchstens einen Fehler der Größe 0,002. Zu größeren Be- 
trägen können sich die Ungenauigkeiten der Anfangswerte nicht summieren, wei) alle sonstigen 
Abrundungsfehler vermieden worden sind: die Sinus- und Kosinuswerte sind einer sechs- 
stelligen Tafel entnommen, und im übrigen wurde die ganze Rechnung achtstellig durchgeführt. 








Zahlentafel 2. Zahlentafel 3. 
Die Koeffizienten der Exponentialreihen. Die Köaeennd FO,CK,HE). 
t I ro | co | ao 
= — 0,2389 g9= — 1.0483 h=-+ 1,3488 | | | 
0.0 | 0.000 | 0.000 | 0.000 
| 

a= + 0,4449 3=-+ 0,8934 05 | — 0.349 | +0.383 | + 0.194 
, 1.0 | — 0.484 | + 0.736 | + 0.624 
HC, —= — 2,0939 AS, = — 0,5080 15 | — 0.450 | +1.014 | +1.099 
2.0 | — 0.329 | |+119% | + 1.493 
FG = 404057 093= —0,0558 #03 = + 0,9118 25 | — 0.176 | | + 1.272 | + 1.746 
3.0 | — 0.0386 | +1.272 | + 1.849 
FC —=— 0,162 605 =+0,1024 HC; = — 0,1553 35 | +0.067 | +1.218 | + 1.830 
FS—=-+00182 61—=+023556 #51 + 0,1311 4.0 | +0123 | +1.141 | + 1.732 


45 | + 0.137 | + 1.064 | + 1.599 


5, Be 3 0, DESEREBRE s H0% — —_ 2 | | 
[0 0,0366 [0 0,0316 1C9 0,041 5.0 | 40.121 | + 1.002 | 11.468 


FG —+00135  6C3—=-+0,0242  #Cr—= + 0,0106 


Jet Ber: malte ur a ger Jbe- ber 
ri +oo e-Foa nc—+oo Be m re 
Pe —+ 0,0203 6 1 —--0,058 AS — + 0,0107 5 |_o0m | 40.969 11275 
FC; = + 0,0013 6C} = +-0,0013 HC = -——- 0,0003 8.0 — 0,027 | | +.0.986 | + 1.300 
rss — — 0,0017  653—= -— 0,0019 483 = — 0,0012 85 | — 0.024 | +0.999 | + 1.326 
90 —0.018 | +1.008 | + 1.347 

FC = +0,0065 005 —=+ 0,051 HCH= + 0,0049 95 | — 0.010 | +1.012 | + 1.361 
FC} = +0,0002  6C}= —0,0059  4C}— — 0,0013 10.0 | — 0.003 | +1.012 | + 1.368 
rS—=— 0,000 655 =— 0,0027 #5 = — 0,021 105 | +0.002 | + 1.010 | + 1.368 
FC = — 0,0001  6C5=+ 0,0000  3C5= — 0,0001 11.0 | + 0.005 | + 1.007 | + 1.365 
rSt = — 0,008  6St—= — 0,0002 HS = — 0,0001 115 | +0006 |, + 1.003 | -+ 1.359 


12.0 | +0.005 | 41.000 | + 1.354 
12.5 | +0.004 | +0.998 | + 1.349 


| 
+ 1.349 


© | 0.000 | + 1.000 

Zu den Werten der Zahlentafel 3 ist noch zu bemerken: Die ersten sechs Reihenglieder 

geben nicht genau den Sollwert 0; ihre Summen sind vielmehr 
Fo = + 0.057 ; Go =+ 0.0035 ; Ha = + 0.0030. 

Daher wurde für das Restglied die negative 6. Teilsumme als Wert angenommen, so daß sich 
durch seine Hinzufügung bei {=0 die Funktionsweıte 0 ergaben. Nun enthalten die Sum- 
manden, aus denen die Restglieder sich zusammensetzen, alle den Exponentialfaktor e-“° in 
7. oder höherer Potenz; daher wurde bei den folgenden Argumentwerten als Ersatzwert für 
das Restglied sein mit e-?7«° vervielfachter Wert bei {=0 hinzugefügt, wodurch aber für 
E22 die vierte Stelle nicht mehr beeinflußt worden ist. Alle Werte der Zahlentafel 3 wurden 
dann auf drei Stellen abgerundet. Die drei Funktionen entfernen sich für £> 12.5 nur noch 
weniger als um 0.005 von ihren Endwerten 0, 1, h. 


$ 8. Erörterung der Ergebnisse. Die horizontale Geschwindigkeitskomponente ist nach 
(1.3), (1.6) durch u=roF(}), v=rwG(£) gegeben, der Winkel zwischen Strömungsrichtung 
und Radius ist also nur von der Höhe z=yrv/w{ abhängig; die gg der Strömung 
sind daher in jeder Höhe logarithmische Spiralen vom Steigungsmaß —- ee = gegen die 
Tangente. Die mit der Höhe sich ändernden Strömungsrichtungen sind A Bild 2 dargestellt. 
Am Boden ist die Abweichung von der Tangentenrichtung einwärts am größten: die Steigung 
beträgt dort Fo —m=-—1,218, entsprechend einem Winkel von 50,6°. Die größte Ab- 
weichung auswärts findet sich bei {=4,63; ihr Steigungsmaß ist + 0,129, was einem Winkel 
von 7,4° entspricht. Demnach sind die Strömungsrichtungen in höheren Schichten bis zu 58° 
gegen die Strömungsrichtung am Boden nach auswärts gedreht. — In noch größerer Höhe 
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dreht die Strömung wieder einwärts und erreicht bei etwa © =8,06 noch einmal einen Extrem- 
wert mit einer Abweichung von 1,6° einwärts. Fortan wechseln mit weiterem Anfsteigen um je 


IT 


Az=yrlio A{=Yrlo B BB 3, Se BERE 
größte Abweichungen auswärts mit solchen einwärts ab, doch nimmt ihr Betrag jedesmal mit 


dem Faktor e "?— 0,21 ab, so daß alle diese weiteren Abweichungen zu vernachlässigen 
sind. In diesen großen Höhen hat die Strömung praktisch schon überall die durch (1.8.2) 
geforderte Richtung der Tangente an die Kreise um die Achse r—=0. 

Auch die senkrechte Geschwindigkeitskomponente ist vom Ort auf der Grundfläche nach 
(1.3), (1.6) nicht abhängig, sondern nur von der Höhe. Die zu einer Zeit i, in einer waage- 
rechten Ebene in der Höhe 2,=yr/»Z, befindlichen Flüssigkeitsteilchen erreichen daher alle 
zu gleicher Zeit t, die Ebene in der Höhe 2,yv/o» £,; die darüber verstrichene Zeit ist durch 


»2 
1( di \ 
1-4, #0 BE FE ET a a a 
1 
gegeben. Die stärkste Aufwärtsbewegung findet in der Höhe {=3,31 statt: dort ist 
w=1853yvo. 
Während des Aufsteigens von der Höhe z, zur Höhe z, bewirkt die Strömungskompo- 
nente ® eine Drehung um die Achse von der Größe 
6a 
[ G ta > 
er |: ee REN 
a 
Hier wie auch in (8.2) kann für Z, nicht der Wert 0 eingesetzt werden, weil am Boden nach 
(1.7.3), (1.8.1) H’(0)=0 ist, die Integrale also divergieren würden. — Setzt man in (8.3) für 
G(£) und H(£) ihre Grenzwerte ein, so erhält man den in großen Höhen bei einem vollen 
Umlauf um die Achse zu überwindenden Höhenunterschied zu 


4dz=yvl» Al=yrloaah=8607 Yılm -. . 2 2 22.2. 89; 
die währenddessen vergangene Zeit ist At=2n/w. Diese während eines Umlaufs um die Achse 
erstiegene Höhe ist zu unterscheiden von der oben angegebenen Höhe T Vrlo® = 7,20 Yrlo, 


um die sich ein Teilchen gehoben hat zwischen zwei gleichsinnigen stärksten Abweichungen 
seiner Bewegungsrichtung von der genauen Tangentenrichtung; unterdessen vergehen also 


0,833 = rund : Umläufe. 185 
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Gasströmungen mit Übergang von Unterschall- 
zu Überschallgeschwindigkeiten. 


Von H. Görtler in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


1. Überblick. Es gibt brauchbare Methoden zur Berechnung von stationären kompres- 
siblen Potentialströmungen mit kleinen Machschen Zahlen M<1, noch bessere Methoden 
zur Berechnung reiner Überschallströmungen, aber es gibt noch keine brauchbaren Verfahren 
zur Berechnung von Gasströmungen, welche die lokale Schallgeschwindigkeit durchschreiten 
[1], [2], [5], [6]. Das liegt einmal daran, daß, sobald irgendwo die lokale Schallgeschwindigkeit 
überschritten worden ist, Verdichtungsstöße auftreten können, worauf die mit analytischen 
Ansätzen operierenden Verfahren versagen; zum anderen liegt der Grund selbst bei völlig 
regulärem Verhalten der Strömung in der mathematischen Schwierigkeit, eine nichtlineare 
partielle Differentialgleichung (die quasilineare Differentialgleichung zweiter Ordnung für 
das Geschwindigkeitspotential der Strömung) zu behandeln, welche beim Durchgang durch 
die Schallgeschwindigkeit vom elliptischen zum hyperbolischen Typus wechselt, und dies 
auch noch an einem Orte, der zunächst unbekannt ist und bei der Rechnung mitbestimmt 
werden muß. 

Selbst die Anzahl der bisher gefundenen speziellen Beispiele solcher Strömungen ist 
sehr gering [#4], [7], [8], [9], [3] und man wird sich zunächst damit begnügen, diese kleine 
Sammlung durch weitere Beispiele zu erweitern und so die geschilderte Lücke notdürftig zu 
überbrücken. 

In der vorliegenden kleinen Mitteilung führt ein Iterationsverfahren zu einer gut kon- 
vergierenden Lösung für den Fall einer Gasströmung an einer welligen Wand. Als Ausgangs- 
näherung für die Iteration wird jene Lösung benutzt, die sich aus der bekannten, zur 
Prandtlschen Regel führenden Linearisierung der Differentialgleichung für das Geschwin- 
digkeitspotential ergibt. 


2. Das Iterationsverfahren. Unter Berücksichtigung der Bedingung der Drehungsfreiheit 
dü Pv 
rg ie heil an a a u 


(&, v bezeichnen die Geschwindigkeitskomponenten in x- bzw. in der dazu senkrechten y-Richtung) 
lautet die einer kompressiblen Strömung zugrunde liegende Differentialgleichung bekanntlich 


u’\o0u vor 2__dU i 
Es bedeutet c die von der jeweiligen Strömungsgeschwindigkeit abhängige örtliche Schall- 
geschwindigkeit. Es ist zweckmäßig, statt mit @, ® mit dem Geschwindigkeitspotential ® zu 
rechnen, welches vermöge a2, ® -% eingeführt wird. Dann wird (1) selbsttätig erfüllt 
und man hat es nur noch mit (2) zu tun. 

Liegt eine Strömung mit &= U=const, #?=0 vor, so gehört zu dieser eine ebenfalls 
konstante Schallgeschwindigkeit c,. Die zu einer beliebigen Strömungsgeschwindigkeit @, v 
gehörige örtliche Schallgeschwindigkeit ce steht mit c, bei Zugrundelegung des adiabatischen 
Gesetzes in der Beziehung 


14°,  all-— 


Fa %:cP 


ee | N) - . 8) 


„= =. Für den späteren Gebrauch merken wir uns insbesondere an, daß demnach die 
v 

Schallgeschwindigkeit ce, dort, wo die Strömungsgeschwindigkeit die zugehörige örtliche Schall- 

geschwindigkeit gerade erreicht, durch 


GR. DI POER „Sonst! Nike 
lt 5 =! ERENTO v 


gegeben ist. 

Weicht eine Strömung nur geringfügig von der gleichförmigen Strömung & = U= const, 
v=0 ab, wofür wir a=U+tu, #=v mit |u)<Z, |vy|<TU schreiben, so erhält man durch 
Linearisierung von (2): 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. 
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Vai: AT a et 


+ 


T 
Vermöge der Variablensubstitution E=x, n = vyı- 1- - kann man (4.1) für —<1 auf die 





inkompressible Potentialgleichung bezüglich &,n en Dies ist die ; ERROR der 
Prandtlschen Analogie [5l. Kennt man also die zugeordnete inkompressible Potential- 
strömung, so kennt man auch die gesuchte kompressible Potentialströmung — solange die zu 
(4.1) führende Linearisierung zulässig ist. Dies trifft für schlanke Strömungskörper, welche 
in Richtung ihrer Längsausdehnung angeströmt werden, zu. Nur in der Umgebung der Stau- 
punkte ist wegen u=— U die Voraussetzung verletzt. Die genannte Zuordnung wird unmög- 
lich, wenn sich die Machsche Zahl u der ungestörten Strömung dem Werte 1 nähert. 
Berücksichtigt man auch die in «, v und deren Ableitungen quadratischen Glieder, so 
gewinnt man aus (2) unter Heranziehung von (3) die Differentialgleichung 
U’\ou Uudu Uu dv U .. 
1-, tn rt 3y+r27 Fady 
während die a auch der kubischen Glieder auf die vollständige und völlig 
strenge Differentialgleichung (2) führt, welche sich auch schreiben läßt: 


1-- 2 tr le+n Dur@tD.,® a2 


(4.2), 


dx 0dy c.9xr 2 





1 er Er | 

+ arriG yuu+ tz w s e) ER EN 
2duf,. \ 

+osgajlo+unl 


Wir haben diese beiden Gl. (4.2) und (4.3) so geschrieben, daß rechts beidemal jene Glieder 
stehen, die gegenüber (4.1) hinzukommen. 

Der Gedanke der Herstellung einer Lösung durch Iteration wird nun so verwirklicht, 
daß zunächst das ins Auge gefaßte Problem auf Grund der linearisierten Gl. (4.1) gelöst wird. 
Die erhaltene Lösung ®, bzw. w,, v, wird als erste Näherung benutzt und in die rechte 
Seite von (4.2) eingeführt. Die Lösung der so entstehenden linearen inhomogenen Differential- 
gleichung liefert die zweite Näherung ®, + ®,. Für den nächsten Iterationsschritt hat man 
bereits die volle Gl. (4.3) zu verwenden, denn die in der rechten Seite von (4.2) entstehenden 
neuen Störungsglieder für eine dritte Näherung wären bereits von derselben Ordnung wie in 
(4.3) hinzukommende kubische Glieder. Alle weiteren Iterationsschritte spielen sich also in 
der vollen Gl. (4.3) ab, wobei jeweils die zuletzt berechnete n-te Näherung in die rechte Seite 
einzusetzen ist, worauf die neue (» + 1)-te Näherung aus der entstehenden linearen inhomo- 
genen Gleichung zu berechnen ist. 

Ist die Größenordnung der auf U bezogenen ersten Korrektur «,, v, der Geschwindigkeits- 
komponenten durch den kleinen Parameter a gegeben, wird also die erste Näherung durch 
einen in a linearen Ausdruck gegeben, so bringt die zweite Näherung zusätzliche in a qua- 
dratische Glieder, die dritte Näherung bringt Terme in a® aber auch zugleich höhere Potenzen 
in.a. Es ist zweckmäßig, für die dritte Näherung nur Glieder bis zur Potenz a’ auf der 
rechten Seite von (4.3) zu berücksichtigen und die höheren Potenzen zu streichen. Die ge- 
strichenen Glieder in a* kommen dann mit allen weiteren Gliedern 4. Grades in «a zusammen 
bei der Berechnung der vierten Näherung herein. Verfährt man entsprechend bei allen 
weiteren Schritten, so erhält man bei jedem Iterationsschritt gerade ein gesamtes weiteres 
Reihenglied in der Potenzreihe der Lösung nach dem Störungsparameter a. In dieser Weise 
soll im folgenden verfahren werden. 


3. Beispiel: Strömung längs einer welligen Wand. Es kam uns bei den im folgenden dar- 
zulegenden Rechnungen lediglich darauf an, irgendein einfaches Beispiel praktisch durch- 
zurechnen und damit das Iterationsverfahren zur Ergänzung der einleitend genannten Beispiel- 
sammlung heranzuziehen. Die Rechnung wurde daher darauf angelegt, einen gut konver- 
gierenden Fall zu finden. Man konnte erwarten, daß eine Strömung, welche um die gleich- 
förmige Strömung %=U=const, v—=0 oszilliert, ohne durch vorhandene Staupunkte an 
einzelnen Stellen allzu stark von dieser abzuweichen, für unser Iterationsverfahren besonders 
geeignet ist, da dann bereits die linearisierte Gl. (4.1) in allen Punkten der Strömung eine 
Näherung liefert. Es wurde daher das Beispiel einer Strömung an einer welligen Wand ge- 
wählt. 
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Eine Lösung von (4.1), welche eine solche Strömung darstellt, ist bekannt und lautet: 


BERN. 9; 
dB=-Uxtae U reos = EN Ne 


U ' 
worin «=y1—M?, M dieMachsche Zahl der ungestörten Strömung: M= <, ! die Wellen- 


o 

länge der Wandwellen und «a, eine Konstante bezeichnen. Für M<1 — und wir beschränken 
uns zunächst hierauf — erhält man eine Stromlinienverteilung mit schwachen Wellen, 
welche mit wachsendem % abklingen und eine gute Näherung darstellen, solange «a, hin- 
reichend klein ist. Man kann etwa die um y=0 oszillierende Stromlinie als Wand wählen 
und sich auf die Strömung oberhalb dieser Wellenlinie beschränken. 

Bevor wir die Näherung (5) als erste Näherung für das oben geschilderte Iterations- 
verfahren verwenden, führen wir folgende dimensionslose Größen ein: 


2 2 
[= I %, n= ı 9° 
(6). 
_2a® _?a 
Be 05° Ba © u 


Darin ist ec, aus (3a) zu entnehmen. Die sich ergebenden Geschwindigkeiten in der weiteren 
Rechnung sind gemäß der Definition von @ auf c, bezogen. 
Die Iteration liefert schrittweise die aufeinanderfolgenden Glieder der Potenzreihe 


E; z - 
-rnrentrent rn re ehe ee 


wobei zur eindeutigen Festlegung der p,(&,n) die Erfüllung von Rand- oder sonstigen die 
gewünschte Lösung charakterisierenden Bedingungen erforderlich ist. Da es uns, wie bereits 
oben bemerkt, nur auf die Herstellung eines gut konvergierenden Beispiels ankam, haben 
wir nicht die Wand in bestimmter Weise vorgegeben, sondern die gesuchte Lösung eindeutig 
durch die Forderungen charakterisiert): 
P&EO=I:+acost, PO) EEE 
1 1 


Diese Bedingungen werden bereits durch die erste Näherung — €+.ap, erfüllt. Man hat also 
1 
3 


Das Iterationsverfahren liefert bis zur dritten Näherung einschließlich das Ergebnis: 


zu fordern, daß 9, (&,0)=0, p,(&,®)=0 für v=2,3,... 
p, = A,e-""cosE, 

Y=A,ne-?""sin2E,  ' 
Ps = rA, (ent — ent) +B,nem?#,cos&+XC,n+ D,n?ye-3#"cos3£ 


Darin sind die Koeffizienten gegeben durch: 


4,=®, 

PEN u Ei 

au Se ET 
oz - 

= at Ba +H PM + HB EM AR +1) u 

8 —3)u+12(& +1) u, > 

co, (2 +1) M* . 

Be Rt -e du, 
e 1)? M® 

ru u? {M’+4u}, 

_&&@+l’M 

a 





2) Die erste Bedingung schreibt eine um U oszillierende Verteilung von %(&,n) längs 7=0 vor: U—%(&,0) - sin. 
Die analoge Vorschrift © (5,0) - cos$ führt nach unseren Rechnungen zu einer schlechter konvergierenden Eutw icklang (7), 
mit der wir nach 3 Schritten (bis 93 einschließlich) keine Ergebnisse erzielen konnten, bei denen die örtliche Sehall- 
geschwindigkeit in mehr als unmittelbar wandnahen Gebieten überschritten wurde. 
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Wie man auf Grund der Ausgangsnäherung erwarten mußte, versagt die Lösung für 
u—=0(M=1). Je größer der Störungsparameter « gewählt wird, desto weiter muß « von 
Null, d. h. die Machsche Zahl M der ungestörten Strömung von 1 entfernt bleiben, damit 
die abgebrochene Reihe eine gute Näherung darstellt. 


Als numerisches Beispiel wurde folgender Fall durchgerechnet: 
a=02, M=09(u=0,4359), »=1,405. 


Als Wand wurde eine um 7—=0 oszillierende Stromlinie angenommen. Es ist e,= 0,9839 ce, 
T \ . 


I 
(somit M=- 09148). 


ce 


r ü 0« v 06 - a i EL n 
Wegen ——= =, ——= ” berechnen sich die Geschwindigkeitskomponenten aus der 


. Be 0 


oben gewonnenen Näherung für @ in der Gestalt: 


(10), 


U Be 4 
= +l(au, nD)+a’u,(n))sinö+ta’u, (n)cos2£-+-a’u, (n)sin3E, | 
1 
— — (av, (n)+ av, (m)) cos&+a®v,(n)sin2£+ av, (n)cos3E | 


also in Form einer Fourier-Reihe | 


bezüglich der Variablen &, wobei = | \ 
höhere Näherungen des Iterations- Rene: au men a 
verfahrens, wie man sieht, noch ent- | 

sprechend kleine Korrekturen der 2 Pr. | 
Koeffizienten bringen können. Um | 
einen Überblick über die gute Kon- In, | | | | 
vergenz dieser Entwicklungen bei -008 LH 


den gewählten Parameterwerten zu N” | | m | 

geben, sind in Bild 1a und 1b die #1 | 0 | 

Kurven für die von n abhängigen - v7 u 

Koeffizienten der obigen abgebroche- / | | | | 
Es 
| 

















nen Fourier-Reihen dargestellt. 


Bild 2 gibt den Verlauf der 
Stromlinien, Bild 3 den Verlauf der Pf 
Kurven gleicher auf die Schall- SWGERER I I 
geschwindigkeit ce, bezogener Ge- Bild 1a. 
schwindigkeitsbeträge 


-016 














Qt ——— + 11 





u 
N ü? +9” = const 
1 


wieder. Innerhalb des von der ge- 7 en ER 
strichelten Linie in Bild 2 ange- | —— 
gebenen Bereiches herrscht Über- 12) | Br 
schallgeschwindigkeit. Einige Mach- - 002 | a 


2. . 1 j ’ av, 
sche Linien sind eingezeichnet. In . nn PR 
Bild 4 ist schließlich die Wandstrom- | (2) &®%, 

4 


linie stark überhöht dargestellt. m 4 





0 
{, 











Die auftretenden symmetri- 
schen, stark in die Strömung hinein- "08 Tu 
greifenden Überschallgebiete erinnern A 
an jene, welche für Düsenströmungen 
an den Wänden in der Gegend des 
engsten Querschnittes berechnet 
worden sind (G. I. Taylors „sym- un 7 < J AR ERE 6 7 
metrische* Lösungen [8], siehe auch Bild 1b. 
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6 [3). Dort konnten jedoch die Über- 
schallgebiete nur in grober Näherung 
angegeben werden. 

5 Bei festem « und mit gegen 1 
wachsenden Machschen Zahlen M 
läßt die Güte der Konvergenz obiger 
Entwicklungen schnell nach. Das 

















e steht insofern in Einklang mit den 

Beobachtungen, als von einer ge- 
7 wissen Machschen Zahl ab Ver- 
34 : dichtungsstöße auftreten (siehe hier- 








zu die in [6] wiedergegebenen Auf- 
nahmen). Inwieweit sich solche Über- 
schallgebiete ausbreiten können, bis 
Verdichtungsstöße auftreten, wird 
sehr wesentlich von der Formgebung 
der Wandwellen abhängen. Dafür 
spricht auch die sehr starke Ab- 
hängigkeit der Konvergenz unserer 
Entwicklungen von der Gestalt der 
Randbedingungen (siehe hierzu Fuß- 
note ?)). 


4. Ein weiteres Beispiel: Strö- 
mung zwischen entgegengesetzt wel- 









































4,7 r . . 
ligen Wänden. Wenn man die in Ab- 
schnitt 3 wiedergegebenen Rechnun- 
gen entsprechend für Machsche 

® Zahlen 

U 
7 M=—>1 
177 

a durchführen will, hat man an Stelle 
von (5) von der Ausgangsnäherung 
(Lösung der linearen Gl. (4.1)) 

! x BR. 

P=Ux+ta, cos ] (uy-+ x) 
(11) 
we 

A ! Pr (d=iu=yM? —1) 

. A . 7 . ig 
Ai $) u - ch nd ad auszugehen. Die Stromlinien stellen 
Bild 3. Kurven konstanter Geschwindigkeitsbeträge. also nach dieser ersten Näherung 
Wellenlinien dar, deren Amplitude 


057 er 7 777 nicht mit dem Wandabstand abklingt; 

| es handelt sich vielmehr um Wellen- 
linien, welche zur Wandwelle kon- 
gruent bleiben, aber gegen diese 
parallel versetzt sind längs Linien, 
die mit der ungestörten Strömungs- 
richtung den Machschen Winkel 
a=are cotg «u einschließen (s. [5] 
S.209). Die entsprechend Abschnitt 3 
durchgeführte Rechnung liefert eine 
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we; “ Entwicklung, welche wegen der auf- 
Bild 4. Wandwelle in lWfacher Überhöhung. tretenden n-Potenzen nun nur in 


nächster Wandnähe konvergiert. 

Um daher einen in n-Richtung begrenzten Fall zu behandeln, werde eine Gasströmung 
zwischen entgegengesetzt welligen Wänden (unendliche Aneinanderreihung von Düsen) be- 
trachtet. Auf Grund der Symmetrie der Anordnung kann man statt dessen auch die Strömung 
zwischen einer welligen und einer ebenen Wand ins Auge fassen, wie sie der Versuchs- 
anordnung in [6] entspricht. 

Die Ausgangsnäherung nach der Prandtlschen Analogie für einen solchen Fall erhalten 
wir durch Überlagerung von (11) und der bezüglich y—0 gespiegelten Lösung. Unter Ver- 
wendung der Dimensionslosen (6) lautet dieselbe: 
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Y- 


fe 


.% 
Sta pP, 
i 
mit (12) 
Co, = 
Y%, > r cos un COSS 
1 
Um eine gut konvergierende Entwicklung bei der nun vorzunehmenden Iteration zu erhalten, 
verlangen wir, daß die gesuchte Lösung eine Geschwindigkeitsverteilung u (£,0) x sin & längs 
der Symmetrielinie 7 =0 liefert. 
Durch Iteration ergibt sich dann für die Entwicklung 


Es 
U or 
. £ v £ 
Tr a @y(5,n) 
1 — 
v=] 


mit einem Schritt: 


9—=kCnsin2un+D(l—cos2iun)\sin 28, 


wobei C=M° nu 2 ; 
c, 16 u Bl a Ser een 272 
u: 4 re 
D=N- {4 — (3— x) M?} 


162‘ 
Für den nächsten Iterationsschritt hat man die Differentialgleichung zu lösen : 


2 12 
En Wa 
a! Ö 7" 


cos& fa, cosun+b,nsinun+te,cos3un+d,nsin3un) \ 




















+ 0083 & {a,cos un +b,nsinun+c,cos3un+ d,nsindnn | 19) 
mit 
1 u . RER EN un | ‚| 
a=—5|M{D(4z—- #2 —-1)M’)— Cu(2x ),+, W 76% PD M, 
-_ 1 > 
B .i 
=; M(a+1)C, 
AR m ojl_ +3), 3@+0 
> 5a u(3—2x)—- D(4 —-(3— x) M’)} + c, \3 4 M?-+ 6 u|, 
1 
d,=>MC{4-(@-x)M’,, 
1 “ = (+1 M?\ 
w—_ 5 |nener 5) Du (7 2) +3@+ 1), +0 °F efı = ). 
2 : 
b,=—3d, 
1 PR „+1 
=; |MICH@=+N) D2—Bx+5) my + .z m “ 
_ 1 > 
d,= 35,. 
Die Lösung lautet 
" (& ie a 2 2320 
PP, —C0S ES 2u 24 .) n sın un] — 41 n’cosun 
c N Er d, u 
+5 a r 39 2) (cos un — cos 3 un) — RyTz n sin 3 u ı ; 
(15). 
Bart d, b, & u b, dä 
+cos3E SITE >) (cos un —eos3un)+ ge nsinun 
en a 0 e 
+ F +3 =) nsindun 12% N’ cos3un 








Von den erhaltenen Ausdrücken gelangt man zu der entsprechenden Lösung für M<1 in 
reeller Schreibweise, wenn man überall « durch i « ersetzt. 
Die Geschwindigkeitskomponenten @, » sind wieder der Gestalt nach durch die Aus- 
drücke (10) gegeben. 
Bevor wir die erhaltene Lösung an einem numerischen Beispiel erörtern, schicken wir 
eine Bemerkung über die zu erwartende Größenordnung der Verbesserungen a’y, und a’y, 
18* 
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im Vergleich zu «ay, voraus. Solange M und «a so gewählt werden, daß die örtliche Schall- 
geschwindigkeit durchschritten wird, muß nach der einfachen Stromfadenvorstellung der 
Strömung an den beiden Durchschreitungsstellen eine Zusammendrängung der Stromlinien er- 
folgen. Das bedeutet aber, daß der Ausdruck a°’y,, der den Faktor sin 2£ enthält, not- 
wendig gegenüber dem Ausdruck ag, mit dem Faktor cos noch stark ins Gewicht fallen 
muß. Dagegen werden die weiteren Glieder der Entwicklung, beginnend mit dem von uns 
berechneten Term a’, klein ausfallen. 
Wir wählen als Beispiel: 


a=092, M=11(u=048), x=14b .....n (16), 


I 
Dann ist e, =1,018e, und M, = = 1,081. 
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0 — — 
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er | 9) 
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-902 

av, 
N) av, rad, 

(2) eve 

£) v 

— 2 
AB. 5a] | 
| 
Bild 5a (oben) | 
Bild 5b (rechts) 006 1 A Be ' 


Entsprechend dem Bild 1 sind in Bild 5 
dieeinzelnen Beiträgezu den Geschwindig- 
keitskomponenten dargestellt. Sie be- 
stätigen unsere oben ausgesprochene 
Erwartung. -Bild 6 gibt das Feld der 
Stromlinien wieder, Bild 7 die Kurven 
gleicher Geschwindigkeitsbeträge. Die 
Geschwindigkeit schwankt von 0,84, 
bis 1,28 ce, um die „ungestörte* Ge- 
schwindigkeit U= 1,08 e,. In Bild 6 sind 
wieder jene Kurven hinzugefügt, welche 
Überschall- und Unterschallgebiete tren- 
nen°®). Schließlich gibt Bild 8 die von 
uns als Wand gewählte Stromlinie stark 
überhöht wieder. Das Verhältnis der 
Querschnitte an der weitesten und eng- 
sten Stelle beträgt 1.07. 

Wir haben damit eine gut konver- 
gierende Lösung gewonnen, in der sich 
in einer völlig regulären Strömung peri- 
odisch Über- und Unterschallgebiete ab- 
lösen. Der einzelne Übergang entspricht 
jenem der Meyerschen Düsenströmung 
[4]. Während jedoch die bisher für die 
Düse angestellten Rechnungen [4], [S], [3] 
Bild 7. Kurven konstanter Geschwindigkeitsbeträge. auf Grund der schlechten Konvergenz der 





Bild 6. Stromlinienverlauf und Grenzen zwischen Unterschall- 
und Überschallgebieten. 




















3) Bemerkung bei der Korrektur: Die Beschriftung von Bild 6 enthält einen Fehler. Es muß heißen: #>1 
bzw. w<1 statt w>c, bzw. w <e;. 
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benutzten Potenzreihenentwicklungen +2 a 
nach £und y nur Aussagen in der Nähe des 
engsten Düsenquerschnitts gestatteten, | 
gilt die hier gewonnene Lösung längs ’# 
der ganzen £-Achse in dem gewählten 
n-Bereich. 

Es sei an dieser Stelle auf ein 
Tlieorem hingewiesen [3] S. 332/333. 
welches in engem Zusammenhang mit 
unserer Problemstellung stehtund welches ER RR , 
wir bereits früher [3] S. 335 auf den | | 
Fall der Strömung zwischen entgegen "asus -2 -7 o 7 2 3a 
gesetzt welligen Wänden angewandt 
haben. In einer engeren Fassung sagt 
dieses Theorem im zweidimensionalen Falle aus: Erreicht in einer zu einer Geraden spiegel- 
bildlichen in einem gewissen Bereich K regulären stationären kompressiblen Potentialströmung 
die Geschwindigkeit auf der Symmetriegeraden in einem Punkte P des Bereichs die örtliche 
Schallgeschwindigkeit als Extremum, so herrscht auf der ganzen Senkrechten zur Symmetrie- 
geraden durch P in K die örtliche Schallgeschwindigkeit. Die Stromlinien durchkreuzen diese 
Senkrechte orthogonal und unter verschwindender Krümmung. 

Denken wir uns bei einer Strömung zwischen vorgegebenen entgegengesetzt welligen 
Wänden oder zwischen einer welligen und einer ebenen Wand — die Machsche Zahl 








Bild 8 Wandwelle in 20 facher Überhöhung. 


M=— von einem kleinen Werte aus anwachsend, so werden sich, ganz wie in Bild 2 


zunächst symmetrische Überschallgebiete an den Wänden in der Umgebung der engsten 
@Querschnitte ausbilden. Es kann aber nicht der Fall schließlich eintreten, daß bei wachsendem M 
durch das Größerwerden dieser Überschallgebiete bis zur gegenseitigen Verschmelzung an Stellen 
der Mittelachse selbst die örtliche Schallgeschwindigkeit analytisch erreicht wird. Denn 
nach obigem Theorem müßte in jenem Stadium, in welchem die örtliche Schallgeschwindigkeit 
auf der Mittelachse gerade (und das heißt: als Maximum) erreicht wird, überall in den 
Ebenen der engsten Querschnitte Schallgeschwindigkeit herrschen und die Stromlinien müßten 
beim Durchschreiten derselben verschwindende Krümmung haben. Dies ist aber nicht mit 
der Gestalt der Wandstromlinien verträglich, solange man ihnen im engsten Querschnitt 
endliche Krümmung gibt. (Wenn sie dagegen hier verschwindende Krümmung haben, bilden 
sich die erstgenannten symmetrischen Überschallgebiete an den Wänden nicht aus — ebenfalls 
zufolge unseres Theorems in einer allgemeineren Fassung, s. [3])) Es muß also, bevor die 
örtliche Schallgeschwindigkeit auf der Mittelachse erreicht wird, der analytische Zusammen- 
hang der Strömung verlorengehen, wobei wir an das Auftreten von Verdichtungsstößen 
denken. Diese Aussage gilt jedoch nur, insoweit die Strömung nicht von der Wand abreißt. 
Tritt jedoch Grenzschichtablösung ein, so versagen natürlich unsere potentialtheoretischen 
Betrachtungen. 

Es sei zum Schluß noch bemerkt, daß kürzlich A. Busemann anläßlich der GAMM- 
Tagung am 4. Mai 1940 in Braunschweig u.a. über Rechnungen berichtet hat, die in ähnlicher 
Weise wie in der vorliegenden zu jenem Zeitpunkt bereits abgeschlossenen Arbeit von der 
nach der Prandtlschen Regel zu gewinnenden Lösung als Ausgangsnäherung für ein 
Iterationsverfahren ausgehen, hier jedoch aufbauend auf der Differentialgleichung für die 
Stromfunktion und mit der anderen Zielsetzung der Berechnung von Unterschallströmungen 
um vorgegebene Berandungen. 


Zusammenfassung. Mit dem Ziele, Gasströmungen zu untersuchen, welche die örtliche 
Schallgeschwindigkeit durchschreiten, wird ein Iterationsverfahren zur genäherten Berechnung 
kompressibler Strömungen entwickelt. Die Ausgangsnäherung wird auf Grund jener Lineari- 
sierung, welche zur Prandtlschen Regel führt, gewonnen. Gegenüber der Ausgangsnäherung 
nach dem Verfahren von Rayleigh und Janzen, nämlich der inkompressiblen Potential- 
strömung, hat die hier verwendete Ausgangsnäherung den Vorzug, bereits für beliebige 
Machsche Zahlen #1 die Kompressibilität zu berücksichtigen, setzt jedoch voraus, daß die 
Strömung nur gering von einer reinen Translationsbewegung abweicht. Durch Iteration 
wird diese Näherung nun so verbessert, daß sie auch für größere Abweichungen gilt. Das 
Verfahren versagt, wenn die Machsche Zahl der „ungestörten“* translatorischen Grund- 
strömung gleich Eins wird, gibt aber sonst gut konvergierende Entwicklungen im wesent- 
lichen Potenzentwieklungen nach einer Maßzahl für die Abweichung der Strömung von der 
Parallelbewegung. 
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Als Beispiel wird zunächst die Strömung längs einer welligen Wand berechnet. Ein 
Zahlenbeispiel liefert geschlossene, symmetrische, stark in die Strömung hineinreichende 
Überschallgebiete an den Wanderhebungen. Als weiterer Fall wird eine Strömung zwischen 
entgegengesetzt welligen Wänden ermittelt und durch ein numerisches Beispiel illustriert, 
in welchem die Strömung mit beträchtlichen periodischen Schwankungen abwechselnd Unter- 
schall- und Überschallgeschwindigkeiten erreicht. 
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Halbebene mit Kreisbogenkerbe. 
Von €. Weber in Dresden. 


l. Einleitung. Das zweidimensionale Problem der Halbebene, die in Richtung parallel 
zur Begrenzungslinie gezogen wird und eine oder mehrere Kerben hat, gehört zu den einfach 
ausschauenden, aber bisher nicht gelösten Problemen der Elastizitätslehre. Dabei wäre es 
praktisch beinahe gleichgültig, ob die Kerben durch Kreisbogen oder durch Ellipsenbogen 
begrenzt sind, da die maximalen Spannungen im Kerbgrunde überwiegend vom Verhältnis 
Grundhalbmesser zur Kerbtiefe abhängt. 

Ich führe im weiteren die Lösung für Kreisbogenkerben durch und nehme hierzu nach- 
einander mehrere konforme Abbildungen vor, wobei die erste die Inversion ist. Hierbei er- 
hält man ein anderes elastisches Problem, dessen Behandlung einen gewissen Vorzug auf- 
weisen wird. 

Für die Halbebene mit periodischen Kerben habe ich eine Näherungslösung gefunden, 
die darin besteht, daß die Halbebene mit gewelltem Rande auf eine Halbebene mit geradem 
Rande abgebildet wird. Diese Lösung soll später veröffentlicht werden. 


ll. Die konformen Abbildungen. Die Halbebene, 
Bild 1, werde in «-Richtung gezogen. Am Rande be- 
findet sich eine Kerbe, die durch einen Kreisbogen be- 
grenzt ist. Das Problem verwandle ich zuerst durch die 
konforme Abbildunz der Inversion in ein Problem des 
Ebenen-Ausschnittes, d. h. eines Stückes der Ebene, das 
durch zwei sich schneidende Geraden begrenzt ist. Ist 
der Winkel zwischen beiden Geraden ein rechter, so er- 
hält man den Sonderfall der Viertelebene. 

Für einen belasteten Ebenen-Ausschnitt gibt es auch 
noch keine Methode zum Auffinden von Lösungen. 

Durch eine zweite konforme Abbildung wird der 
Ebenen- Ausschnitt in einen Streifen verwandelt (ein Bild 1. 
durch zwei parallele Geraden begrenztes Gebiet). An 
Stelle der Airyschen Spannungsfunktion wähle ich hierbei eine andere geeignete Funktion. 





[a 18# 1) 
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Die Airysche Spannungsfunktion der Halbebene mit Kreisbogenkerbe 
FER[g (re +i)+lr — iy)g(e + iy)] 


gehe im Unendlichen gegen Fz =y’, so daß ox==2 wird. Am Rande müssen F=0 und 


die Ableitung senkrecht zum Rande 2 —0 werden. Der Koordinatenanfangspunkt liegt in 
A, Bild 1. Durch «@+iy'=1/(z-+iy) bilden wir das Gebiet konform auf den Ebenen- 
Ausschnitt mit dem Eckwinkel a> ab. Die Funktion F’(#’,y')=Fi(x,y)/(e +?) ist 
wieder eine Airysche Spannungsfunktion. Am Rande müssen 


a, 


“ OF E ö a 
F'=0 und 7 —=(0 werden. Die unendlich entfernten Punkte der 


Halbebene gehen in den Punkt U’ mit «’=0, y"=0 über (Bild 2). 
In diesem Punkte erhält man eine singuläre Stelle. Die Airysche 
EIER 2 
et N 
Punkte des «-y-Gebietes und des x’-y'-Gebietes ergibt sich aus den 

Bildern 1 und 2. 

Nun verlegen wir den Koordinatenanfangspunkt nach dem Eck- 
punkt B’ dureh Einführung des Koordinaten vr’ d”" =" +iy’—(@’+N)+tiy. Die 
Spannungsfunktion geht über in: 


Spannungsfunktion wird dort Die Zuordnung der 








F’ (x, y’) = F’' (x —1,y”) = F” (&’, y”). 


An Stelle dieser Funktion F’ untersuchen wir im weiteren die Funktion F*=F”/r". Die 
Singularität von F* in der singulären Stelle bleibt dieselbe wie von F”, da dort r” =1 ist; 
auch die Randbedingungen bleiben erhalten. 


Es ist von besonderer Wichtigkeit, daß F* wie folgt geschrieben werden kann: 


> 


*_ B gr a ER. gr | T\ „ 
F*=sınd”o, +sin|9”’+az Ip, » 


wo 9, und 9,” zweidimensionale Potentialfunktionen sind und die Koeffizienten von 9,” 
und 9,” kein r’” enthalten. 


Nun bilden wir durch eine zweite konforme Abbildung den Ebenen- Ausschnitt mit dem 


Eckwinkel a auf einen Streifen von der Breite 5 durch Einführung der Veränderlichen A 


und 9 ab: 
ee (#4 i9) — eia(t ii) M ei iln r’’) ; 1 iy”. 
Dann wird 
’ 967 
. sinald —+ 
Bat nad. - wi | > Ko 
F* (x2”,y") _ Rfld -il)+ Rf(id ir). 
a 


Hierbei ist gesetzt: 
 . a > 082 2 20 
Po en Rd UA), Yı a Rf,(d th). 
Für a=0 werden die angeführten Abbildungen hinfällig; statt des Grenzüberganges 
sei eine unmittelbar zum Ziele führende Abbildung angegeben. Wir setzen in Bild 4 


stHiy=o ‚ und erhalten damit die konforme Abbildung auf den Streifen von der 


1 
+iy 


Breite 5. Die Airysche Spannungsfunktion wird: 


7 ’ ’ ’ o ’< . ’ _ f ’ 7 ıpr ’ ze 
Fi, y)=Fle, I =yRf(y — ie’) + mM +5 IEIVAN, i). 


Diese Gleichung stimmt mit dem Grenzfalle «>0 überein. Auch die Randbedingungen und 
die singuläre Stelle sind dieselben, 
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'ir kehren wieder zum allgemeinen Falle mit a=+ ück. en Rändern des 
Wir kel d allgemeinen Falle mit OÖ zurück. An den Ränd d 
Streifens ist F* gleich null, ebenfalls a. Im Koordinatenanfangspunkte haben wir die 


Singularität 


Y g . gr x t > gr x i 
os sn sn HR —, ui Rz ; 
((e”’ — 1? +9" ?]r” (.«”— h)+iy elnr +19” _ 
und da nur die Umgebung von 9’ —=0, r"—1 zu beachten ist: 
1 sinad .- 1 
- &) 


: _ grg Pr 
sin 9’R : == t —. 
(#” —ilnr”) a F—il 
Im weiteren benötigen wir zwecks kürzerer Schreib- 
weise die verschobenen Koordinatensysteme nach Bild 3. 
Wir nennen die Koordinaten im unverschobenen Koor- Fi 


. fi . 2 P>r . T 
dinatensysteme statt 9,4 nunmehr #,,4, für dieum k 5 Ne S 





verschobenen Koordinatensysteme %;, A; hierbei ist _— 





9 =9d,—kZ mit k=...—2,—1,0,1,2,... Weiter 


führen wir die komplexen Veränderlichen I; = #, —i4 





; ra ; 1 1 \ \ 
ein. Die Funktion R-—=NR 3.277 wird häufiger vor- 
sk ' Su ei 


kommen; wir schreiben hierfür 97. 
Die Singularität im Koordinatenanfangspunkte ist a3) A 





in dieser Schreibweise gegeben durch sin a Do: Po: Bild 3. 


Ill. Aufbau der Lösung. Wir setzen: 


sin ad sina® 


0ıy 119 
F* = + —P, 


a 
wo Y, und Y, zwei noch unbekannte Potentialfunktionen sind. Wir bilden beide aus 
Näherungsgliedern 

JR pn 

/ — Wo T Yon + VINEEN) rt.» 

ara ‚ ’ i 

P,= Y otYyıntWYınt--- 
Als nullte Näherung für Y, nehmen wir y,.,=9,, für Y, nehmen wir y,.»=0. Dann ist 
die nullte Näherung: 


A ER .c1 
Fo rrz sın a d, "Po 


Am Rande 9,=0 sind alle Randbedingungen erfüllt; am Rande #,=0 werden weder F* 
> F.* 
Of. . Fe u . Ba . i R 
noch 39 gleich null. Wir fügen die ersten Ergänzungsglieder y,,, und y,,, hinzu, die 
ı# 
wir so wählen, daß am Oberrande die beiden Randbedingungen F*=0 und 35 —=( erfüllt 
di 


? 2 i i ’ - a i 0g £ 
werden. Dabei beachten wir, daß 9, eine antimetrische Funktion in 9, und cz eine 
ot 
symmetrische Funktion in 9, ist. 


Für y,., und y,., erhalten wir folgende Bedingungen: 


[sin ad, (9 + You) +sinad,-w, «lg, ER ee 
FR 
" ( | r .. ( n . 
5 „sin dt r+snad, Ynlgz, 0: - (ib). 
Hieraus 
ER 
sinagz(%+Wwnga.,=0: - - - - (2a) 
T . 7 'O Pos ı yon 
BR a a a ja r = 9 
a Cc0s a 2 (Pot Won)g, „rsın a 2 | 2% bu. 2% )+ ay 109,9 oe Br ER: (2b). 
Aus (2a) folgt Yu, —g,, di 1 9, für 9,—=0, also 9, = und , = — 5 si 
Aus (2: gt Yun =%;,, da p, und 9, für 9, —=0, also d,=— und 9, = — sich aufheben. 


ed 
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Weiter folgt aus (2b) 
») 


= R IT Ö Ya 
Yo > sın a 


2 9% 


Jetzt fügen wir die nächsten Glieder y,., und %y,., hinzu, die so gewählt werden, daß am 
anderen Rande, also für 9, —=0 die Randbedingungen erfüllt werden und erhalten 
‚9 7 2 d° q ) 


j A» ’ 2 ung „T Ö v 
FT .Ttrit gagG“ ar 5 q 
10m u a 2 hi) u? 


Mom Maz ag: 

So geht es weiter; bei der nächsten Näherung kommt 9, und dessen Ableitungen hinzu. 
Das Ergebnis ist in folgender Tafel 1 zusammengestellt. Hierbei gibt die Summe der unge- 
> 


s BEN or R FE n 
raden Zeilen Y,, die Summe der geraden Zeilen Y,; zur Abkürzung ist „Sina 5 =s gesetzt. 
Tafel 1. 
T-6 F-4 P-2 Yo 73 Ir I 
0g dc dc m Öc = HE dd f 
3, Pr 7 ‚92 0 so Pa 2, IP 3: 00 
DE) 0%, 0% RE), 09 0% 
Po %”g_, 0’p_, op, Ö’ pe 
6s? — 38 — ” Fu 0 0 Ss? > 38 a 
ou 00° 00° O0" 09? 
; 0° « 4 z 0° f; Au 0° q, 
we ae 0 0 0 Fe IT 
0% REVE Ö V, ou’ 


Die Zahlenkoeffizienten in der »n-ten Zeile sind rechts von den Nullen: 


N n n(n-+1) n(n + )(n—+2) 
‚  ° .-4 ° 1-3:3 
oder 
(n — 1)! n! (n +1)! (n +2)! 


DB. e—- Br  e—- Diz! (n — 1)!3t' 


Die bei jeder Näherung hinzukommenden Glieder sind durch die senkrechten Säulen 
gegeben. Für die nullte Näherung ist y, zu nehmen, für die erste Näherung kommt 9, und 
die Ableitung davon hinzu, für die zweite Näherung 9_,, dann @,, g_, und gleichzeitig die 
darunter stehenden Ableitungen und so fort. Die in jeder horizontalen Zeile auftretenden 
Glieder lassen sich zusammenfassen. Die erste Zeile gibt 


D-RF, (= Reotg (dd: -ÄNVENEE (Id, 4, iA). 


ne a : r ls Ö ’ ' ' 
Für die nächste Zeile können wir schreiben w NR go P), da dieser Ausdruck sich wie 
.d ( 
folgt umformen läßt: 
ls Ö 


"2 R, 5 Flo) 


2 


iin 980 


3 = - 2a 7’ 


(,ntmn): 


' ( 
— > ii 7 om 


18,09 \'mn u % OYam 
T ' ‘ f > 77 ms - 
1 . AT Ö ) |— „sem 1 “ |— Ö ) 
Die Zusammenfassung der nächsten. Zeilen gibt: 
1 ( s\. 0° Liavy. 
+ B -rr- m (r ee Ben: Se > jr 
2!\r n 0: 10 >: Fo) si = " da 6062 6-: Fo). 


Der Beweis ist derselbe wie bei der vorherigen Zusammenfassung. 


Nun können wir als Näherungen gleich die zusammengefaßten Zeilen verwenden, d.h. als 
nullte Näherung 


Z ek Ar 
A - sina®,RF,(d, 
( 


als zweite Näherung 
er. A ae ran 
IR „sin a I,RF,(L) 


 __ 9 8 A B. Pie 
- nad, 1,7 ° 39 ol). 
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Für die 2»-te Näherung erhält man: 


mn u 2 
Br .; 1. av 0° AR R ne 
F* = —sina y er asm(Co 5 S-2''"&-am-n» Sam F, (£)) 
2m) a ® 4 (2 m)! 7 Ö 3? 2 2 (2 2 2 { 
m=] 
(3). 
a m m 
1 _ 1 g\em-ı d° 1 
. ’ ; q . . ı fr» 
E£* sın a v, 7 2m—1)! 7 R gem 60 626-2° Sem-» +-em-n» FO) 
m—1l 
Für die (2» +1)-te Näherung erhält man: 
mn : 
1 \’ 1 g\em am 
Fnsn=—sinad, _ ER DEE ET ER 
en 1) a o >) m)! A d gen 0>2 2 @am-—- 2») 2m ol 
wi 
(4). 
mn 
1 \ 1 s 2m-+ı +: 
—Sı — B RP: ER); r a Fr 
a sinad, Om-+Tmilr R gem Fil6o 626-2 Sam &-am Fılö)) 
wi . 


Bei jeder Näherung tritt jeweils ein Summenglied hinzu. 


Für 9,=0 und #,=0 wird am Rande bei jeder Näherung F*=0 (bis auf die singu- 
läre Stelle), da entweder sina® oder der Faktor dahinter gleich null wird. 


ı% %* 


N 
8% 
der (2n»—+1)-ten Näherung für 9, =0. Die Berechnung folgt, wobei Glieder, die unmittelbar 
null werden, fortgelassen sind: 


i a C u. Re . 
Die Ableitung 3» der 2n-ten Näherung wird für 9, =0 gleich null, ebenso 


mn 
OF‘, .. “ 1 Tg Po ude (£ 2 2 F 2 
RW) u 2 m)! 7, gm >o Ss-(2m-» sam 0 (5)) 
m—=1 


,—=V; 


u; Pr 1 g \m-ı am A 
+7 sin (a 2, (2m — ie R gem (o'''S-em-n» F, 0) 


/ 


mn 


S. l Ss v ar . . . 2) Pe 
ng: / et Rgeml&o 2 am »llam tma)F,()]. 


m=] 


Die Klammer vor F,(£) gibt aber £,; F,(£) wird mit einer geraden Funktion in £ 
multipliziert, so daß alle Glieder der Summe gleich null werden. 


0 


k ; ’ d Frun4 : . 
In gleicher Weise erfolgt der Beweis, daß für d, =0 - 3 gleich null wird. 


Auf der entgegengesetzten Seite des Streifens bleibt für die Ableitung nach # jeweils 
das zuletzt hinzugefügte Glied der Summe nach, da die anderen sich, wie eben bewiesen, 
wegheben. 


Also wird für 9, =0: 


d Ft, .. pn 1 g\rn gen+ı R . A 
na) Rgenri (lo: Can Feld) 


(öl 


1 1 s2n+1 d?r+1 


2 (2n)! a?” Bern 





ZT 
> 





)) (a en (3) )ies 





und für 9, —=0 


Fan+ 1 ; FL l s\ntı gan+2 BR 
a znlag)e Dil RN enrslle Lanka Fl) 
0} a 2 2n+1)! 7 g?n+2\0 en en#o 


1 1 s?n+t?2 d?n+2 
2» ı_.nf+ti Euren e 7°) ++ -(d,° —n?a®)etg£,]. 
2 @n+1)!n 0, z 


Zur späteren Berechnung der Spannung im Punkte M des Bildes 1 bzw. im Punkte M’ 


)2 F* 
des Bildes 2 benötigen wir IR für den Punkt 9,=0, i=0. Schließen wir mit der 


(2n+-1)-ten Näherung ab, so erhalten wir aus (4) unter Fortlassung der Glieder, die null werden; 
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mn 
> em 2m. 
1 ( Ss 2 0° 1 


22 E* ) 
c «ıntn)] 2 - m jr 
RR gem+illo Lam FO) 


2cosad, 











00° Hl, i=0 (2m)! zT 
m—U 
mz=Nn m N Ir 
y 1 s\? +1 d? +? | 
— »)E0HS < 4 SR RE 5 Fi | 
2cosa#d | RN m ll ER ld) 
ı  , @m+1)!\n, ee ac 9,0, a=0 
m—U 
m n . Pr £ 
„J ie SE ff Eu ud (»; {a \ r [2=\) (0, Kuh MM in 
= — zıc0Sa — 1 re Bat 4 (i ki )-- u a) )tgi 
| 2 —_ (@m)!\n 7 Kusäai 2 e 2 2 er 
m=i 
mn : R e 

B 1 u ee ((» r (=) (9, 2m—1 "ls 2m+I1 r y\\ 

mi: 5 asmr31Id, > |, mn z are. —— aitgı . 
dat Em+nilz 89,7 M a\l-ı > N 2 | 1 2 | 1 ls, j 
m—U 

\ ; z .. 20 

In der eckigen Klammer können wir 5 7 fortlassen, da die daraus entstehenden 

Glieder 9%, in ungerader Potenz enthalten. Den zu differentierenden Ausdruck fassen wir 
/ 2 
weiter als Produkt auf von 9, und ((»;: -(3 )- tg 9) und differentieren nach Fak- 
toren. Dann erhält man: 
mn 
u na nn l gi" 
nt] —=2 (2m + 2)—(2m+1l)cosaz: - ‚| 
09° Nh=0,i—=0 PR 2 s/)2m+1)!\r 
m—U 


(5). 
AEaL N 


a\® 2 [ \ a 
gm ‚(9 (5) |- j (9? | “ a) |tg 9) | 


IV. Bedeutung der Singularitäten und der Näherungen. Für die tangentiale Kraft von der 
Größe z nach Bild 5 erhält man den Spannungszustand 


3 2 2 
Fr xy x 
= FR 2 y ‚ Tey° 2 #, 
r | r x! r 


u 
0. = _ 


Lassen wir zwei entgegengesetzt gerichtete Kräfte z/e im Abstande e nach Bild 6 
wirken, so erhalten wir den Spannungszustand durch Differentiation nach x, also 
Ö ( x 22° —6x’y’ 
) 


06. = USW. 
x dx r® 


Dieser Spannungszustand stimmt überein mit dem Spannungszustande der Singularität 
im Punkte U’ des Ebenen-Ausschnittes: 


























0 y? 22-6’? 
a y 2\r’? L y 2) r'® USW. 
2 D 
r Aids) 7 
IABc a jan Bild 5. 
€ 
LLLLLLLLLLNNLLLLLLLLLL., 
| NE 
= 2a [Aiss$) Re A184.78.d 
Bild 4. Bild 6. Bild 7a und b. 


Nach der nullten Näherung erhalten wir einen Ebenen-Ausschnitt mit dieser singulären 


Belastung im Punkte U’ und einer tangentialen Randbelastung der Kante 9" — «a 5, die im 


Punkte r"—=1 das Vorzeichen wechselt, Bild 7a. Bei der nächsten Näherung wird ein 

Spannungszustand hinzugefügt, der diese Randbelastung beseitigt. Dafür tritt jedoch am 
n - - zz : 

Rande 9’ =0 eine tangentiale Randbelastung auf, Bild 7b. Die Randbelastungen werden 

mit jeder Näherung immer geringer und verschieben sich immer mehr vom Punkte r" —1. 
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‚ 


Die gesamte von der Ecke nach außen wirkende Tangentialkraft von r” —=0 bis r”—=1 wird 
nach der (2? rn + 1)-ten Näherung an der Kante 9” —=0 
1 ) 
- Kia u s 
ir \i T.,)dx \ 2. d.x' e R 
2 + er “ N ; m 
fie nr ‚er 0dy Oyleı 


(° E [0 (F” r")] 1 ? Br n+n 
09 Ie_ı, sr —o | 09 = en q 0% 0, 90 
1 1 1 s?r+?2 (+2 
- = 2 (9,192? —7?)---(M? "n)ctg W,)- 
a 2 (2 n +-1)'n?"t!9 FRURER o o ( u „ 7T ) tg U) 
a 2 r gr T 
Ebenso nach der 2n-ten Näherung an der Kante 9’ =a 5: 
me) 
[62737 Ri) 9” r— 1,9 a - 
N 1 | sr+tı gen+i ((o er j (» Z PEN 2 
=. ze a, ee 1 VRR FOREN EHEN  y FERN 
ei ae! er et . 2 ) "| 2 en 


Die Auswertung gibt: 


N > 
PEERENR,. C 5 s 
=gan,t8 v, Ia 
z ERENR 08, 2 ss er 
Ten 2a | a ) 1!09? u 205 Ban 2:08 
ee EEE}, 4 N Br nn 
ar) 5 (sehn, 30178 
-. s (s\ 
E 15 ) +4-0,114 
m S [ Ss \ - > 
& tz -5: 0,081 
m s (Ss \ . ö 
as == r> -65:0.061. 
Für a=1, also für die Viertelebene, erhält man: 7,=1, T,=—-0,425, T, = 0,220, 
T,=— 0,119, T,,=0,065, Ta, = — 0,037. Das Verhältnis der aufeinander folgenden Zahlen- 


r e 2 s San + 7T . R 
werte bei der Viertelebene nähert sich allmählich dem Werte —, wie folgende Überlegung 


zeigt. Die erste Näherung stellt den Spannungszustand dar nach Bild Sa, wo die Viertel- 
ebene zur Halbebene y >0O ergänzt wird; hierbei wird der hinzugekommene Rand y’=0 so 
belastet, daß auf der Linie «”—=0 nur Schubspannungen entstehen. Um diese Schub- 
spannungen zu beseitigen wird nach Bild Sb ein Spannungszustand hinzugefügt, der dadurch 
entsteht, daß man die Halbebene «’”’ >0 so belastet, daß jetzt auf der Linie y„’"=0 nur 
Schubspannungen entstehen; und so fort. 





_— Ai —_— 
Aids Ba. c 





Bild Sa. Bild Sb. Bild 8e. 


Y 


Vollführt man einen solchen Prozeß, wobei auf jeder Viertelebene nach Bild Se nur 
eine tangentiale Einzelkraft wirkt, so erhält man auf dem anderen Rande der \ iertelebene 
Schubspannungen nur nach einer Seite, deren Gesamtheit eine Tangentialkraft gibt, die das 
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) 

z fache der belastenden Einzelkraft ist. Dasselbe Verhältnis erhält man auch, wenn auf 

dem ersten Rande eine verteilte Tangentialkraft wirkt, wobei die Schubspannungen gleich- 

gerichtet sein müssen. Sind die Schubspannungen auf dem ersten Rande zum Teil nach 

außen, zum Teil nach innen gerichtet, so heben sich auf dem anderen Rande die Schub- 

> 

spannungen teilweise auf, so daß das Verhältnis kleiner als — wird; hierbei ist auf jedem 
” T 

Rande die Summe aller Absolutwerte der Kräfte zu nehmen. Teilt man die erzeugende 

u . “ . 

5 ist. Dieser Fall liegt 


bei der Untersuchung nach Bild Sa und b vor. Da der Hauptanteil der Kräfte bei jedem 
Schritte sich immer mehr nach der Ecke, bzw. ins Unendliche verschiebt, so nähert sich das 


Kraft durch die erzeugte, so erhält man einen Wert, der größer als 


r = . en . r T 
Verhältnis der Kräfte immer mehr dem Werte —. 


V. Spannung im Punkte WM, Bild 1. Im Punkte M’, Bild 2, tritt nur die Spannung o,r 
auf; diese ist 
0° lo Er er 


\ = 
3 > 3a Er). 
..: 3 Treat) 


Die beiden ersten Glieder des Differentiationsoperators werden gleich null; 0,» vereinfacht 

EEE 

aut 09 

nach einer ungeraden Näherung, so daß dann der Rand, auf dem M’ liegt, ohne Schub- 

belastung ist. Auf dem anderen Rande #9—=0 hingegen tritt eine Schubspannung auf, die 

von 2” —=0 bis «”=1 die tangentiale Kraft T,n+,», von @&”’=1 bis 2&’”—=» die Kraft 
Tan+n gibt. 


sich zu Diesen Ausdruck haben wir, siehe G]. (5), berechnet und zwar 


Die Spannung 0,” ist numerisch leicht zu berechnen; die Näherungswerte hierfür seien 
Os Os On, usw. Man erhält 


REN. 
cos az \ zsmaz 
1 aa 2 2 
Ou=-—12 2 
: a? A B na 
sin a5 
a u \3 
cos a zsına 
FAR. „na 2 2 aa: 
9,0, tr514-35- - 0,356 
; a” Re JT Ta 
sınayz 
na u 
coSa5Z ssnay 
@. 3 - —} .0,162 
o 0, +16 —D -V,162 
-. = = PL Ta 
sın a5 


und so.fort. 


Nun soll noch eine Verbesserung vorgenommen werden, welche die Schubkräfte am 
Rande 9’ =0 berücksichtigt. Wir bringen hierzu an der Ecke des Ebenen-Ausschnittes die 
Kraft — T,,„+,, an, so daß wir eine Belastung nach Bild 9 hinzu- 
fügen. Der Spannungszustand für eine solche Scheibenbelastung 


ist bekannt. Auf der Linie 9’ —=a-- erhält man in der Ent- 


[7 


fernung r I die radiale Normalspannung 





. 7 7 : 
snay —-az(0Say 
Tr m - Ö 
a\ ME 
Bild 9. (a 5) sın’ a5 


D) 


fo » . vr . . Jı E 
Ist o,» für Punkt M’ gefunden, so wird für M, Bild 1: 0x, = (2sina 5) or", wobei 


0%,» —2 war. Für a=1, also für die Viertelebene, wird die Verbesserung gleich 137 T. 
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Wir erhalten im Punkte M’ eine gesamte Spaunung, die für «—=1 und für verschiedene 
Näherungen weiter berechnet sind. Die Spannungen ohne — 1,37 T sind kleiner, mit — 1,37 7 
größer als die wirklichen; bei den höheren Näherungen nimmt der Fehler schnell ab, be- 
sonders mit der Verbesserung: 


9 — 1,37 T.., = 2,546 + 0,582 = 3,128 
0% — 1,37 Ta, = 2,915 + 0,162 = 3,077 
6, — 1,37 T,, = 3,018 + 0,052 = 3,070. 


Im Punkte M im Grunde der Halbkreiskerbe wird mit der letzten Näherung o,, y = 6,14 
bei ox,,&—=2. Die Spannungserhöhung ist eine 3,07-fache. Zum Vergleich: Bei der Vollebene 
mit Kreisloch ist die Spannungserhöhung eine 3-fache. 


Für «=0, s—=n, Vollkreiskerbe nach Bild 4, erhält man: 


® L: 
I„=-.52088; 1T,=-, 54.0145 17,=—-—5:6:0061. 


r . 3 ‚2 
Die Spannung durch — T,,„+, wird durch Grenzübergang — > BERR. 


Für die Spannung in M’ erhalten wir für die einzelnen Näherungen: 


4 T, 1 h 


n a ) 


2,000 + 1,333 = 3,333 


9 
a” (0. 


4 To,\ x 
a? (0. F” . ) — 2,356 + 0,912 —= 3,268 


4 To, 


T a 


I 


\= 2,518 + 0,732 = 3,350. 


. Yo . Pi} - IT y h} - 
Mit dem letzten Näherungswert wird a’ oyr = 3,25 und o,,y = (5) :3,25—=8,0 gegen o,,. —2. 


Statt des Grenzüberganges «>0 kann man diesen Fall auch unmittelbar lösen. Die Zahlen- 
rechnung bleibt genau dieselbe. Vergleichen wir dieses 
[ Ergebnis mit der Spannungserhöhung der Vollebene 
RE mit elliptischem Loch. Hierbei soll der Krümmungs- 


y, u) halbmesser am Ende der langen Achse gleich der 
8 
Dr) 





Hälfte der langen Halbachse a sein und die Ebene in 

Richtung der kurzen Achse gezogen werden, Bild 10. 
%,»-2 Das Achsverhältnis ist dann gleich Y2 und die größte 
ER | Spannung im Punkte M wird, wenn 0, —=2 ist: 








. 5 L — f D) ” —2(1 2YV2)—2.38 

BR 7 ou=0.(1427)=2(14+2y2)=2-38. 
:w 

em Die Spannungserhöhung wird eine 3,8-fache; dies 
Bild 10. ist ein recht brauchbarer Näherungswert. 




















VI. Schlußbemerkungen. Wichtig erscheinen mir nicht die zahlenmäßigen Ergebnisse, 
sondern die durch diese Untersuchung erschlossene Möglichkeit, weitere Scheibenprobleme 
zu behandeln. Selbstverständlich kann man die 9,/ als angewandte Koordinaten im ur- 
sprünglichen Gebiete einführen. Die Anschaulichkeit wird hierdurch verringert und die Ent- 
stehung des Lösungsansatzes unverständlich. Auch die Inversion erweist sich von Vorteil, 
da hierdurch die Verbesserung der Näherungslösungen durch die Kräfte — T,n+, möglich 
geworden ist. 184 
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The Stress Distribution in a Semi-infinite Domain Having 
a Plane Boundary and Compressed by a Rigid Body. 
By Hajimu Ökubo in Yonezawa (Japam). 


Zusammenfassung. Es wird das zweidimensionale Problem des elastischen Körpers be- 
handelt, auf den ein starrer Körper beliebiger Form einen Druck ausübt. Ist ersterer groß 
gegen den starren Körper, kann man dessen Begrenzung angenähert durch eine Ebene er- 
setzen. An der Berührungsfläche werden nicht wie sonst üblich die Drücke, sondern die 
Verschiebungen vorgeschrieben, wie sie der Oberflächenform des starren Körpers entsprechen. 
Dabei wird zur Vereinfachung der Rechnung angenommen, daß in der Berührungsfläche kein 
Gleiten eintreten kann. Die dadurch entstehenden Scherspannungen sind, wie die Rechnung 
ergibt, klein gegen die Normalspannungen. — Zahlenmäßig wird der Fall durchgerechnet, daß 
auf einen Körper aus Gußeisen von einem starren Körper mit ebener, 2 mm breiter Fläche 
ein Druck ausgeübt wird. 


$ 1. Introduction. We will treat in this paper, the two-dimensional problem of the 
equilibrium of an elastic solid body compressed by a rigid body of any shape. Such a problem 
has hitherto been treated by using potential theory, under the assumption that the shear 
vanishes on the contact surface. And the normal tension on the contact surface is determined 
by giving the normal displacement of the contact surface without any consideration of the 
tangential displacement. But since the shape of the contact surface after deformation is not 
determined by the normal displacement exclusively, we will treat the problem in another way. 

If the linear dimensions of the elastice body are large compared with those of a rigid 
body, we can assume it to be bounded by a plane. In the case where the boundary is the 
axis of «, and the axis y penetrates into the body, in the ordinary way the normal stress Y,, 
at the contact surface is assumed and the problem is solved under the corresponding boun- 
dary conditions. But such an assumption does not generally have a theoretical basis, and 
the results are not capable of exact solution. 

So we will not assume the stress distribution at the contact surface but will give the 
displacements which are determined by the magnitude of the force operated and by the shape 
of the rigid body. 

To simplify the caleulation it may be assumed that between the contact parts of two 
bodies, a slip does not oceur and therefore the displacement of the x direction vanishes at 
the part. As will be seen in the later example, at the contact surface the shearing stress is 
small compared with the normal stress, so the above assumption will hold generally in 
practice. 


$ 2. Equations of displacements and stress components. When 
the elastie body is compressed by a rigid body with a vertical force 
P and the two bodies are in contact in the interval —e<Sr<se, 


f 
as is shown in Fig. 1, the boundary conditions are given by displace- KH 
ments at —e<x=c and by surface tractions at z >c and z=< —.c. Z x 





For the caleulation, it is not eonvenient to use the cartesian 





. ” . * — 
coordinates under such boundary conditions, so we will take the 2\c 
ceurvilinear coordinates defined by TER 

5 J 

z=cCosw, where z=er-+tiy and w=a-tiß Fig. ı. 
or 

= .cÜosacosß 

EN (N, 

y=eSinasinß 

where 
. eX—e-* ö ex e—* 
Cose = = . Sn - 


If we denote by « and v the displacements in x and y directions respeetively, then «, 
v can be expressed in the forms’) 


ı) H. Love: The Mathematical Theory of Elastieity, 4th, edition, p. 205. 
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E A+u 0 x 
u: + of N + u, 
2 u 2ul+2u)Fdy 
1 (2 
7 it+u 07 er 
V = = gr ! " 
2(Ä+2u) 2u(@+2 u) 9 Yy 
where E+ in and ®’ + iu’ are two funetions of z. 
By a slight calculation, we have 
vr +Hiu=il+3 NP, + WM yyY,(2)+ 9; (2) (3), 


where , and g, are two arbitrary functions of z and $} denotes the imaginary part of a 
complex funetion divided by i. 

If us and «,„ are the tangential and the normal displacements with respect to the 
curve a= constant, then we have 


N\z 
R 0 ö 
ustiu, hr, eriW. ae Er EEE a © 


-where 
l (/0 2 oy ® "a ; 
+| : e YSin’a-+ sin? pP. 
h | da) de * P 
Since y, and g, are functions of z, they are also function of w and if we denote them 
by f, and f, respectively, from equations (3) and (4) we have 
(A+u)y 
ec Sin w 


Ö a 


Sin w | 
y Sin? a + sin? $ | 


HOESAUE 7: 


ustiu,: +30) Sf (w)] 


Now the strain components are’) 


EB. 6 
mus) arm). (6), 


0 do /1\ 3 Oo /1 

OU, ( Üu3 r ( 
ee +h?’uss; ) e33—=h +h’u I; ), Ca3 = 
=. da Ö Al h e ey u “Yalh 3 0 





| N 
or 
BOOTE |. ’ 
ET u |(Sin‘ a+ sin? pP) I - +sıinpcospu N A 
de 
TE PURE EBEN ME: \ . 
e3—=ch Sin’«a + sin? f) Ta + Sın a Cosa Ha I 
ep 
TE [93 , du ” i 
e.3—e?h? Sin’ a + sin? pl + a) — Sin a Cosa us — sin cos u 
\ \da oß J 


and the stress components are 


aa=(i+2u)euat fe33, pB=lÄ+2u) es3+ Aegas afzuas.. . © 

$ 3. One type of solution. The vertical displacement of the contact surface is deter- 
mined by the magnitude of the resultant of the surface tractions at the part and by the 
shape of the rigid body and is expanded in general in Fourier’s cosine series in the interval 
0=/#=7. Hence the boundary conditions are 


.) 


I’ —0 at P=0, N 


BB —y 


— 


Y 


Ir; 


U„eos2nP, us), BR. rt 


Now putting 


® 2 
Klw)=i 3 A,e?"”+ SD B,e er tDw 
n—1 n—\ 
© e2 . . . . . R e (11), 
f,(w)= 2’ C„e RR ı Eid tw 


n=]1 n—0 
where A,„, B„, C„ and D,„ are real constants, to satisfy the boundary conditions (9), from 
equations (5), (6), (7) and (8), we have 


2) H. Love: loe. eit. p. 54. 
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a4 > hi) u3 5 a Ö u 
(A+ 2 (Sin a dB +Cosa u.) +4Sin a 3 a —0, | 
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En e > 2 )@ 5) 5) „- Zn+1)a 
öß (1 Be Due 2u > (@en+1)B,e , 
#8 n—( 
a 
7) I _0 : 
ı=— Y @n+1)D,er+De 
0 a — 
n—(0 
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oO us \ I Pr ® FR 
> =- (+3 Y 2n- Ane """— \ ET 
ea — — 
n—=1 n—] 
Hence the constants must satisfy the following equations 
ui+2W2n +1) B,„+tni—-W)D,„=0 .. et A 
2uA,-+Cn: N ee Eee ee 
By putting «=0 in equation (5), the displacements at the contact surface are 
E a 
tu Pat +3 u) A„+ C„)cos?n RAT, +3 u) B„—D,)sin@n +1), | 
n—=]1 n—=0 
(15) 
ee) - 
= )C, sin2nß— }D, cos(2n+1)Pß | 
„=ıi n—0 
$ 4. Another type of solution. In this article, we shall obtain I Kir 
the stress components and displacements in a semi-infinite domain * 
having a plane boundary and compressed by a force P which is 2Ic 
distributed with the uniform pressure on the interval -—e<x<se of 
the boundary as is shown in Fig. 2. Such a case is readily solved 
by a slight calculation, we have’) 2 
Fig. 2. x 
Y 
e P j (e+a)y (e—x)y c+x c—x| 
.=sneletartp C-a+ arctg y — arctg y 
P (+2)y (e—x)y c—x c+x| 
Y,=- | . -: arctg arc 6). 
y Irelicta®+y (ap yrar tg y +arctg yf (16) 
win. Py’ l. 1 1 
I 2nelce+2®’+yP (—-a’+Yy 
3) H. Love: loe. eit. p. 192. 
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The displacements of the case are 


Fr (i+2u) (+2) + y° c+x c— | 
u = naar 3 ylog,, Zr etranig y + ule — x) arctg ri 
P i+2u) & ” (A+2u) a. R en 
v= "Tau (ce+ x) log ((e+ x)?+ ?) + > (e— z)log((e— a)” +y’) (17). 
c+x c | 
+ u yarctg uyarctg 
uy y +uy y| 





At the contact surface, we have 
P 


u aa Terra 
(+2 u) P önk 
4 su 
ce FE ‘ ’ » o ir -r » y .\‘ 
v= Incnürmte trakg x) + (ce — x) log (ce — x)} 
or at a=0, they are 
u = "FürmeeP: 
IL (19). 
+2 R 
v Er ((1+ cos A) log (1 + cos P) + (1 cos B) log (1 — cos ß) +2 log c} 


Anuli+u) 


$ 5. Solution satisfying the boundary conditions. The solutions obtained in artieles 3 
and 4 both satisfy the boundary condition (9), so it is evident that the sum of the two 
solutions also satisfies the boundary condition (9). And since, on the contact surface, the 
resultant of the normal tension aa given in article 3 vanishes *), the sum of the two solutions 
is the solution of the case where the elastic body is compressed by the vertical force P. 

Next, we will determine the arbitrary constants A,, B„, C„ and D,„ to satisfy the 
boundary condition (10). 


At a=0, the displacements are 








,+2u) P 
ir tl + cos P)log (1 +e0sß)+( - cos B) log (1 — cos #)+2log c} 
u Y ((A+3u) A„+ C,) cos2n pP — \ (+3 u) B„— D,)sin2n+1)P, 
n=1 n—( 
»=— cosß + \'o sin2nß— Sn cos(2n+1)Pß 
"7 Aätu) TA ec: ES er 
n=1 n = 
In the interval O<P=<r, we can put 
2 4 = 2m+1 
 ., n N u.‘ 
ART: Amt+la’a > 2m+1) BEER P 
n—=] 

Aue Weg n . 2 
oe (U EEE } u Hi ) Va | 
cs? m+1V)ß= 2. Em-+ir_gmsin2nP, 

n—=1 
(1-+cos ß) log (1 + cosß) + (1 - cos ß) log (1 — cos $) +2 loge - Bi b„cos2nß 
n—=Uü 
where 


b,= 0.615 +2loge, b,—0.667, b,—0.067, 5,—0.0%0, b,—=0.07, b,— 0.005, 
b,=0.0M1, b,—=0.003, b,=0.01, b,=— 0.01, d,=0.02, b,,= — 0.001 °). 


10 


ü a 2 i an r ? 00% 5 

4) Consider a semi-eirele with its centre at the origin, then the radius of the eirele is 5 when « is taken large 

enough, while the stresses given by (8) are the infinitesimals of the order e”2«@, the resultant of the stresses on the 
— 

semi-eirele vanishes when the radius becomes infinity. Therefore the resultant of the normal tension @« on the eontaet 
surface must vanish. 

5) b„ taken here is caleulated by Zipperer’s method (see W. Hort: Technische Schwingungslehre, 2. Aufl., 
S. 136) and in this case gives the approximate value with suffiecient aceuraey. 
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To satisfy the boundary condition (10), from the equation (20) we have the following 





relations 
ee 
Bi 2nP 8n \’ Ds = 
"Taitu)dn?—1) nn ku Zm+1) tn? (22), 
m—U( 
N 2 |ranDm _ A+2WbP 
r.. wc St taruli+ u) 
m v 
(23), 
rt S_mem+ „ _iA+n, _ Ü+2Mb„P . 
Be ir 2u " Aruliu) 
m—\ 
where 
A433) (mi — u) N 


Un = —- E= 
" ulk+2W2m-+l1) 


From the equations (13), (14), (22) and (23) we can determine the constants A,„, B„, ©, 
and D,„ by solving the four systems of simultaneous equations. Inserting them in equations 
(5) and (8) and superposing the values given by equations (19) and (18) respectively, we can 
know the displacements and the stress components at any point in the elastic body. 


$ 6. Numerical calculations. As an example, we will take a body of cast iron which 
is compressed by a rigid body with a flat face whose width is 2 mm with a vertical force 
of P kg weight. The displacements and the stress components on the contact surface are 
caleulated as follows: 


Taking the values of the elastie constants as A=8,664, u=8,204 (unit kg weight 
per mm?)°®), the simultaneous equations to determine the constants A„, B„, C„ and D, are 


ı — 0.947 








” - Y . . Y 
An = — 0.609 %xX 10". C,, B„ = — 0.4212 X 10" X wa DD; 
n P \ s Du 
C,,= 0.3177 0". — 2,547 - -r 
was Tax ı 4n?—1 . 2m+1)?’— 4n? 
y m 
“> 3.402 m + 2.056 n — 0.327 - EIN: n P 
h - D.. = 0.1132 0*.t: > +0.30: 0. — 2 
ee 2m+1®? —4n? Din 1182 X 1 nP+ xı 4n?—1 
m 
where m =0,1,2, 10. 1:53.38, 11, 
Table 1. 
n ® Bu Cu Du 
0 — 0.974 x 10°? P 2.441 X 10° P 
1 0.596 x 10 P 0.002 0.978 X 10° P 0.288 
2 0.114 0.009 — 0.187 0.102 
3 0.060 - 0.007 0.098 0.054 
4 0.038 0.004 0.063 0.031 
5 0.029 — 0.003 0.047 0.022 
6 0.023 0.002 0.038 0.014 
7 0.019 - 0.002 — 0.031 0.010 
8 0.015 — 0.001 0.025 0.008 
9 0.014 0.001 0.023 0.003 
10 0.010 — 0,001 0.017 0.004 
11 0.014 — 0.023 


*) Landolt: Physikalisch-Chemische Tabellen. 
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Since in this case UU=U,=-U,=..... =l„=0. The vertical displacement of the 
‚ase IS 
\ '2.166 m — 0.205 
1. 2m+1) 


m 


U,= 


Din — 0.887 X 10° P. 


The values of the constants caleulated from the above simultaneous equations are 
tabulated in table 1. 


In this case, the vertical displacement of the' contact surface is 
UV =—-0.285xX10°P. 


To compare the displacement of the contact surface in the case treated in article 4 
with the result obtained here, the displacements of both cases are plotted in Fig. 3. 
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Next, we will caleulate the stress components on the contact surface. From equations 
(7) and (8), at «—=0, we have 
an 
ıP P 
4 +u) 2’ 


A+2u |’. s EEE RR 
sinB u (Zn+1l)(2uB„— D,)sin(2n+1)ß- 


n—U 


aa 


—  ulA+3u) . i ı S ? 
aß= sin B N 2n Ans 2 tn >) (Zn+1l)f(Aä+u) B„+ D„} cos (2n-+1)P 


n—=U == 


(+2W)P, Z. a 3 
"ia r+ mtl +cosp) log (1 — cos p)} 


The values of stress components calculated from the above equation are plotted in Fig. 4°). 


(24). 





$ 7. Concelusion. The vertical displacement of the contact surface is not uniform and 
the horizontal displacement does not vanish, when the body is compressed with uniform 
pressure, as is shown in Fig. 3. 

Likewise, the normal traction on the contact surface is not uniform and the shear does 
not vanish, when the contact surface is displaced vertically with uniform magnitude, as is 
shown in Fig. 4 Thus, we can not assume the stress distribution on the eontact surface 
even in approximation simply by the shape of the rigid body such as is used in ordinary 
practice. 

In conelusion, we wish to express our thanks to Prof. Iwao Kobayashi and to 
Prof. Seiichi Higuchi of the Töhoku Imperial University for their kind guidance in 
this study. 204 





?) A similar case is treated by M. Sadowsky under the assumption that the shearing stress on the contaet 
surface does not exist (see Z. angew. Math. Mech., Bd. $ (1928), S. 107). 

But such an assumption is far from the actual state, since the frietion between the two bodies gives the 
shearing stress on the contact surface. In so far as the distribution of the normal tension on the eontaet surface is 


eoncerned, the result Yy=— v >, obtained by him is rather near in value to the result obtained here under the 
ne x? 


more actual boundary condition. 
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Tafel der acht ersten Kugelfunktionen zweiter Art. 
Von F. Vandrey in Göttingen. 
Bericht der Aerodynamischen Versuchsanstalt Göttingen E. V., Abteilung für theoretische Aerodynamik. 


ie folgenden Tafeln enthalten die Werte der Kugelfunktionen zweiter Art von Q, bis Q, im 
Intervall O<x<1. Für 1<x=0 ergeben sich die Werte dieser Funktionen nach 
der Beziehung 
9, -- 2) =(- 1)" r!Q,(&). 


Zur Berechnung der Tafeln wurden die Werte der Funktion 
9, (2) = Artgx 
neunstellig der Tafel von K. Hayashi') entnommen. Die Berechnung der weiteren Funktions- 


werte erfolgte dann für @, bis Q, direkt mit Hilfe der bekannten geschlossenen Ausdrücke 
für diese Funktionen), für Q, bis 9, mit Hilfe der Rekursionsformel der Kugelfunktionen 


n dn+(n — Dan-.—- (2nr -YD)x n-ı=0. 
In den Tafeln bedeutet ein Punkt oben an der letzten Ziffer, daß der betreffende Wert 


erhöht wurde. Lineare Interpolation ist bei Beschränkung auf drei Stellen möglich im Intervall 
von O0 bis ungefähr 0,86 bei Q, bis Q, und von O bis 0,69 bei Q.. 











E % Qı Q: Q; Q Q; % Q: 
0,00 ' —+- 0,00000 1.0000 0.00000 | —+ 0,66667' 0,00000 0,53333 0.00000 | 1 0.45714 
01 | + 0,01000  — 0,99990' 0,02000° | —- 0,66627° | + 0,02666' 0,53253 0,03198° | —+- 0,45586 
02 | + 0,02000 — 0,99960° 0.03999° | —+- 0,66507 | —+- 0,05327 0,53014° 0,06383° | —+- 0,45203 
03 | —0,03001 0,99910' 0,05996 | —0,66307 | —+- 0,07978 0,52615° | — 0,09542 | —+ 0,44567 
04 | + 0,04002 0,99840" 0.07991 | + 0,66027 | —+-0,10615 0.52057 0,12664° | -+- 0,43680' 
0,05 |, + 0,05004 | — 0,99750° | — 0,09983 | + 0,65668° | + 0,13233 0,51343 0,15734 | + 0,42548° 
06 | + 0,06007 0.99640' 0,11971 | + 0,65229 | + 0,15827 0,50474 0,18742° | + 0,41175 
07 | —+0,07011 0,99509 0,13954 | + 0,64711 —+- 0,18393' 0,49452' 0.21674° | + 0.39570° 
08 | —+ 0,08017 0.99359" 0,15932° | + 0,64115° | + 0,20925° 0,48279 0.24518 + 0,37739 
09 | + 0,09024 0.99188° 0.17903° | + 0,63440° | + 0,23 119° 0,46958 0,27264° | + 0,35693° 
0,10 | + 0,10934° 0,98997' 0,19866 | + 0,62687° | + 0,25870' 0,45493° 0,29899° | + 0,33441 
11 | + 0,11045° 0,98785 0,21822° | + 0.61856 | + 0,28274' 0,43887 0,32412° | + 0,30996° 
12 | + 0,12058 — 0,98553 0,23769° | + 0,60948 | —- 0,30626° 0,42144° - 0,34793° | + 0,28369 
13 | + 0.130744 0,98300 0.25706° | + 0,59964 H 0,32921' 0,40268° | — 0,37031 + 0,25575° 
14 | + 0,14093° | — 0,98027 0.27632° | + 0,58904° | + 0,35155 0.38264° | — 0,39117 | + 0,22627 
0.15 | + 0,15114 0,97733° 0,29547° | + 0,57769° | + 0,37324 0,36137 0,41041 | + 0,19542° 
16 | + 0,16139° 0,97418° | — 0,31450° | + 0,56559° | + 0,39424° 0,33893° 0,42795° | —+ 0,16335' 
17 | + 0,17167° - 0,97082° | — 0,33339 | + 0,55275 | + 0.41449' 0,31537' 0,44370° | —+- 0,13023 
18 | + 0,18198 0,96724 0,35215° | + 0,53918 | + 0,43395 0,29075° | — 0,45757 | —+ 0.09625 
19 | + 0,19234° 0,96346° 0,37075 | + 0,52490° | + 0,45259 | — 0,26513 | — 0,46952° | + 0,06158 
0,20 | + 0,20273 0.95945 0.38920 | + 0,50990 | + 0,47037° | — 0,23859' 0,47946° | + 0.02642 
21 |: + 0,21317 0,95523 | — 0,40748 | + 0,49420 | + 0,48723 0,21119° | — 0,48734° | — 0,00904 
22 | +0,22366° | — 0,95080° 0,42559 | + 0,47781 | + 0,50315 0,18300 0,49310 0,04461° 
23 | + 0,23419' 0,94614° | — 0.44351 | + 0,46074 | —+- 0,51808 0,15411° | — 0,49672° 0,08008° 
24 + 0,24476 0,94126° | — 0.46123 | + 0,44301 | + 0,53199 - 0,12459' 0,49814 0,11524 
0,25 + 0,25541 0,93615° 0.47876 | + 0,42461 | + 0,54484' 0,09451 - 0,49735 0,14990 


1) K. Hayashi: Sieben- und mehrstellige Tafeln der Kreis- und Hyperbelfunktionen. Berlin 1926. 
2) Vgl. z.B.: E. Jahnke u. F. Emde: Funktionentafeln. 2. Aufl. Leipzig und Berlin 1933, S. 176. 
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ZT #7 Q, Q; Q; Q, Q; u Q; 
0,25 | + 0,25541 | — 0,93615 0,47876 | + 0,42461 | —+ 0,54484 — 0,09451 ı — 0,49735 | — 0,14990 
26 ' —+0.,26611° | — 0,93081 | — 0,49607 | + 0,40558° + 0,55659 0.06398° , — 0,49432 0,18385° 
27 | + 0,27686 | — 0,92525° 0,51316° | + 0.38591 | —+ 0,56721 - 0.03306 - 0,48904 0,21688° 
28 | + 0,28768 | — 0,91945° | — 0,53001° | + 0,36563 | + 0,57667° | — 0,00186 — 0,48151  — 0,24879' 
29 | + 0,29857° | — 0,91342° | — 0,54662° | + 0,34474 | + 0,58492 | —+-0,02953 | — 0,47173 0,27938° 
0,30 | + 0,30952° | — 0,90714 | — 0,56297 | ++ 0,32328° | +0,59195 | + 0,06103 0,5972  — 0.30845° 
31 | + 0,32055° | — 0.900638 ' — 0,57907° | + 0,30124° : + 0,59772 | + 0,09254° 0,44551°  — 0,33580 
32 | + 0,33165' 0,89387 0,59488 | + 0,27864 | + 0,60220 | +0,12395 | — 0,42912° 0,36126 
33 | + 0,34283° | — 0,88687°  — 0,61041 | + 0,25552° | + 0,60537 | + 0,15518' 0,41059 0,38464 
34 | + 0,35409 | — 0,87961' | — 0,62565° | + 0,23187 | + 0,60720° | + 0,18611° | — 0,38999 | — 0,40577 
0,35 | + 0,36544 0,87209 0,64057 | + 0.20773° | + 0,60766 | + 0.21664 | — 0,36737 0.42449° 
36 | + 0,37689' 0,86432 , — 0,65518° | + 0,18311° | + 0,60674 | + 0,24668 0,34281° | — 0,44063 
37 | + 0,38842 | — 0,85628 | — 0,66945° | + 0,15803° | + 0,60441 | —+ 0,27611 0,31638° | — 0.45407° 
38 | + 0,40006° ı — 0,84798° | — 0,68338° | + 0,13251° | —+- 0,60065 | -—+- 0,30484' - 0. 28818° | — 0,46466° 
39 | + 0,41180 | — 0,83940° | — 0,69695° | + 0,10658 | + 0,59545 | + 0.33274 | — 0,25830 | — 0.47229 
0,40 | -+ 0,42865° | — 0,83054 0,71015° | + 0.080%6 | 4 0,58879 | + 0,35972 | — 0,22687' 0,47686 
41 | + 0,43561 | — 0,82140° | — 0,72297° | + 0,05357 | + 0,58066 | —+ 0,38567 0,19399° 0.47828 


42 | +0,44769 | — 0,81197° | — 0,73539° | + 0,02654 | + 0,57105° | + 0,41048° | — 0,15981°) — 0.47649° 
43 | +0,45990° | — 0,80224 | — 0,74740° | — 0,00080 | -+0,55994 | + 0.43404° | — 0,12445 0,47142° 
44 | + 0,47223 | — 0,79222° | — 0,75898° | —0,02844° | + 0,54734° | -+ 0,45624 | — 0,08808° | — 0,46304' 


0,45 | + 0,48470 | — 0,78188 0,77012 | — 0,05634° | + 0,53323° | + 0,47698 | — 0,05085° | — 0,45133 
46 | + 0,49731 | — 0,77124° 0,78081° | — 0,08446 | +0,51761 | + 0,49615 | — 0,01292 | — 0,43631 
47 | +0,51007 | — 0,76027 | — 0,79102 | — 0,11279 | +0,50050° | +0,51365 | + 0,02552° | — 0,41800 
48 | + 0,52298 | — 0,74897° | — 0,80075' 0.14129° | + 0,48188 | + 0,52937 | + 0.06428 0,39645' 
49 | + 0,53606 0,73733 | — 0,80997 | — 0,16992 | + 0,46177 | + 0,54322° | + 0,10318 | — 0,37172° 

0,50 | + 0,54931° | — 0,72535° | — 0,81866 | — 0,19865 | + 0,44017 | + 0,55508 | + 0,14201 0,34392° 


51 | + 0,56273° | — 0,71301° | — 0,82682° | — 0,22745 | +0,41711° | + 0,56487° | + 0,18056 | — 0,31312° 
52 | + 0,57634° | — 0,70030 | — 0,83441° | — 0,25628 | + 0,39259° | + 0,57249° | + 0,21862° | — 0,27958 


53 | +0,59015° | — 0,68722 | — 0,84141 | — 0,28510 | + 0,36663° | + 0,57785° + 0,25595° | — 0,24337° 
54 | + 0,60416° | — 0,67376° | — 0,84782 | — 0,31387° | + 0.33926 | + 0,58086° | + 0,29233° | — 0,20471 
0,55 | +0,61838 | — 0,65989 | — 0,85360 | — 0,34254° | + 0,31051° | + 0,58143 | + 0,32752 0,16383° 
56 | + 0,63283 | — 0,64561 | — 0,85873 | — 0,37107 | + 0,28040° | + 0,57950° | + 0,36129° | — 0,12097 
57 | +0,64752 | — 0,63091 | — 0,86319 | — 0,39942 | +0,24897 | + 0,57498 | + 0,39338 0,07642° 
58 | +0,66246 | — 0,61577 | — 0,86695 | — 0,42754° | + 0,216%6 | + 9,56781° | + 0,42355 | — 0,03047 


+ 0,1823 | + 0,55792 | + 0,45155 | + 0,01655 


59 | + 0,67767° | — 0,60018° | — 0,86999° | — 0,45537 


0;60 | + 0,69315° | — 0,58411 | — 0,87227 | — 0,48287° | +0,14720° | + 0,5456 | + 0,47713° | + 0,06429° 
61 | + 0,70892 | — 0,56756° | — 0,87378° | — 0,50997° | + 0,11094 | + 0,52979° | + 0,50008 | + 0,11236 
62 | + 0,72501° | — 0,55050° | — 0,87446 | — 0,53662° | + 0,07362 | + 0,51145 | -+0,52000 | -+ 0.160836’ 
63 | + 0,7414 | — 0,53291° | — 0,87431' 0,56275° | + 0,03530° | + 0,49023° | + 0,53680° | + 0,20786° 
64 | + 0,75817 | — 0,51477° | — 0,87326 | — 0,58830 | —0,00395 | + 0,46609 | + 0,55017 | + 0,25441 


0,65 | + 0,77530° 0,49606° | — 0,87130 | — 0,61321° | — 0,04405° | + 0,43903 | + 0,55989° | + 0,29955° 
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t % Q, Q; Q; Q Q%; % Q: 
0,65 ı + 0,77530°  — 0,49606° | — 0,87130 | — 0,61321° | — 0,04405° | + 0,43903 | + 0,55989°  —- 0,29955° 
66 | + 0,79281 | — 0,47674 0,86838 | — 0,63739 | — 0,08490 | + 0,40905 | + 0,56570 | + 0,34277 
67 | —+-0,81074 | — 0,45680 | — 0,86446° | — 0,66078' 0,12642° | + 0,37616 | + 0,56740° | + 0,38358 
68 | + 0,82911  — 0,43620 | — 0,85948 | — 0,68328° | — 0,16849 | + 0,34039 | + 0,56476 | + 0,42145 
69 | + 0,8476  — 0,41491 | — 0,85341 0,70481 | — 0.21101° | + 0,30178 | + 0,55759 | —- 0.45585' 
0,70 |; + 0,86730 0,39289° | — 0,84618 0.72529° | — 0,25384 | + 0,26039 | + 0,54570 | + 0,48622 
71 +0,88718  — 0,37010' 0,83775° | — 0,74460 0,29686° | —- 0,21630° | + 0,52893° + 0,51203 
72 | + 0,90764 | — 0,34650° | — 0,82804' 0,76265' 0.33991° | + 0.16960° | + 0,50713° | + 0,53273 
73 | + 0,92873' 0,32203° | — 0,81699° 0,77931 — 0,38283 | + 0,12041° | + 0,48017 | + 0,54777 
74 | +0,95048° | — 0,29665° | — 0,80452' 0,79447 0,42546° | + 0,06887 | + 0,44798 | + 0,55662 
0,75 ' + 0,97296° | — 0,27028 | — 0,79055° | — 0,80799 | — 0,46758 | + 0,01516° | + 0,41050° | + 0,55877 
76 | + 0,99622° ° — 0,24288° | — 0,77499° | — 0,81973 | — 0,50900 | — 0,04053 | + 0,36770° | —+ 0,55372' 
77, + 1.02033° 0,21435° | — 0,75774° 0,82953' 0.54949° 0,09797° | + 0,31961° + 0,54101 
78 | +1,04537  — 0,18461 0,73868 | — 0,83721 0.,58878 0,15688 | + 0.26631 | —+ 0,52024 
79 | + 1,07143 | — 0,15357° | — 0,71770° 0,84259 0,62660 | — 0,21696 | —+ 0,20794° | —+ 0,49104 
0,80 | + 1,09861 | — 0,12111 0,69464° | — 0,84544 0,66264 0,27785 | + 0,14469° —+- 0,45312 
81 | + 1,12703° | — 0,0871 - 0,66935°  — 0,84555 | — 0,69656° 0,33914° | + 0,07684 | —+ 0,40628 
82 | + 1,15682' 0.05141° | — 0,64134 0,84264° 0,72795 0,40035° | —+ 0,00477 | + 0,35042 
83 | + 1,18814° | — 0,01385' 0,61132 | — 0,83641 0,75641° | — 0,46094 | — 0,07106 | + 0,28556' 
84 | + 1,22117 | + 0,02579° | — 0,57810° | — 0,82653' 0,78112 - 0,52029° | — 0,15006° | + 0,21187 
0,85 | + 1,25615 | + 0,06773° | — 0,54172 | — 0,81259 | — 0,80244° | — 0,57766' 0,23149° | + 0,12972° 
86 | + 1.29334 | + 0,11228° | — 0,50184° 0,79415° 0,81882° | — 0,63221° | — 0,31444° | + 0,03969 
87 + 1,33308° | + 0,15978° | — 0,45808 | — 0,77066° | — 0,82981°  — 0,68295 0,39780 0,05735° 
88 + 1,37577 | + 0,21068° | — 0,40979 | — 0,74148° | — 0,83453 0,72872° | — 0,48022 0,16020 
89  -+-1.42193° | + 0,26551 | — 0,35650 | — 0,70582 0,83194° | — 0,76811° | — 0,56002' 0,26725' 
0,90 | + 1,47222° | + 0,32500° | — 0,29736 | — 0,66271° 0,82075° | — 0,79944° | — 0,63512 0,37633° 
91 | + 1,52752 | + 0,39005' 0,23135° 0,61091° 0,79936 | — 0,82063' 0,70295° | — 0,48458 
92 | + 1,58903° | + 0.46190 | — 0,15708 0,54880 | — 0,76575 | — 0,82905° | — 0,76020 | — 0,58824 
93 | + 1,65839 | + 0,54230 ' — 0,07268 - 0,47419 0,71724° | — 0,82130 | — 0,80262 0,68227 
94 | + 1,73805° | + 0,63377° | + 0,02459' - 0,38399 | — 0,65011° | — 0,79279' 0,82448 | — 0,75978° 
0,95 + 1,83178 | + 0,74019 | + 0,13888 | — 0,27356 | — 0,55896 0,73697 0,81776 | — 0,81107 
96 | + 1,94591 | + 0,86807 | + 0,27707 | — 0,13540 0,43528' 0,64384 | — 0,77043° 0,82170° 
97 | + 2,09230° | + 1,02953° | + 0,45181 | + 0,04408 - 0,26403 | — 0,49627° | — 0,66250° | — 0,76807 


98 | + 2,29756° | + 1,25161° | + 0,69108 | +0,29437° | — 0,01348° | — 0,25926 | — 0,45458 0,6051 
99 | +2,64665 | + 1,62019° | + 1,08265 | + 0,70625° | + 0,41159 | + 0,16845 | + 0,03735 0,07590° 


1,00 | + + 00 + 4% +0 + +x + 


\ \ 1 


Abgeschlossen im August 1939. 156 








x y : r 2. . Math.Mech. 
280 Sauer, Über Interpolation von Kurvenscharen I eg 





Über Interpolation von Kurvenscharen mit Anwendung 
auf die Berechnung von Geschoßflugbahnen. 


Von R. Sauer in Aachen. 


ei Problemen der angewandten Analysis liegt häufig die Aufgabe vor, nicht nur eine ein- 
B zelne Integralkurve einer Differentialgleichung, sondern eine ganze Integralkurvenschar zu 
ermitteln. Man versucht dann aus möglichst wenig durchgerechneten Einzelkurven die übrigen 
Kurven der Schar durch Interpolation zu erhalten. Im folgenden wird gezeigt, wie durch 
gewisse affıne Transformationen die Interpolationsgenauigkeit gesteigert und dementsprechend 
die Anzahl der unmittelbar zu berechnenden Integralkurven verringert werden kann. Das Ver- 
fahren läßt sich in besonders einfacher Weise auf die Berechnung von Geschoßflugbahnen 
anwenden. Es spezialisiert sich hierbei zu der Interpolationsmethode, die kürzlich H. Athen‘) 
aus einem Ansatz von K. Popoff?) entwickelt hat. 


1. Voraussetzungen. Die zu interpolierende einparametrige Kurvenschar (Parameter A) 
denken wir uns in rechtwinkligen Koordinaten x, y mit Benützung einer Hilfsvariablen t dar- 
gestellt durch 
z=z{t, A), a 5 er 


Die Gl. (1) definieren neben den zu interpolierenden Kurven = const (kurz „A4-Kurven“) noch 
eine zweite Schar von Kurven t=const (kurz „t-Kurven“), welche auf den 4-Kurven nach t 
bezifferte Skalen ausschneiden (Bild 1; die -Kurven sind durchgezogen, die t-Kurven gestrichelt). 
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Wir setzen ein für allemal voraus: 

a) Das betrachtete x, y-Gebiet (@) ist einfach zusammenhängend und wird von 2 4-Kurven 
(A=4,4=4,) sowie von 2 t-Kurven (t=1t,,t=t,) begrenzt. Die Kurven {=1,=0 und t=t, 
dürfen auch zu Punkten einschrumpfen. 

b) Die 4-Kurven und t-Kurven bilden in dem Gebiet (@) ein Kurvennetz, d. h. durch 
jeden Innenpunkt von (@) geht genau eine 4-Kurve und eine t-Kurve. 

c) Die 4-Kurven haben stetig veränderliche Tangenten und Krümmungskreise. Sie sollen 
bei ?=0 mit nicht verschwindender Krümmung und nicht zur x-Achse senkrechter Tangente 
beginnen. Für die Anfangswerte von x, bei t=0 und ihre ersten und zweiten Ableitungen 
führen wir die Bezeichnungen ein 


2 (0,4)=a,(), y0,A)=b,(); 


d’x 


u ERRN .. et d’y . (2). 
Can W740, Frei} 4; (Tr 


E a,4), ( b, () 


ee 


Für den Tangentenwinkel # (4) und die Krümmung - bei {=0 hat man dann 
ER SR Ei 1 ,b5—aob 
= art tg a, 3% £ DR P == (a, +b,)® > 0. 


!) Z. angew. Math. Mech. Bd. 19 (1939), S 361 bis 371. 
2)C.Cranz: Ballistik, Ergänzungsband, Berlin 1936, S. 47 bis 49. 
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Für das Folgende ist es belanglos, ob die zu interpolierende Kurvenschar (1) durch 
Differentialgleichungen 1. Ordnung 


dx dy 
Fra) de v9) 


oder durch Differentialgleichungen höherer Ordnung vorgegeben ist und ob der Parameter A 
lediglich in den Anfangsbedingungen a, (4), b, (4) usf. oder auch in den Differentialgleichungen 
selbst explizit auftritt. Die Differentialgleichungen werden beim Interpolieren nur zur Be- 
rechnung der Anfangswerte a,,b, usf. verwendet. 


2. Affintransformation der zu interpolierenden Kurven. Vor der Interpolation wird das 
Gebiet (@) der xzy-Ebene durch eine Transformation (T) folgendermaßen in ein Gebiet (I) 
einer &,n-Ebene verzerrt (Bild 2): 


a) Jede 4-Kurve von (G) wird affin in eine 4-Kurve von (/') transformiert. 
b) Alle t-Skalen der 4-Kurven von (/') stimmen für =0 in den Anfangswerten 


/ 


dE 


a ERDE ; ns dn BI _.n. e2 d’n RR i 
AREA (ur EL, ar FT Ey (3) 


überein; die von A unabhängigen Konstanten a, 7:0, a,,#==0 können willkürlich gewählt 
werden. Die Bedingungen (3) bedeuten geometrisch, daß die t-Skalen aller 4-Kurven des Be- 
reiches (/') die 3 ersten „unendlich benachbarten“ Skalenpunkte gemeinsam haben. Außer- 
dem folgt aus (3): 

c) Alle 4-Kurven von (/') beginnen im Koordinatenanfangspunkt &=n=0, haben dort 
die &Achse als gemeinsame Tangente sowie die nämliche nicht verschwindende und endliche 


Krümmung 
A Y 
\ 7 
EA _ 
7. 
0 a, 
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Die affın verzerr- 
ten 4-Kurven erhält man 
durch folgende einfache 
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3 N Konstruktion (Bild 3 

i und Bild 4): 
5 re Die Punkte einer 
N 4Kurve seien in (@) 
nn durch die schiefwink- 

Sr . WG BA 8 Ginin). ligen Koordinaten 
—dq 

ud AR at y—b). » 2200. 


festgelegt. Diese Koordinaten s,g werden bei der affinen Deformation der 4-Kurve in 


o=ms, ng BEE ls ea 


verzerrt, während der Koordinatenwinkel 5— —® in w übergeht. Hierbei ist 


_ 


2 | a,ß 


a, > 2 22 2 . 
m? — 4 "—_ ee a +l(a,a,- a,a)),, Bo=- (6) 
a,? + b,?’ (a, b, - 6, d,); | ı P u aan „“ g ’ 


a,0,— a,4, 
wie sich aus (4) und (5) durch zweimalige Differentiation nach t bei Berücksichtigung der 
Anfangsbedingungen (2) und (3) sowie der Gleichungen 
E=06—x%C0Sw, n=xsin® 
nach kurzer Rechnung ergibt. 

Man beachte, daß m, n und ® Funktionen von A sind, die 4-Kurven also durch ver- 
schiedene Affintransformationen verzerrt werden. Die Transformation (T) des Gebietes (@) 
in das Gebiet (/’) ist infolgedessen selbst keine affine Transformation. 


3. Interpolationsvorschrift. Unter Beschränkung auf einen nicht zu großen t-Bereich 
werden durch die Transformation (T) des Gebietes (@) in das Gebiet (/') die 4-Kurven samt 
den ihnen aufgeprägten f-Skalen einander angenähert (vgl. Bild 1 und Bild 2) und die Bogen- 
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längen der #-Kurven verkürzt; denn auf Grund der Bedingungen (3) hat man bei t=0 für 
die Bogenlängen 7; (f) der t-Kurven des Bereiches (/') 
2 
Lo. d L_o. d L_ ) 
di d#? 

d. h. die 3 ersten „unendlich benachbarten“ f-Kurven von (/') schrumpfen auf Punkte ein. 
Infolgedessen liefert die Interpolation der 4-Kurven im deformierten Gebiet (/') genauere Er- 
gebnisse als im Gebiet (G@) vor der Deformation. Dies gibt Anlaß zu folgender Interpolations- 
vorschrift: 

a) Von der gesuchten 4-Kurvenschar des Gebietes (G) werden einige Kurven, etwa die 
Kurven A=4,,4=4,, A=4,, durch Integration der vorgelegten Differentialgleichungen un- 
mittelbar berechnet. 

b) Die Kurven Ä,, #,, 4, werden auf Grund der Gl. (5) von (@) in das Gebiet (/') affin 
transformiert. 

c) Zwischen entsprechenden Skalenpunkten der f-Skalen der Kurven 4,, 4,, 4, des Ge- 
bietes (/') wird linear oder quadratisch interpoliert, wodurch sich in (/') zwischengeschaltete 
4-Kurven samt t-Skalen ergeben. 

d) Die nach e) in (/’) zwischengeschalteten 4-Kurven werden gemäß (5) affın in das 
Gebiet (@) zurückübertragen. 

Sowohl die affinen Verzerrungen als auch die lineare oder quadratische Interpolation 
kann rechnerisch oder zeichnerisch durchgeführt werden. 

Falls die erzielte Genauigkeit noch nicht ausreicht, ist die Anzahl der berechneten 
4-Kurven in a) zu vergrößern oder der t-Bereich in entsprechend kleinere Teilbereiche zu 
zerlegen. 


4. Abschätzung der Genauigkeit. Bei der Abschätzung der Genauigkeit handelt es sich 
erstens um die bei der Interpolation der 4-Kurven in (/') auftretenden Fehler und zweitens 
um die Vergrößerung oder Verkleinerung dieser Fehler bei der affinen Rückübertragung der 
zwischengeschalteten 4-Kurven aus dem Gebiet (/') in das Gebiet (@). 

a) Die Interpolationsfehler 40, 4x in (/') werden näherungsweise gleich der Differenz 
der quadratisch interpolierten Werte o,,*, gegenüber den linear interpolierten Werten 07, x; 
gesetzt: 

do=o,— a, Axr=m, —xı-. 


b) Aus den Fehlern 40,4x in (/') ergeben sich mittels (5) bei der affinen Rücküber- 
tragung nach (@) die Fehler 


1 
As=—4o, Ag=—Ax. 
m n 


Als Längenfehler der interpolierten Punkte in (G@) hat man also 


+ RI RE Fee r 
f=Vidar +4 =) mA + dr a, 5inddoAx er 5° 


5. Anwendung auf die Interpolation von Flugbahnscharen. Besonders einfach gestaltet 
sich das Interpolationsverfahren bei der Anwendung auf die Berechnung von Luftschuß- 
tafeln. Hierbei sind für ein vorgegebenes Geschütz und Geschoß und ein empirisch be- 
stimmtes Luftwiderstandsgesetz bei normierten atmosphärischen Bedingungen die von einem 
festen Anfangspunkt <=y=(0 mit einer festen Anfangsgeschwindigkeit v, ausgehenden Ge- 
schoßflugbahnen zu bestimmen. Sie ergeben sich aus den Differentialgleichungen’) 


z=—cyvr&Klh,y, Yy=—-cy)vryKb,)W)-gaW) Wert) .. © 


für die Anfangsbedingungen z=y=0, v=»,, (Wi .) —=tg» bei t=0. Die Punkte bedeuten 
Differentiationen nach der Zeit #. i 

Als Parameter der Flugbahnschar nehmen wir den Abgangswinkel 9=4. Die 4-Kurven 
sind die gesuchten Flugbahnen, die t-Kurven sind die Kurven gleicher Geschoßflugzeit. Die 
Kurve {=0 entartet zum Anfangspunkt =y=0, das Gebiet (@) wird begrenzt von 2 Flug- 
bahnen 9=4,, 9=4, und einer Geschoßflugzeitkurve t=1, (Bild 5). 

Mit Benutzung der Differentialgleichungen (8) erhält man für die Anfangswerte (2) 


a,=b,—=0; a, =r, cos oo #5 ), b,=v,sind; a,=—c,v, K,cos®#, b,=— c,v, K,sin# — 9. 


3) ©. Cranz: Ballistik, Ergänzungsband, Berlin 1936, S. 58. 
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Durch geeignete Wahl der noch willkürlichen Anfangswerte (3), nämlich 


0,:=0,, en A a * 
spezialisieren sich die Gl. (6) zu 
n=n=1, w=nj2 
und die Affinitäten (5) zu 
EEE Re De 
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Bild 6 (rechts). 
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Die Konstruktionsvorschrift für die affinen Verzerrungen wird also denkbar einfach 
(Bild 5 und Bild 6): Die schiefwinkligen Koordinaten s,q der nichtverzerrten Flugbahnen 
in (@) sind gleich den rechtwinkligen Koordinaten o,x der verzerrten Flugbahnen in (/'). 


Aus (7) und (9) folgt 
=(4’ +49 —2AsAgsind<(As/+/dg’=s2[(dAs’+(Ag], 


d. h.: Bei der affinen Rückübertragung der in (/') durch lineare oder quadratische Interpolation 
zwischengeschalteten A-Kurven in das unverzerrte Schußtafelgebiet (@) ist die etwaige. Ver- 
größerung der Längenfehler stets kleiner als y 2. 


Das hier durch Spezialisierung von (6) hergeleitete Interpolationsverfahren deckt sich 
mit der eingangs erwähnten Methode!) von H. Athen. 

Nach Nr. 2 werden bei 
den affinen Deformationen (9) 
alle Flugbahnen in Kurven 
verzerrt, welche im Anfangs- 
punkt den nämlichen Krüm- 3 
mungskreis besitzen. Dies 
ergibt sich auch unmittelbar y 


durch folgende einfache geo- 
metrische Überlegung (Bild’7): 
02 


Für die Krümmung 4 nd 


>«w 
—3\ 





Js 





k-44 | 


7 

















= >X 
’ . Ae7) 2 
der nichtverzerrten Flug- : Bil 7 
. di. 
bahnen im Anfangspunkt hat 
man 
. 2Ah .  2Agcosd 
— lim 4 gest 


= — lim 
r, ss>0olds— Agqsnd’ > (As? ° 


Ne 1 
für die Krümmungen = der verzerrten Flugbahnen dagegen 
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Hieraus und aus dem Ansatz der Normalbeschleunigung im Anfangspunkt 
U 
= g,c0sd 
Y, 
folgt sofort 
1 _9 
9% % 

6. Verallgemeinerung und Anwendung auf ein weiteres ballistisches Problem. Das hier 
beschriebene Interpolationsverfahren geht davon aus, daß die Kurvenschar vor der Inter- 
polation in eine Schar von Kurven affin deformiert wird, welche 3 „unendlich benachbarte“ 
Punkte gemeinsam haben. Wenn von der Kurvenschar nicht nur Anfangswerte (Anfangs- 
punkte, Anfangstangenten usw.) sondern auch Endwerte bekannt sind, läßt sich das Verfahren 
verallgemeinern, derart, daß an Stelle der 3 „unendlich benachbarten“ gemeinsamen Punkte 2 
„unendlich benachbarte“ und ein dritter nicht benachbarter Punkt treten. 

Als Beispiel erwähnen wir folgendes ballistische Problem: 

Für ein vorgegebenes Geschütz und Geschoß und eine vorgegebene Anfangsgeschwindig- 
keit v, sind die Schußweiten & in Abhängigkeit vom Abgangswinkel #9 bekannt („Erdschuß- 
tafel“); gesucht sind die zugehörigen Flugbahnen („Luftschußtafel“). Da man jetzt für die 
Flugbahnen neben Anfangspunkt und Anfangstangenten auch die Endpunkte kennt, kann man 
die Flugbahnen vor der Interpolation affın so deformieren, daß sie alle mit gemeinsamer 
Tangente vom Anfangspunkt abgehen und in einem festen Endpunkt einmünden. Die so ver- 
zerrten Kurven werden unter Umständen noch weniger voneinander abweichen als bei den 
affinen Deformationen (5) und infolgedessen noch bessere Interpolationsergebnisse liefern. 178 


= const. 


Messung der Krümmung von flachen Kurven. 
Von Klaus Katterbach in Göttingen. 
(Aerodynamische Versuchsanstalt Göttingen, Institut für instationäre Vorgänge.) 


I. Einleitung. Die vorliegende Arbeit entstand im Rahmen von Versuchen, welche die 
Ermittlung der an bewegten Körpern wirksamen Kräfte zur Aufgabe hatten. Da die direkte 
Messung dieser Kräfte nicht möglich war, wurden die Körperbahnen photographisch festge- 
halten mit dem Ziel, ihre Bahnkrümmung in einer Reihe von Punkten zu bestimmen und 
daraus die Normalkräfte zu berechnen. Zur Bestimmung der Bahnkrümmung wurden zunächst 
die bekannten Meßverfahren herangezogen. Die Ergebnisse verschiedener Meßreihen wichen 
aber zu stark voneinander ab und konnten daher keineswegs befriedigen. Die Ursache dafür 
lag in der flachen Form der Versuchskurven. 

Fast alle bisher bekannten Verfahren beruhen nämlich auf der zweimaligen Differen- 
tiation und der Anwendung der Formel 





1 [23 t4 a. 
.,y .d+y®)®. 


Sie sind daher mit denselben Schwierigkeiten behaftet wie allgemein die graphische Differen- 
tiation. Diese sind vornehmlich einmal die genaue Konstruktion der Tangente, zum andern 
Interpolationsschwierigkeiten beim Zeichnen der Kurven y und y’, da die gegebenen Punkte 
von y selten und von y’ nie exakt sind. Bei Kurven mit großen Krümmungsradien kommt 
ein weiterer Umstand hinzu, der die Anwendung dieser Verfahren unmöglich macht. Für 
solche flache Kurven ist nämlich die Änderung der ersten Ableitung innerhalb eines nicht 


. [77 A ' 2 y= x D 
zu großen Intervalls gering und, da y’’ = =% p ist, der Winkel @ von der Größenordnung 


der absoluten Meßgenauigkeit. Es kann daher nicht mehr mit Sicherheit angegeben werden, 
inwieweit eine gemessene Krümmungsänderung den gegebenen Verhältnissen entspricht 
oder auf einen Meßfehler zurückzuführen ist. 

Die übrigen Krümmungsmeßverfahren, nach denen der Krümmungsradius durch Anlegen 
eines starren Lineals mit bekannter Krümmung ermittelt wird, fallen für genaue Messungen 
aus, da ihre Genauigkeit bestenfalls 10—20°/ beträgt. Ein besonderer Vorteil dieser „Anlege- 
verfahren“ ist es, daß sie im Gegensatz zu den Differentiationsverfahren bei ihrer Durch- 
führung unabhängig von der Wahl eines Koordinatensystems sind und hierin der Invarianz 
der Krümmung gegenüber orthogonalen Koordinatentransformationen entsprechen. Daher lag 
der Gedanke einer Weiterentwicklung dieser Verfahren nahe, bei der das starre Kurven- 
lineal durch ein biegsames ersetzt wurde, dessen Krümmung dann besonders zu be- 
stimmen war. 
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Die ersten Versuche befaßten sich mit einem Verfahren, das zwar nicht zum gewünschten 
Erfolg geführt werden konnte, in seiner Art aber sehr interessant ist und zudem erst den 
Weg zur endgültigen Lösung vorgezeichnet hat. Aus diesen Gründen sei es hier unter dem 
Namen des „elastomechanischen“ Verfahrens aufgeführt. Der Vorschlag dieses Verfahrens 
und seine erste Ausführung stammt von Dr. Heinz Billing (z. Zt. bei der Wehrmacht). 
Daraus wurde das „elastooptische“ Verfahren entwickelt, welches sich in der KONNEN 
und hinsichtlich der Meßgenauigkeit allen andern Verfahren als überlegen erwies 


II. Das elastomechanische Krümmungsmeßverfahren. Bei der Durchbiegung eines Balkens 
\ ii ' ! ; 1 A i ’ : 
gilt für den Krümmungsradius die Beziehung = en, worin r den Radius, M das Biege- 


moment, E den Elastizitätsmodul, / das Trägheitsmoment, also E: I die Biegesteifigkeit des 
Balkens bedeuten. Auf dieser Beziehung wurden die Versuche aufgebaut. 

Eine Blattfeder aus Stahl oder Tombak von etwa 0,5 mm Stärke ist in einem Rahmen 
aus Isolierstoff so aufgehängt, daß ihre Enden frei beweglich sind (Bild 1). Zu beiden Seiten 
der Feder sind eine Reihe Stellschrauben angebracht, mit denen der Feder die der Bahn- 
kurve entsprechende Krümmung gegeben wird. Zur Einstellung wird ein Diapositiv der zu 
messenden Kurve über die Feder gelegt. Durch Lenker an den Enden der Feder und Angriff 
der Stellschrauben in der Symmetrieebene ist dafür gesorgt, daß die Feder nicht tordiert wird. 

Zu messen sind gemäß obiger Beziehung bei bekannter Biegesteifigkeit die Momente 
bzw. die Biegekräfte, d. h. die Druckkräfte zwischen Schrauben und Feder. Zu diesem Zwecke 
wird über die Feder im Berührungspunkt einer Schraube ein Gehänge gelegt, das mit einer 
Waagschale verbunden ist. Die Waagschale wird so weit belastet, bis die Feder von der 
Schraube gerade abhebt. Zur genannten Feststellung des Moments der Abhebung sind Feder 
und Schraube in einen Stromkreis geringer Spannung gelegt, der beim Abheben unterbrochen 
wird. Zur besseren Isolation ist die Oberfläche geölt. 





Bild 2. Elastomechanisches Verfahren. Durchbiegung 
der Feder durch fünf Schrauben und Darstellung der 
Kräfte, die durch das Gewicht G zu überwinden sind. 
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Bild 1 (links). Aufhängung und Lenkerführung der Feder 
bei den elastischen Verfahren. 


Die Meßergebnisse mit dieser Anordnung sind für Kurven mit konstanter Krümmung 
befriedigend. Die Abweichungen betragen, wenn die Feder in jeder Beziehung möglichst 
gleichmäßig ist, rd. 5°. Bei andersartigen Kurven aber weichen die einzelnen Messungen 
voneinander stark ab. Zur Erklärung dieser Tatsache betrachten wir Bild 2. Der Feder ist 
durch eine Anzahl Schrauben die gewünschte Krümmung gegeben. Nach obigem sind die 
von den Schrauben aufgenommenen Druckkräfte zu messen, was durch Belastung der Waag- 
schale mit einem P, entsprechenden Gewicht @ geschieht. Ist @ genau P,, so wird zwar 
die Schraube entlastet, aber Feder und Schraube bleiben weiterhin in Kontakt. Zum Ab- 
heben der Feder ist ein Übergewicht @,„ erforderlich, das einmal die elastische Kraft der 
Feder, zum andern die Reibung an den übrigen Schraubenspitzen (dargestellt durch die kleinen 
Pfeile) zu überwinden hat. G,„ wird demnach größer 1. mit wachsender Krümmung und 
2. mit wachsender Anzahl der zur Einstellung erforderlichen Schrauben. Der Meßfehler ist 


aber gleich diesem Übergewicht, das von der Gesamtkraft nicht getrennt werden kann 
(6=P,+ 6). 
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Da die Kräfte von der Durchbiegung der Feder abhängig sind und zwischen zwei 
Federeinstellungen immer kleine Differenzen sein werden, ist auch die Erklärung gegeben, 
warum zwei Meßreihen schlecht übereinstimmen. 

Also auch dieses Verfahren befriedigt genau so wenig wie die anderen. Es würde 
vielleicht zum Erfolge führen, wenn man sich von dem Übergewicht freimachen könnte, 
indem man die Druckkräfte direkt mißt, etwa durch Deformation kleiner Meßfedern. Diese 
Möglichkeit ist nicht weiter untersucht worden, da das folgende Verfahren mehr Aussicht 


auf Erfolg hatte. 


III. Elastooptisches Verfahren. Bei diesem Verfahren wird eine Verbindung zwischen dem 
elastomechanischen und einem aus der praktischen Optik bekannten Krümmungsmeßverfahren 
für optische Flächen hergestellt. Wie oben wird einer elastischen Feder mittels Schrauben 
die gewünschte Krümmung gegeben. Die Feder ist jedoch in diesem Falle auf einer Seite 
mit einer optischen Ansprüchen genügenden Verspiegelung versehen. Gemäß Bild 3 fallen die 
von einem Objekt 0 ausgehenden Lichtstrahlen durch eine Zylinderlinse Z, auf die als Zylinder- 
spiegel wirkende Feder, werden hier reflektiert und ergeben nach nochmaligem Durchgang 
durch die Linse ein Bild 0’ von O, das zwar in Richtung der parallelen Achsen von Linse 
und Spiegel verzerrt, senkrecht dazu aber unverzerrt ist. Bei geeigneter Linsenstellung 
kommen Objekt und Bild zur Deckung, bzw. sie fallen in eine Ebene senkrecht zur optischen 
Achse der Anordnung. Mit den Bezeichnungen des Bildes 3 folgt dann: 

E r=xc+b 
ee | 
oder nach ET : fy 


u 
fy 
3 y=c+,—. 
Also +7 -Y 
Für konvexe Krümmungen folgt analog nach Bild 4: 
r=b—ı, 





Da f die Brennweite der Linse ist, sind also zur Bestimmung von r nur die beiden Abstände 
x und y zu messen'). 
























































Bild3. Elastooptisches Verfahren. Strahlen- | 
gang für konkave Krümmungen. v0 
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Bild 4. Elastooptisches Verfahren. Strahlen- 





gang für konvexe Krümmungen. “r — 
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Bild 5. Elastooptisches Verfahren. Schema u, 




















der Anordnung und Strahlengang. as] 
Bild 5. 


') Bei der Messung ist die Lage der Linsenhauptpunkte, die experimentell festgelegt werden müssen, zu be- 
riicksichtigen. 
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Die Anordnung des Apparate$ zeigt schematisch Bild 5. Ein Doppelspalt S wird mit 
einer intensiven Lichtquelle @ beleuchtet. Die Opalglasscheibe O zwischen Lichtquelle und 
Spalt gibt das gewünschte, gleichmäßige Streulicht. Das Spaltbild entsteht bei entsprechender 
Linsenstellung Z auf der Mattscheibe M und wird über ein Prisma P mit einem Mikroskop M; 
von etwa 10—12facher Vergrößerung betrachtet. Doppelspalt und Mattscheibe liegen in einer 
Ebene. Zur Messung wird durch eine Spaltblende B von der Feder F ein schmaler Bereich 
von etwa 2-3 mm ausgeblendet und senkrecht zur optischen Achse eingestellt. Zu diesem 
Zwecke ruht der Rahmen, in dem die Feder hängt, in einer Aufnahmevorrichtung, die sowohl 
gedreht als auch in der Längsrichtung der Feder verschoben werden kann, so daß jede ge- 
wünschte Einstellung der Spiegelfläche möglich ist. Sodann wird die Linse Z, die zweck- 
mäßig chromatisch korrigiert ist, so lange verschoben, bis das Bild auf der Mattscheibe die 
größte Schärfe hat. Schließlich folgt die Messung der Abstände. Die Genauigkeit, mit der 
sie gemessen werden müssen, ergibt sich aus der Bestimmungsgleichung für r. Danach ist 
für x eine Messung bis auf einige Millimeter vollkommen ausreichend. 4 dagegen muß recht 
sorgfältig gemessen werden, vornehmlich wenn flache Kurven vorliegen; denn für y gegen f 
geht der Ausdruck rY 

f-y 
die Abweichungen bei Radien von etwa 5 m 3°, bei Radien von 10 m unter denselben 
Voraussetzungen schon 10°. Zur Erreichung so genauer Einstellungen und Ablesungen ist 
daher die Verstelleinrichtung der Linse mit Skala und Trommelteilung versehen. 

Eine Schwierigkeit, auf die noch hinzuweisen wäre, liegt in der Einstellung der 
Krümmung der Feder gemäß der Vorlage. Auf sie sind vielfach die Abweichungen in den 
Ergebnissen zurückzuführen. Zur Einstellung benutzt man zweckmäßig eine schwach ver- 
größernde Aufsetzlupe mit Diopter. Damit erreicht man bei einiger Übung und entsprechender 
Sorgfalt gute Übereinstimmung. 

Manchmal erhält man Einzelwerte, die aus der Meßreihe erheblich herausfallen. Solche 
Abweichungen rühren von kleinen Buckeln oder Dellen der Feder her. Es ist darauf zu 
achten, daß die Feder in dieser Beziehung einwandfrei ist. Sollten trotzdem größere Ab- 
weichungen auftreten, so muß man zwischen den Nachbarwerten interpolieren. 

Die Einstellung des Federrandes auf die gegebene Kurve wird nie völlig genau sein 
können, schon aus dem Grunde, weil zwischen zwei Stellschrauben das Moment und daher 
nach der Balkentheorie auch die Krümmung eine lineare Funktion der Längenkoordinate sein 
muß. Das Auge des Einstellenden besorgt hier jedoch eine Ausgleichoperation gefühlsmäßig, 
die auch der Mathematiker nach der Methode der kleinsten Fehlerquadrate vornehmen würde, 
wenn ihm 15 Punkte der Kurve genau gegeben wären. ° Denn er würde sicher versuchen, 
diese 15 Punkte nicht durch eine ganze rationale Funktion 14. Grades, sondern durch eine 
einfachere Funktion niederen Grades darzustellen, welche an sich durch 15 Punkte über- 
bestimmt ist. 

Mitunter interessiert weniger die Krümmung als die beiden ersten Ableitungen, die bis 
heute nicht immer mit hinreichender Genauigkeit bestimmt werden können. Es ist Jaher 
vorgesehen, den Krümmungsmesser so auszubauen, daß gleichzeitig mit der Krümmung auch 
die erste Ableitung gemessen und dann die zweite Ableitung nach der Formel 


gegen Unendlich. Beträgt z. B. der Ablesefehler + 0,1 mm, so sind 


FR in? 
Mack a a Ei 


berechnet werden kann. Die diesbezüglichen Versuche sind zur Zeit noch nicht so weit durch- 
geführt, daß Aussagen über die Genauigkeit dieser Methode gemacht werden können. Sie 
hat aber vor den Differentiationsverfahren neben einer vielleicht größeren Genauigkeit auch 
den Vorteil, daß die Interpolation zwischen experimentell gegebenen Kurvenpunkten mit Hilfe 
des biegsamen Lineals verhältnismäßig schnell und genau ausgeführt werden kann. 


IV. Meßergebnisse und Auswertung. 


Zahlentafel 1. R 





Y T r 





Messungen an einem Kreis mit dem Radius r = 40 cm. 
7,6 | 22,2 | 39,8 
7,3 | 23.1 | 39,4 
7,2 | 254 | 41,2 
8.7 | 16.0 | 41,5 
82 | 17,5 | 39,1 
8.6 | 14,2 | 38,9 


79 | 20,3 | 40,0 Den Zahlen liegen drei Federeinstellungen zugrunde. Bei jeder 
79 | 20,3 | 40,0 Einstellung wurden drei Messungen an verschiedenen Punkten der 
78 | 20,4 | 39,5 Feder gemacht. 
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Zahlentafel 2. * 
Messungen an einer experimentell erhaltenen Kurve, Bild 6. (z£+%y)=31,3 cm. 

















%« | x r | Ir 








a ı 9 %» | U 





| 
0,0 1220 | 12,38 12,47 | 12,30 19,00 | 226,9 | 0,0044 
05 11,90 11,58 11,61 | 11,0 19,60 | 131,4 | 0,0076 
1,0 11.00 11,12 11,20 | 11,15 | 20,15 | 96,1 | 0,0104 


1,5 10,88 10,66 | 1085 | 1061 | 20,9 772 | 0,0189 
2,0 9,99 988 | 1040 | 10,06 2124 | 649 | 0,0154 
2,5 9,70 99 | 2m | 572 ' 0,0175 
3,0 _ 9 | — | 98 2197 | 546 | 0,0188 
3,5 9,38 932 | 952 9,35 2,95 | 548 | 0,0182 
40 9,52 SE 950 | 21,80 | 573 | 0,0174 
4,5 00 | — | 91 | 95 2155 | 60,0 | 0,0167 
5,0 26 9 — 9,97 21,33 | 634 | 0,0157 
5,5 10,33 10,15 10,15 10,25 21.05 | 685 | 0,0146 
6,0 10,42 10,20 10,40 10,58 20,77 | 7,1 | 0,0183 
65 | 1088 | 10,7 10,77 | 1080 | 2050 | 82,7 | 0,0121 
287 1: SH 11,00 | 11,10 | 2020 | 93,8 | 0,0107 
75.) 1161 | 1134 | 11,16 11,32 | 19,98 | 104,6 | 0,0096 
8,0 11,60 | 11,55 11,55 11,60 | 19,70 | 12287 | 0,0081 





85 11,98 | 11,62 11,80 11,87 | 1943 | 148,3 | 0,0067 
9,0 12,02 12,05 1220 1215 | 1,15 192,0 | 0,0051 
95 | 12,59 12,43 12,31 12,42 | 18,88 268,9 | 0,0037 
100 | 12374 | 18738 | 12369 | 12868 | 1862 | 443,7 | 0,0028 
10,5 12,81 12,87 | 18,00 12,90 18,40 | 1011,4 | 0,0010 





Zahlentafel 3. 
Vergleichsmessungen zu Zahlentafel 2 an der Kurve, Bild 6. x2-+%y=25 cm. 











ee A " | 9% REN RR, 1/r 
| | | | | | | 
00 | 1230| 22 12,36 | 12,64 | 240.1 | 0,0042 
05 | 1179 | 11900 | 1185 | 1815 | 1899 | 0,0071 
00 | ne | am | 1130 | 1370 | 97.1 | 00108 
15 | 10,70 | 10,89 | | 1082 | 14,18 77,1 0,0130 
2,0 10,21 | 10,53 1085| 14,75 | 61,8 0,0162 
2,5 10,02 | 10,00 | 997 | 2508 | 572 | 0015 
z 3,0 ‚1 | — | 980 | 1520 | 546 | 0,0183 
3,5 | 9,78 | 9,80 | 980 | 1520 | 54,6 0,0183 
4,0 9,95 | 9,98 | 992 15,08 | 56,2 0,0178 
4,5 10,10 | 10,12 ı 10,12 14,88 | 59,7 0,0167 
5,0 1023 | 10,34 10,34 14,66 | 64,1 0,0156 
5,5 10,49 | 10,50 | 10,57 14.43 | 69,6 | 0.0144 
6,0 10,76 10,84 | 10,80 1,20 | 764 | 0.0131 
6,5 11,00 | 11,18 | 11.08. | 1897 | 84,6 0,0118 
70 | 11,23 | 11,16 | 11,24 13,76 | 94,1 0,0106 
75 | 1141 | 11.88 ' 11,45 1355 | 105,9 | 0,0094 
80 | 11.70 | 11,80 | 11,65 13,35 | 1205 0,0083 





85 | 1195 | 11,90 11,90 13,10 | 145,1 | 0,0069 
.90 | 1220 | zur 12,16 12,84 187,6 | 0,0053 
95 | 1240 | 12,37 1241 | 1259 | 2586 | 0,0039 
100 | 271 | 12374 1268 | 1232 | 4373 | 0,0028 
10,5 | 12,8 | 123,92 ı 12,88 ' 1212 899,1 | 0,0011 

In den Zahlentafeln 2 und 3 sind mit a die Meßpunkte auf der Feder bezeichnet. Die y; 
bedeuten verschiedene Meßreihen, gewonnen aus ebenso vielen neuen Federeinstellungen. 
Die 4; sind als Funktionen von a graphisch dargestellt und von diesen Kurven die y, abge- 
lesen worden, die bei der Berechnung des Radius r eingesetzt sind.. Um vergleichbare Me&- 
reihen y; zu erhalten, ist bei den Messungen die Summe der Abstände («+,) konstant ge- 
halten worden. Sie beträgt für die Werte der Zahlentafel 2 31,3 cm, für die der Zahlentafel 3 
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Bild 6. Die Versuchskurve, 

die den Zahlentafeln 2 und 3 

und den graphischen Dar- 

stellungen der Bilder 7 und 8 
zugrunde liegt. 

Die Ku ist im Maßstab 1 2 9015|- 
wiedergegeben. 
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Bild 7. Darstellung der rezi- 0010 
proken Krümmungsradien der 
Versuchskurve Bild 6, ge- 
messen mit dem elastoop- 
tischen Krümmungsmesser. 
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Bild 8. Darstellung der rezi- 
proken Krümmungsradien 
der Versuchskurve Bild &, 
erhalten nach zweimaliger 
Differentiation mit dem 
Prismenderivator. f') 
an 
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Bild 7. Bild 8. 
25 cm. Alle Größen sind in cm gemessen. Beide Zahlentafeln beziehen sich auf dieselbe 
Kurve, Bild 6. Die Ergebnisse stimmen also gut überein. Die Abweichungen betragen im 
Mittel etwa 2°J. Stärkere Abweichungen treten nur bei den größeren Radien von einigen Metern 


auf. In Bild 7 ist z als Funktion von a graphisch dargestellt. 


Zum Vergleich ist die Krümmung der Kurve Bild 6 auch mit dem Prismenderivator bestimmt 
worden. Aus einer Anzahl von Messungen ergaben sich zwei Kurventypen für die Darstellung 


= —=f(a), die in Bild 8 durch die Kurven 1 und 2 wiedergegeben sind. Sie weichen unter- 


einander und auch von der Kurve Bild 7, die in Bild 8 als Kurve 3 nochmals eingezeichnet 
ist, teilweise erheblich ab. Die Ursache für diese schlechte Übereinstimmung ist, wie schon 


oben gesagt, in der flachen Ausgangskurve und in Verbindung damit in der Unzulänglichkeit 
Diese Tatsache sei an Hand des vorliegenden Beispiels 


des Primenderivators zu suchen. 
nochmals kurz erläutert. 

Die Neigungswinkel der Tangenten am Anfang und Ende der Kurve unterscheiden 
sich trotz stetiger Zunahme des Winkels längs der Kurve nur um etwa 6°. Das bedeutet 
bei etwa 20 Me&punkten im Mittel eine Winkeländerung von Punkt zu Punkt von etwa 0,3°. 
Die Ablesegenauigkeit des Prismenderivators beträgt höchstens + 0,2°, d.h. Winkeländerung 
und Meßgenauigkeit sind von derselben Größenordnung. Es ist also nicht einwandfrei fest- 
zustellen, inwieweit eine zwischen zwei Punkten gemessene Winkeldifferenz auf einen Meß- 
fehler oder eine Winkeländerung zurückgeht’). 

Infolge der Streuung von + 0,2° erhält man bei mehreren Messungen zur graphischen 
Darstellung von y’ Punkthaufen, durch die die Kurve y’ zu legen ist. Selbst bei Berück- 
sichtigung der Mittelwerte ist y’ niemals eindeutig. Stets wird man mit gleicher Berechtigung 
mehrere Kurven ziehen können, die zwar in den Einzelwerten wenig abweichen, deren Neigung 
aber beträchtlich verschieden sein kann. Die Folge davon sind große Abweichungen bei der 

zweiten Ableitung und der Krümmung. Die in Bild 8 dargestellten Kurven 1 und 2 sind aus 
denselben Meßreihen für die erste Ableitung hervorgegangen. Ihre Verschiedenheit ergibt sich 
aus zwei verschiedenen, in gleicher Weise möglichen Kurven y'. 

2) Man könnte die gegebene Kurve durch geeignete Zentralprojektion so verzerren, daß ihre Neigungen und 
Krümmungen groß werden, um diese Meßfehler herabzudrücken. Von einer praktischen Erprobung dieses Verfahrens 
ist nichts bekannt. 

20 
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Die größten Abweichungen von etwa 30°, bezogen auf den Höchstwert von Kurve 2, 
liegen im Gebiet stärkster Krümmung. Die Erklärung dafür gibt uns wieder die erste Ab- 
leitung, die an dieser Stelle einen Wendepunkt hat, für den die Tangentenrichtung infolge der 
geringen Meßgenauigkeit des Prismenderivators nicht sicher genug festgestellt werden kann. 
Aus dem Grunde sind in Bild 8 die beiden Kurven 1 und 2 in diesem Bereich gestrichelt, 
um die Unsicherheit der Messung anzudeuten. 


Anderseits wird die Messung mit dem elastooptischen Verfahren durch mehrere Tat- 
sachen erhärtet. 


1. Es sind im ganzen zehn unabhängige Federeinstellungen mit anschließender Messung 
gemacht worden, die alle, von geringen Abweichungen abgesehen, dieselben Werte 
wie die beiden aufgetragenen Meßreihen ergeben haben. 

fy 

ie . 
optischen Verfahrens um so größer, je größer die Differenz f—y wird, weil dann 
Ablesefehler nicht so stark ins Gewicht fallen. In dem strittigen Bereich ist f— y 
verhältnismäßig groß. 

3. Durch Anlegen von Kreisen an die Kurve wurde an der Stelle stärkster Krümmung 
ein Radius von etwa 60 cm gefunden. 





2. Gemäß der Bestimmungsgleichung: r=b+ ist die Genauigkeit des elasto- 


V. Zusammenfassung. Von den beiden hier zum erstenmal beschriebenen Krümmungsmeß- 
verfahren ist das elastooptische vor allem zur Messung an flachen Kurven mit gutem Erfolg 
anwendbar und den Differentiationsverfahren weit überlegen. Die einzige Schwierigkeit und 
Ursache größerer Fehler liegt in der Einstellung der Blattfeder gemäß der vorgegebenen 


Kurve. 


an Zeit ist nicht größer als bei der Differentiation. 
bauen, daß gleichzeitig mit der Krümmung auch die erste Ableitung erhalten wird. 


Bei mehrfacher Messung beträgt die erzielte Genauigkeit etwa 3°Jo. 


Der Aufwand 
Es ist vorgesehen, das Gerät so auszu- 
191 
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Bemerkung über den Einfluß des Seiten- 
verhältnisses auf den Zusammenhang 
von Normalkrait und Tangentialkraft 
eines Tragflügels. Mit dieser kleinen Be- 
merkung soll die Aufmerksamkeit auf eine einfache 
Tatsache der klassischen Tragflügeltheorie gelenkt 
werden, die vielleicht bei manchen Überschlags- 
rechnungen von Nutzen sein kann: Die Funktion 
Cn (c) ist praktisch vom Seitenverhältnis des 
Flügels unabhängig. 


Man hat frühzeitig bemerkt, daß die übliche 
Polare c„ (Cu) vom Seitenverhältnis!) A=F/b? 
des Flügels abhängig ist und unter der Voraus- 
setzung elliptischer Auftriebsverteilung die be- 
kannten Umrechnungsformeln entwickelt ?), die in 
den meisten praktisch vorkommenden Fällen aus- 
reichen, wie Messungen zeigen. Wir wollen nun 
beweisen, daß unter denselben Voraussetzungen in 
den Zusammerhang zwischen dem Normalkraft- 
beiwert c„ und dem Tangentialkraftbeiwert c; das 
Seitenverhältnis nicht eingeht. Zwischen dem Auf- 
trieb und dem induzierten Widerstand besteht die 
bekannte Beziehung c„;=Aca?/r, in welche das 
Seitenverhältnis als Faktor eingeht. Außerdem ist 
bei elliptischer Auftriebsverteilung 


ERERE 8 
Title 


wobei 
ee: (Se) 
Bra \da)i=o 


..» Wir rechnen mit A statt mit A=b2/F, weil sich so 
übersichtlichere Formeln ergeben. 


2) Siehe Ergebnisse der Aerodynamischen Versuchs- 
anstalt Göttingen, I. Lieferung, S 50 und Lehrbücher über 
Aeromechanik. 





gesetzt ist. (Bei der ebenen Platte ist o—=1, für 
Profile im allgemeinen etwas kleiner als 1.) Wir 
wollen nun in diese Beziehungen c„ und c; ein- 
führen. Unter den üblichen Beschränkungen auf 
kleine Anstellwinkel gilt: 


C9HA=Clat cv a 
G=—Ccaa+Cy- 


In der ersten Gleichung können wir c„« gegen Ca 
vernachlässigen, so daß sich ergibt: 


23 
n=a=17975° EEE ER 


Für c; erhalten wir mit „=cwi+t Cwp: 


2no 4n)o? 
NER de ER TE RE 11 rk x 
[7 i+27*® tar teen 


2no 
9 


= apart... . 


Eliminieren wir a aus (1) und (2), erhalten wir 
endgültig 


G=— — C?+cuwpr . (3), 
eine Beziehung, in welcher das Seitenverhältnis 
nicht mehr auftritt. 

Dieses Ergebnis können wir uns in einfacher 
Weise veranschaulichen. Wir denken an die 
Überlegung, welche zu den Umrechnungsformeln 
für die c„(c„)-Polare führt. Dabei wird unter 


der Voraussetzung elliptischer Auftriebsverteilung 
(a;=const. über die ganze Spannweite) für kon- 
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stantes c„ und damit für konstanten effektiven 
Anstellwinkel &,=ag— ca; die Änderung des 
geometrischen Anstellwinkels a, und die von Cwi 
mit dem Seitenverhältnis bestimmt. Da für zwei 
Vergleichsflügel von verschiedenem Seitenverhältnis 
der gleiche c„-Wert angenommen wird, so ist wegen 
des gleichen effektiven Anstellwinkels auch die 
Luftkraftresultierende nach Größe und Richtung 
die gleiche. Das heißt aber, daß auch c„ und c; 
sich nicht ändern. 
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Bild 1. 


Die bekannten Messungen aus der „I. Lieferung“, 
nach denen in Bild 1 c„ über c; für 7 Rechteckflügel 
verschiedenen Seitenverhältnisses aufgetragen ist, 
bestätigen die Richtigkeit der Beziehung (3) bis 
zum Seitenverhältnis 1:3 herunter. Der Gültig- 
keitsbereich unseres Satzes ist also offenbar der 
gleiche wie der der Umrechnungsformeln für die 
Ca (Cw)-Polare. 
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Göttingen. Dietrich Küchemann. 


teilung von Skalen anderer Funktionen. 


1. Die projektive Skala ist wegen ihrer großen 
Anpassungsfähigkeit außerordentlich brauchbar für 
die Unterteilung größerer Intervalle vorliegender 
Skalen, insbesondere wird man sie dazu beim Ent- 
wurf der Skalen solcher Funktionen verwenden, für 
die die Berechnung der Einzelwerte eine größere 
Rechenarbeit erfordert. Notwendig für eine der- 
artige Verwendung ist, daß man eine möglichst 
gute Fehlerabschätzung hat. Nun ist aber die 
übliche Fehlerabschätzung !) sehr wenig genau und 
führt dadurch zu einer Unterschätzung der Brauch- 
barkeit der Annäherung durch projektive Skalen. 
Zu einer besseren Abschätzung kommt man nach 
einer Vermutung, die zuerst Herr Tollmien aus- 
gesprochen hat, "mittels des Restgliedes der Thie- 
le schen Interpolationsformeln, die die reziproken 
diviıdierten Differenzen benutzen ?). Diese Form 
des Restgliedes läßt sich, wie im folgenden gezeigt 
werden soll, sehr einfach, ohne auf die Thiele- 
schen Formeln Bezug zu nehmen, ableiten; die 
numerische Auswertung der neuen Form des Rest- 


1) Z.B. na 
lin 1924, S. 


2) Z.B. y Milne- ngsss: The Caleulus of finit 
Differenees, London 1933, S. 


Lehrbuch der Nomographie. Ber- 


gliedes fordert allerdings etwas mehr Rechenarbeit, 
wie die der üblichen Form. 

2. Es seien die drei Konstanten a, b und d in 
- azx+b . ö Rp 
y= Fe so bestimmt, daß die Werte von % mit 
denen der gegebenen Funktion 4=/ (x) in den drei 
Punkten P, (2, %ı), Ps (2, Y,) und P, (25, %) über- 
einstimmen, d. h. diese Konstanten seien aus den 
drei linearen Gleichungen 


az +b—dy=zY 
a, +5 - u nn)... . 
az, +b—dy,= < 
berechnet. Setzt man dann 
az+b 
== + ee - 
(x) 435" R (x) (2) 


so wird das zur Abschätzung des Unterschiedes 
der Funktionswerte dienende Restglied R (x) Null 
werden, wenn x die Werte z,, z; oder z, annimmt, 
d. h. es wird den Faktor (2 — x,) (£ — 2) (e — z;) 
enthalten. Es erweist sich als praktisch für R (x) 
den Ansatz zu machen 


z— 2)(2 — 2) (r — r,) 


ER u & 
R (x) = (za) Pe 3 


Der Nenner x +.d ist für die hier in Frage kom- 
mende Verwendung stets von Null verschieden. 
Dann hat die Funktion 


= at+b (d—-z)(tt - 2)t—z,) 
eÜ)=F(t) {ta dLar x (2) 
die Eigenschaft, daß sie für die Werte t=z,,2», 
T gleich Null wird, denn für diese Werte sind die 
ersten beiden Summanden gleich und der letzte 
wird Null. Sie wird aber auch für =x Null, wie 
aus den Gl. (2) und (3) zu ersehen ist. Die gleiche 
Eigenschaft, vier Nullstellen in dem durch die vier 
Werte 2,, 25,2; und x bestimmten Bereich zu haben, 
hat auch die Funktion 


B(t)=(t+AM-p(Ü)=FE)(t + A)? 
—(t—z)(t - z,)(t 
Nach dem Rolleschen Satz wird also die dritte Ab- 


leitung ®”’ (£) eine Nullstelle £ in diesem Gebiet 
haben. Es wird also 


-(at+b)(t+d) 
2) 2 (8) 


2.(2) = 3; FAOTC + 4)?] 
sein. Damit ist das gewünschte Restglied ge- 
wonnen. 
(2-2) 2) (2 —2,) d® ,, BEN 
Ria)= Herta) Frau (4). 


Beträgt die Längeneinheit der Skala von /(z)I!cm. 

so hat man also mit diesem Wert von R (x) die 

Gleichung 

ac+b 
z+d 


und hierin ist Z2-R (x) so zu bestimmen, daß es die 
Zeichengenauigkeit nicht überschreitet, also etwa 


En ee 63 


Daraus läßt sich entweder das Intervall bestimmen, 
in dem eine Unterteilung mittels einer solchen pro- 
jektiven Skala erlaubt ist, oder es läßt sich danach 
die Zeicheneinheit Z! so wählen, daß in einem be- 
stimmten Intervall eine solche Unterteilung er- 
laubt ist. 


ID) =1- +88 .:,..® 
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3. Als Beispiel sei die Größe von / für eine Skala 
von y=tgzr bestimmt, in der man das Intervall 
zwischen 40° und 50° mittels projektiver Skala 
unterteilen kann. Für diese Skala erhält man aus 
den Gleichungen (1), falls man x in Grad mißt: 


— 2 — 12,156 
z — 102,156 ’ 


Si 


während sich nach der Formel von Schwerdt, 
der die Funktionen nebst erster und zweiter Ab- 
leitung für z=45° zur Übereinstimmung bringt 


— 2 — 12,295 
z — 102,295 


I | 


3 
ergibt. Der FehlernachSchwerdt R, =; FEN 


wo jetzt x in Bogenmaß zu messen ist, wird 





5) 16 -4cos?50° 


<[% -4cos’ 50° _ 50035. 
R<(ge 6 cost500 ws 


d. h. 1,-0,0025 < 0,01 cm oder /,=4 cm. 


Nach unserer Formel wird 
(x z,) (€ 2) (2 — zT) 
3! (2 — 1,785)? 


6-4c0L ,_ mgnp j 16sin? ._ . | 
a E-17B0° + (£-1,785) + 


R< 





cos? & 


cos?£) 


Die eckige Klammer ist negativ. Sie hat absolut 
genommen zwischen 40° und 50° ihren größten 
Wert bei 50° nämlich 6.5, während der vor der 
Klammer stehende Ausdruck etwa bei 48° also für 
ı=0,838 sein Maximum hat. Daraus ergibt sich 


R < 0,0003, 


womit sich für 230 cm berechnet. Die Werte 
der eckigen Klammer nähern sich mit abnehmen- 
dem Winkel zunächst der Null, werden dann positiv 
und zwar wird der größte Wert für {=0 ange- 
nommen; dieser ist etwas größer als 12. In dem 
ganzen Gebiet ergibt sich kein kleinerer Wert für 
l als 30 cm. Würde man dagegen noch das Stück 
von 50° bis 60° hinzunehmen, so würde hier der 
Fehler schon bei entsprechender Interpolation 
maximal auf über 0,5 mm wachsen; in diesem In- 
tervall müßte also eine Unterteilung in zwei Ab- 
schnitte vorgenommen werden, und jeder der 
beiden müßte mittels projektiver Skala weiter 
unterteilt werden. 


Rechnet man die Werte für tg z und % im Inter- 
vall 40° bis 50° aus und damit die genauen Werte 
für die Unterschiede, so ergibt sich aus diesen mit 
unserer Formel != 50 cm, mit der von Schwerdt 
dagegen /, =20 cm. Man sieht daraus, daß hier 
die Fehlerabschätzung mittels R wesentlich besser 
ist, als die mittels R,, abgesehen davon, daß die 
maximale Abweichung nach unserer Formel kleiner 
sein muß als die nach der Formel von Schwerdt. 


Dresden. Willers. 157 


BUCHBESPRECHUNGEN 


Dr. B. von JUHOS, Erkenntnisformenin 
Natur- und Geisteswissenschaften 
(Pan-Bücherei, Gruppe: Philosophie, Nr. 20). 57 8. 
Leipzig 1940, Rudolf Birnbach, Abt. Pan-Verlag. 
Preis brosch. 2,50 M. 


Die Schrift sucht die Frage nach Wesen und Ziel 
der natur- sowie der geisteswissenschaftlichen Er- 
kenntnisse zu beantworten; für die Naturwissen- 
schaften werden die Ausführungen an die physi- 
kalische, für die Geisteswissenschaften an die histo- 
rische Forschung angeschlossen, so daß die Ergeb- 
nisse der Überlegungen streng genommen zunächst 
nur für diese beiden Wissenschaften gelten. Verf. 
geht von der Physik aus und zeigt, daß man hier 
zwei Stufen von Gesetzen und dementsprechend 
auch zwei Stufen von Erkenntnissen zu unterschei- 
den hat. Wenn unter denselben Umständen auf ein 
bestimmtes Ereignis A immer das Ereignis B folgt, 
so spricht er von einem Gesetz erster Stufe. Eine 
Voraussage ist hier nur für eine ganz bestimmte, 
bereits beobachtete Ereignisfolge möglich. Kann 
man aber aus einer Reihe solcher Ereignisfolgen 
eine funktionelle Beziehung herauslesen, die dann 
auch Voraussagen von Ereignisfolgen ermöglicht, 
die nicht zur Gewinnung jener Beziehung benutzt 
worden sind, dann spricht er von Gesetzen — und 
Erkenntnissen — zweiter Stufe. Diese Darlegungen 
werden an Beispielen aus den beiden wichtigsten 
modernen physikalischen Theorien, der Quanten- 
und der Relativitätstheorie, erläutert. Die beiden 
für die Physik nachgewiesenen Erkenntnisstufen 
benutzt Verf. nun auch für die Geisteswissen- 
schaften. Das Problem der Willensfreiheit führt 
ihn zur Annahme auch von zwei Stufen des Deter- 
minismus und des Indeterminismus. Von hier aus 
gelangt Verf. zu dem Schluß, daß die Geschichts- 
wissenschaft nur Gesetze erster Stufe kennt, auf 
Gesetze zweiter Stufe und Voraussagen allgemeinen 
Charakters aber verzichten muß. Sie hätte das 


Ziel, durch Nacherleben einmaliger Fälle mensch- 
liches Erleben und Handeln zu verstehen und so 
durch Bewußtmachen wirksamer bedeutender seeli- 
scher Kräfte die Seele zu bereichern. Das Buch 
stellt einen interessanten Beitrag zu der Frage 
nach dem Unterschied zwischen natur- und geistes- 
wissenschaftlicher Forschung und Erkenntnis dar. 


Berlin. E. Mosch. 169 


Prof. WILHELM SPANNHAKE und Dipl.-Ing. 
HERMANN DOSENBACH, Geschwindigkeit 
und Druck in Schaufelkanälen hydro- 
dynamischer Getriebe. 1. Teil. (Deut- 
sche Kraftfahrtforschung, Heft 38.) 16 S. m. 26 Bil- 
dern. Berlin 1940, VDI-Verlag G.m.b.H. Preis 
brosch. 1,40 M. 


Über die in hydrodynamischen Getrieben tat- 
sächlich auftretenden Strömungsvorgänge sind bis- 
her nur wenige einwandfreie Untersuchungen durch- 
geführt worden. Deshalb ist es zu begrüßen, daß 
das Institut für Strömungsmaschinen an der Tech- 
nischen Hochschule Karlsruhe die hier vorliegende 
umfangreiche Aufgabe in Angriff genommen hat. 
Berichtet wird über Messungen an 3 geometrisch 
ähnlichen Holzmodellen verschiedener Größe von 
einem Leitapparat für Strömungsgetriebe. Als 
Strömungsmittel wurde Luft benutzt, weil damit 
die Messungen am leichtesten durchzuführen sind. 
Zur Messung wurden Staukugeln bekannter und 
neuer Bauart mit Erfolg verwandt. Das Ergebnis 
ist, daß die Geschwindigkeitsverteilung in dem 
untersuchten Bereich von der Reynoldsschen Zahl 
und von der relativen Rauhigkeit der Kanalober- 
fläche praktisch nicht beeinflußt wird. Die in den 
Schaufelkanälen auftretenden Verluste werden 


noch weiter untersucht und die Ergebnisse sollen 
in einer zweiten Arbeit veröffentlicht werden. 


Augsburg. 


E. Sörensen. 1% 
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Dr.-Ing. habil. RICHARD DZIALLAS, Unter- 
suchungen an einer Kreiselpumpe 
mit labiler Kennlinie. 62 S. m. 26 Bildern. 
4} 1940, VDI-Verlag G.m.b.H. Preis brosch. 
6,50 M. 


Es ist bekannt, daß Kreiselpumpen mit teilweise 
labiler Kennlinie bei bestimmten Betriebsverhält- 
nissen in den anschließenden Rohrleitungen starke 
Druckschwingungen hervorrufen können. Der Ver- 
fasser hat zum erstenmal die hier auftretenden 
Vorgänge eingehend rechnerisch und praktisch 
untersucht und dabei alle Einflußgrößen berück- 
sichtigt, die von Wichtigkeit sind. Die Rechnung 
führt zu einer nicht linearen Differentialgleichung 
„weiter Ordnung, für die Näherungslösungen ange- 
zeben werden, die in gewissen Bereichen mit den 
Versuchsergebnissen gut übereinstimmen. Für den 
Praktiker ist der durch Messung festgestellte Ver- 
lauf der Kennlinie für negative Fördermengen be- 
deutungsvoll. Ebenso wichtig ist aber der Hinweis 
darauf, daß die bei stationärer Strömung gemessene 
Kennlinie nicht mehr gilt, wenn die Fördermenge 
pulsiert. — Die Versuchsergebnisse können nicht 
den Anspruch erheben, allgemeine Gültigkeit zu 
besitzen. Sie sind aber dennoch sehr wertvoll, 
weil sie einen Einblick in bisher wenig bekannte 
Vorgänge vermitteln. Für die Untersuchung von 
konkreten Fällen empfiehlt sich die zeichnerische 
Behandlung der Schwingungsvorgänge, auf die der 
Verfasser verschiedentlich hinweist, wenn er sie 
auch selbst nicht angewandt hat. 


Augsburg. E. Sörensen. 200 


Dr.-Ing. HEINRICH PRESS, Der Boden als 
Baugrund mit Ergebnissen eigener 
Versuche. 2. neubearb. Aufl. (Mitteilungen 
aus dem Gebiete des Wasserbaues und der Bau- 
erundforschung, Heft 9). V-+49 S. m. 64 Text- 
abb. Berlin 1940, Verlag Wilhelm Ernst & Sohn. 
Preis brosch. 5,80 M. 


Das Heft gibt auf 49 Seiten eine kurze Darstel- 
lung der neueren Baugrundlehre. In dem Be- 
streben, eine gewisse Vollständigkeit zu erreichen, 
mußte die telegrammstilartige Darstellung der 
ersten Auflage, auf deren Besprechung in dieser 
Zeitschrift 1940, S. 126, hingewiesen werden darf, 
beibehalten werden, die dem Anfänger, Nichtfach- 
mann oder Grenzgebietler die Benutzung nicht ge- 
rade erleichtern wird. Durch die Fülle des aus- 
zugsweise angeführten fremden Stoffes werden die 
eigenen Versuche des Verfassers leider etwas er- 
drückt. Preß hat eine große Anzahl von Ver- 
suchen im kleinen und großen Maßstab, sowie 
Messungen am Bauwerk durchgeführt und verfügt 
über zahlreiche Erfahrungen und Beobachtungen, 
die es hätten gerechtfertigt erscheinen lassen, 
wenn er sich auf die Darstellung seiner Arbeiten 
und die unmittelbar benachbarten Arbeiten anderer 
Forscher beschränkt hätte. Der Wert des Heftes 
hätte u. E. dadurch nur gewinnen können. 

Die neue Auflage beseitigt hinsichtlich der An- 
führung fremder Arbeiten die Mängel, die bei der 
Besprechung der ersten Auflage ausgestellt werden 
mußten. Es ist zu begrüßen, daß der Verfasser auf 
jene Anregungen weitmöglichst eingegangen ist. 

Naumburg (Saale). Dr. Scheidig. 194 


Dr.-Ing. E. h. WOLMAR FELLENIUS, Prof. a. d. 
Kgl. Techn. Hochschule zu Stockholm, Erdsta- 
tische Berechnungen mit Reibung 
und Kohäsion (Adhäsion) und unter 
Annahme kreiszylindrischer Gleit- 
flächen. 2, ergänzte Aufl. IV+48 S. m. 
38 Textabb. Berlin 1940, Verlag Wilhelm Ernst & 
Sohn. Preis brosch. 4,80 M. 


Seit dem Erscheinen der 1. Auflage im Jahre 1926 
gehört das Büchlein von Fellenius zum unentbehr- 
lichen Rüstzeug jedes Erd- und Tiefbauingenieurs, 
der mit erdstatischen Berechnungen für Böschun- 
gen großer Dämme, Rutschungsfragen, Kaibauten 
u. dgl. zu tun hat. 

Die Annahme kreiszylindrischer Gleitflächen bei 
erdstatischen Gleichgewichtsberechnungen, die sich 
schon seit langem im Schrifttum findet, ist durch 
die systematische Behandlung von Fellenius prak- 
tisch erst fruchtbar geworden. Die in der Praxis 
meist vorkommenden verwickelten Fälle gestatten 
keine rein mathematischen Lösungen und die un- 
günstigste Lage der Gleitfläche muß durch Pro- 
bieren gefunden werden. Hierfür gibt die Arbeit 
wertvolle Fingerzeige durch systematische Vor- 
arbeit. Wenn auch neuerdings an Stelle des Gleit- 
kreises die logarithmische Spirale vorgeschlagen 
worden ist, werden die Methoden von Fellenius 
doch ihre Bedeutung nicht verlieren. Dasselbe gilt 
von der Annahme über das Zusammenwirken von 
Reibung und Kohäsion, die trotz der Erkenntnisse 
der neuen Erdbauforschung noch nicht durch eine 
praktisch brauchbare neue Grundlage ersetzt wer- 
den konnte. 

Durch das Fehlen der Arbeit wäre nach dem Ab- 
satz der ersten Auflage eine empfindliche Lücke 
entstanden. So ist es zu begrüßen, daß sich der 
Verlag zu einer Neuauflage in ergänzter Form ent- 
schlossen hat. 


Naumburg (Saale). Dr. Scheidig. 180 


HANS ERTEL, Elemente der Operatoren- 
rechnung mit geophysikalischen An- 
wendungen. 133 S. m. 8 Textabb. Berlin 1940, 
Verlag Julius Springer. Preis 7,50 M. 


Die kleine Schrift von Ertel gliedert sich in vier 
Hauptabschnitte. Der erste (20 S.) gibt allgemeine 
mathematische Ausführungen über Differential- 
gleichungen, enthaltend die Definitionen und einige 
wichtige Sätze, meistens ohne Beweis, aber an Bei- 
spielen erläutert. Der zweite Teil (39 S.) bringt die 
Differentialgleichungen der mathematischen Physik. 
Unter den gewöhnlichen D.-Gl. finden wir die Be- 
wegungsgleichungen der Punktdynamik, bezogen 
auf ein Inertialsystem oder auf ein rotierendes 
Koordinatensystem, die  Bewegungsgleichungen 
eines Elektrons im elektromagnetischen Feld und 
die Schwingungsgleichungen. Von den partiellen 
D.-Gl. werden die Potentialgleichungen, die Glei- 
chungen der Hydrodynamik, der Elektrodynamik, 
der Wärme und Diffusion behandelt. Auf dem vor- 
gesehenen knappen Raum ist eine ausführliche Dar- 
stellung und Herleitung freilich nicht möglich; es 
handelt sich eher um ein kurzes Repetitorium aus 
der mathematischen Physik. 

Der dritte Teil (37 S.) behandelt das eigentliche 
Thema, die Elemente der Operatorenrechnung. Hier 
werden die Hauptregeln für das Rechnen mit dem 
Differentialoperator formal hergeleitet; für eine 
eigentliche Begründung und Darlegung der Gültig- 
keitsgrenzen fehlt natürlich der Raum. 

Im vierten Teil (31 S.) werden folgende geophysi- 
kalische Aufgaben mittels Operatorenrechnung ge- 
löst: Die Elektronenbewegung und die Ausbreitung 
der elektrischen Wellen in der Ionosphäre; harmo- 
nische Wasserwellen; die Eigenschwingungen von 
Seen; die durch den Wind hervorgerufenen statio- 
nären und nichtstationären Driftströme im Ozean; 
der Ausgleich von Salzgehaltsstörungen im Ozean 
durch Turbulenz, und die Wärmeleitung im Erd- 
boden. Ein Verzeichnis des Schrifttums, eine 
ÖOperatorentabele, Namen- und Sachregister 
schließen das Buch ab. 

Es bringt in gedrängter Kürze einen umfang- 
reichen Stoff und stellt an die Konzentration des 
Lernenden keine geringen Ansprüche. Der fort- 
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aber die Fehlerabschätzung war in dieser Unter- 
suchung bis zu einem gewissen Grade durch den 
günstig ausgefallenen Vergleich mit dem Experi- 
ment ersetzt. Durch meine gasdynamische Arbeit 
(diese Zschr., Bd. 17 (1937), S. 117 bis 136), auf die 
Herr Ringleb in diesem Zusammenhang über Ge- 
bühr hinweist, ist die Ungefährlichkeit des Durch- 
gangs durch die Schallgeschwindigkeit für adiaba- 
tische Potentialströmungen allenfalls an Hand von 
exakten Beispielen nur noch bekräftigt worden. 


2. Dagegen sind zum ersten Male in der er- 
wähnten Arbeit von mir Grenzlinien adiaba- 
tischer Potentialströmungen ausführlich diskutiert 
worden, über die hinaus eine Fortsetzung von Strö- 
mungen der angegebenen Art nicht möglich er- 
scheint. Herr Ringleb bezeichnet diese Grenzlinien 
in seiner Arbeit als Stoßlinien, was vielleicht 
zu Verwechselungen mit Verdichtungsstößen führen 
kann, die zunächst etwas anderes sind. Die Strö- 
nung hat an diesen Grenzlinien gewisse singuläre 
Eigenschaften, und auf die singulären Eigen- 
schaften allein kommt es an. Diese sind, so wie 
sie Herr Ringleb formuliert, bereits vollständig an 
den in meiner Arbeit behandelten Strömungen von 
mir seinerzeit festgestellt. Das Druckgefälle bzw. 
die Beschleunigung senkrecht zur Grenzlinie ist un- 
endlich, ebenso die Krümmung der Stromlinien 
oder anders ausgedrückt, die Stromlinien haben 
dort eine Rückkehrkante. Weiter hat Herr Ringleb 
und vor ihm schon W. G. Bickley *) bemerkt, daß 
die Grenzlinie eine Machsche Welle ist; denn die 
Geschwindigkeitskomponente senkrecht zur Grenz- 
linie ist gleich der Schallgeschwindigkeit. Dies 
letztere bedeutet freilich kein singuläres Verhalten. 
Ich hatte aber seinerzeit über die von Herrn Ring- 
leb wiedergefundenen Eigenschaften hinaus die 
m. E. grundlegende Tatsache festgestellt, daß die 
Grenzlinie eine Enveloppe Machscher Wellen ist. 
Da Herr Ringleb nicht auf meine Untersuchungen 
zur Frage der Grenzlinien eingeht, scheint mir diese 
Erinnerung angebracht. 

Es ist sehr begrüßenswert, daß Herr Ringleb zu 
dem interessanten und wichtigen Phänomen einer 
natürlichen Grenze adiabatischer Potentialströmun- 
gen neue Beispiele beigesteuert hat. Besonders 
erwünscht wäre der Nachweis, daß auch bei den 
Strömungen von Herrn Ringleb die Grenzlinien 
Enveloppen Machscher Wellen sind. 

Dresden. W. Tollmien. 213 





1) Vgl. G. I. Taylor: Recent work on the flow of ecom- 
pressible fluids. Journal of London Mathem. Soe. Bd. 5 
(1930), S. 224, Deutsche Übersetzung in Z. angew, Math. 
Mech., Bd. 10 (1930), S. 334. Taylor hat der Frage einer 
natürlichen Grenze adiabatischer Potentialströmungen viel 
Nachdenken gewidmet, ohne jedoch zu einem klaren Ergeb- 
nis zu kommen. 


Erwiderung. Herr Tollmien bemängelt in den 
vorstehenden Ausführungen, daßichinmeiner 
Arbeit nicht genügend auf sein Bei- 
spiel von Grenzlinien hingewiesen 
habe, welche sich begrifflich mit den Stoß- 
linien in meinen Beispielen decken und die 
gleichen Eigenschaften zeigen. 

Das Beispiel von Herrn Tollmien hielt ich für so 
allgemein bekannt, daß es mir zu genügen schien, 
wenn ich die Tollmiensche Arbeit gleich am An- 
fange meines Aufsatzes in der üblichen Weise 
zitierte. Ich hätte allerdings wohl besonders her- 
vorheben sollen, warum ich die Bezeichnung 
„Stoßlinie“ wählte. Das möchte ich hier nach- 
bolen und überhaupt in diesem Zusammenhange 
einige Feststellungen von allgemei- 
nerem Interesse treffen: 


In meiner Arbeit handelt es sich viel weniger 
darum, einzelne neue Beispiele für adiaba- 
tische Strömungen mit Stoßlinien (Grenzlinien) an- 
zugeben, als vielmehr darum, ein weitreichen- 
des Integrationsverfahren zur Lösung 
der Differentialgleichungen einer adiabatischen 
Gasströmung zu entwickeln. Ich möchte daher zu- 
nächst besonders hervorheben, daß es sich bei 
meinen Lösungen nicht um Einzelbeispiele oder 
spezielle Lösungsscharen — wie bei Herrn Toll- 
mien—handelt, sondern um eine Lösungs- 
gesamtheit, welche die gleiche Mäch- 
tigkeit besitzt wie die Gesamtheit 
allerinkompressiblen Potentialströ- 
mungen. 

Die Beispiele meiner Arbeit lassen sich dement- 
sprechend beliebig vermehren. Worauf ich nun 
durch die kurzen Angaben in der Einleitung meiner 
Arbeit hinweisen wollte, war die Tatsache, daß 
die Stoßlinie —- die also Herr Tollmien erst- 
mals in einem Spezialfall erhielt — eine ganz 
allgemein auftretende Erscheinung 
ist. Das zeigten zunächst noch zahlreiche andere 
aus dem allgemeinen Verfahren sich ergebende 
Beispiele, die ich nicht in meine Arbeit aufge 
nommen habe. Ich habe überdies in einer weiteren 
Arbeit, die seit längerer Zeit im Druck ist und in 
der „Deutschen Mathematik“ erscheinen wird, die 
Theorie der Stoßlinie als natürliche Grenze einer 
adiabatischen Gasströmung in voller Allgemeinheit 
entwickelt. Auf diese Arbeit habe ich übrigens auf 
S. 190 in einer Fußnote hingewiesen, die Herr 
Tollmien vielleicht übersehen hat. Die Taylorsche 
Frage nach den natürlichen Grenzen einer 
adiabatischen Strömung und die Anregung von 
Herrn Tollmien betreffend die Stoßlinien als 
Einhüllende der Machschen Linien haben 
in dieser meiner zweiten Arbeit, auf die ich hier noch- 
mals verweise, bereits ihre volle Klärung erfahren. 

Ich möchte nun noch begründen, warum ich die 
Bezeichnung Strömungsstoß bzw.Stoßlinie 
gewählt habe und nicht Grenzlinie oder sin- 
guläre Linie. Die letzte Bezeichnung ist nicht 
zu empfehlen, weil die Bezeichnung singulär schon 
allzuviel benutzt wird und hier zu mancherlei Irr- 
tümern führen könnte, z. B. zu einer Verwechse- 
lung mit Strömungssingularitäten. Gegen die Be- 
zeichnung „Grenzlinie‘ kann man einwenden, daß 
es fraglich ist, ob es sich dabei wirklich um die 
Grenze einer realen Strömung handelt. Es ist nur 
bewiesen, daß eine analytische Fortsetzung 
einer vorliegenden Lösung über ihre „Grenzlinie“ 
nicht möglich ist. Demgegenüber ist die physika- 
lische Bezeichnung Stoßlinie vollkommen ein- 
wandfrei, da eben das Auftreten eines Stoßes eine 
die betreffende Linie eindeutig charakterisierende 
physikalische Eigenschaft ist. Eine 
andere Bezeichnungsmöglichkeit wäre noch „Stau- 
linie“, weil die allgemeine Bedingung für eine 
solche (S. 190 Fußnote) im Unterschallgebiet die 
Staupunkte liefert, die Staulinie also als deren Ver- 
allgemeinerung im Überschallgebiet aufgefaßt wer- 
den kann. In einem Staupunkt tritt aber kein Stoß 
ein, so daß doch ein wesentlicher Unterschied vor- 
liegt. Es bleibt also als treffendste Bezeichnung 
die der „Stoßlinie“. 

Schließlich möchte ich noch bemerken, daß ich 
es nur deshalb für angebracht hielt, nochmals aut 
das im allgemeinen ungefährliche Verhalten einer 
adiabatischen Strömung bei Durchgang durch die 
Schallgeschwindigkeit hinzuweisen, weil trotz der 
Existenz einer Lavaldüse in dieser Hinsicht öfters 
Unklarheiten auftauchen. 
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